Background: This is the analysis script for Study 1 of the 2025 Cogsci Paper: “Adults hold two parallel causal frameworks for reasoning about people’s minds, actions and bodies” by Joseph Outa and Shari Liu.

Study 1 presented two tasks in a random order in a within-subjects design: a Freesorting Task and a Causal Task. Subjects sorted 15 items in a sorting canvas, and rated causal relatedness between all 210 pairs of the 15 items (excluding self pairs).

The N was 102 adults (before exclusions), recruited online via Cloud Research Connect. Study demo: https://liulaboratory-experiments.github.io/mbc/study1

The two research questions are: (1) How do subjects sort items? Do they sort based on three latent causes, 2 latent causes, 6 latent causes, or feature-based similarity. (2) Do subjects’ direct cause representations (from the Causal Task) reflect latent cause boundaries or do they go beyond, and if so, what is their structure?

The manuscript and project materials for the pilot study can be found at: https://github.com/liulaboratory-experiments/mind-body-action-causation

Script description: This script:-
- proceeds from the output of the preprocessing script.
- takes as input data_processed_step1.csv, and data_model_prediction.csv - plots Sorting and Causal Task rdms, gets correlations between empirical and models RDM, gets within subject correlations.
- outputs data_freesort_mins_causals.csv for use in study 2 script for item selection.

1 Setup

1.1 load packages

1.2 set variables

# variables
desired_order <- c("see something", "hear something", "choose what to do", 
                   "remember something", "think about something", "reach for something", 
                   "sit down", "jump up and down", "kick something", "take a walk", 
                   "get tired", "become hungry", "feel scared", "experience pain", "get sick")
domain_order <- c("mind", "action", "bio")
freesort_columns <- c("subject_id", "x", "y", "item")
freesort_columns_with_trial_order <- c("subject_id", "x", "y", "item", "trial_type", "trial_type_order")
causation_columns <- c("subject_id", "itemA", "itemB", "response")

# paths
## script inputs
base_path <- here("code", "pilot", "analysis")
snapshots_input_path <- paste0(base_path, "/study1-measure/snapshots_study1_preprocessing_analysis") 
data_input_path <- paste0(snapshots_input_path, "/data_processed_step1.csv")
model_predictions_input_path <- paste0(snapshots_input_path, "/data_model_predictions.csv") 
rdms_input_path <- paste0(base_path, "/study1-measure/figures/theoretical-rdms")

## script outputs
snapshots_output_path <- paste0(base_path, "/study1-measure/snapshots_study1_pilot_analysis")
rdms_outputs_path <- paste0(base_path, "/study1-measure/figures")
freesort_minus_causal_distances_output_path <- paste0(snapshots_output_path, "/freesort_minus_causals.csv")

## saved dataframes (not to be used in subsequent scripts, but saving for efficiency)
data_processed_output_path <- paste0(snapshots_output_path, "/data_processed_step2.csv") #not used in a subsequent script
correlations_output_path <- paste0(snapshots_output_path, "/correlations.csv") #not used in this script

1.3 read data

d <- read_csv(data_input_path)
## Rows: 12200 Columns: 19
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (8): subject_id, trial_type, trial_type_order, itemA, itemB, categoryA, ...
## dbl (9): trial_index, question_index, response, correct_response, original_t...
## lgl (2): attention_check, correct
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
#how many subjects
d %>% 
  slice(length(d$subject_id)) %>% 
  pull(subject_id)
## [1] "subj_50"

1.4 reshape data

Create separate data frames for Sorting Task, Causal Task and a combined data frame

df_causal <- d %>%
  filter(trial_type == "causation", !attention_check, itemA != itemB) %>%
  select(all_of(causation_columns)) %>%
  
  #normalize causal distances and flip scale
  mutate(
    causal_distance = (1 - response / 100 - min(1 - response / 100, na.rm = TRUE)) /
                            (max(1 - response / 100, na.rm = TRUE) - min(1 - response / 100, na.rm = TRUE))
  ) %>%
  select(subject_id, itemA, itemB, causal_distance) %>% 
  mutate(domain_itemA = case_when(itemA == "hear something" ~ "mind", #labels for first item of pair
                             itemA == "choose what to do" ~ "mind",
                             itemA == "remember something" ~ "mind",
                             itemA == "think about something" ~ "mind",
                             itemA == "see something" ~ "mind",
                             itemA == "reach for something" ~ "action",
                             itemA == "sit down" ~ "action", 
                             itemA == "jump up and down" ~ "action",
                             itemA == "kick something" ~ "action",
                             itemA == "take a walk" ~ "action",
                             itemA == "get tired" ~ "bio", 
                             itemA == "become hungry" ~ "bio",
                             itemA == "feel scared"~ "bio",
                             itemA == "experience pain"~ "bio",
                             itemA == "get sick"~ "bio"),
         domain_itemB = case_when(itemB == "hear something" ~ "mind", 
                               itemB == "choose what to do" ~ "mind",
                               itemB == "remember something" ~ "mind",
                               itemB == "think about something" ~ "mind",
                               itemB == "see something" ~ "mind",
                               itemB == "reach for something" ~ "action",
                               itemB == "sit down" ~ "action", 
                               itemB == "jump up and down" ~ "action",
                               itemB == "kick something" ~ "action",
                               itemB == "take a walk" ~ "action",
                               itemB == "get tired" ~ "bio", 
                               itemB == "become hungry" ~ "bio",
                               itemB == "feel scared"~ "bio",
                               itemB == "experience pain"~ "bio",
                               itemB == "get sick"~ "bio")) %>% 
  relocate(domain_itemA, .after = "itemB") %>% 
  relocate(domain_itemB, .after = "domain_itemA") %>% 
  mutate(itemA = factor(itemA, levels = desired_order),
         itemB = factor(itemB, levels = desired_order)) 


df_freesort <- d %>%
  filter(trial_type == "free_sort") %>%
  select(all_of(freesort_columns_with_trial_order)) %>%
  rename(item_freesort = item) %>%
  
  #expand pairs within each subject
  group_by(subject_id, trial_type_order) %>%
  summarize(pairs = list(expand.grid(itemA = item_freesort, itemB = item_freesort, stringsAsFactors = FALSE)), .groups = "drop") %>%
  unnest(pairs) %>%
  filter(itemA != itemB) %>%
  
  #join in freesort coordinates
  left_join(d %>% filter(trial_type == "free_sort") %>%
              select(subject_id, itemA = item, x, y),
            by = c("subject_id", "itemA")) %>%
  rename(x_itemA = x, y_itemA = y) %>%
  
  left_join(d %>% filter(trial_type == "free_sort") %>%
              select(subject_id, itemB = item, x, y),
            by = c("subject_id", "itemB")) %>%
  rename(x_itemB = x, y_itemB = y) %>%
  
  group_by(subject_id) %>% 
  
  #compute and normalize Euclidean distances
  mutate(
    euclidean_distance = sqrt((x_itemA - x_itemB)^2 + (y_itemA - y_itemB)^2),
    freesort_distance = round( (euclidean_distance - min(euclidean_distance, na.rm = TRUE)) /
                    (max(euclidean_distance, na.rm = TRUE) - min(euclidean_distance, na.rm = TRUE)),2),
    
    trial_type_order = ifelse(trial_type_order == "first", "causation_second", "causation_first"),
    
    itemA = factor(itemA, levels = desired_order),
    itemB = factor(itemB, levels = desired_order)) %>%
  select(subject_id, itemA, itemB, freesort_distance, trial_type_order) %>%
  ungroup() %>% 
  mutate(domain_itemA = case_when(itemA == "hear something" ~ "mind", #labels for first item of pair
                             itemA == "choose what to do" ~ "mind",
                             itemA == "remember something" ~ "mind",
                             itemA == "think about something" ~ "mind",
                             itemA == "see something" ~ "mind",
                             itemA == "reach for something" ~ "action",
                             itemA == "sit down" ~ "action", 
                             itemA == "jump up and down" ~ "action",
                             itemA == "kick something" ~ "action",
                             itemA == "take a walk" ~ "action",
                             itemA == "get tired" ~ "bio", 
                             itemA == "become hungry" ~ "bio",
                             itemA == "feel scared"~ "bio",
                             itemA == "experience pain"~ "bio",
                             itemA == "get sick"~ "bio"),
         domain_itemB = case_when(itemB == "hear something" ~ "mind", 
                               itemB == "choose what to do" ~ "mind",
                               itemB == "remember something" ~ "mind",
                               itemB == "think about something" ~ "mind",
                               itemB == "see something" ~ "mind",
                               itemB == "reach for something" ~ "action",
                               itemB == "sit down" ~ "action", 
                               itemB == "jump up and down" ~ "action",
                               itemB == "kick something" ~ "action",
                               itemB == "take a walk" ~ "action",
                               itemB == "get tired" ~ "bio", 
                               itemB == "become hungry" ~ "bio",
                               itemB == "feel scared"~ "bio",
                               itemB == "experience pain"~ "bio",
                               itemB == "get sick"~ "bio")) 

df_combined <- df_freesort %>%
  # left_join(df_causal, by = c("subject_id", "itemA", "itemB")) # %>% 
  left_join(df_causal, by = c("subject_id", "itemA", "itemB", "domain_itemA", "domain_itemB")) %>%
  mutate(domain_itemA = factor(domain_itemA, levels = domain_order),
         domain_itemB = factor(domain_itemB, levels = domain_order)) %>% 
  relocate(freesort_distance, .before = "causal_distance") %>%
  relocate(trial_type_order, .after = "domain_itemB")

paged_table(df_combined)
write.csv(df_combined, data_processed_output_path, row.names = F)

2 Freesort RDMs

2.1 plot subject-specific raw responses

d %>% 
  filter(trial_type == "free_sort") %>%
  # filter(subject_id == "subj_03") %>% 
  mutate(y = 500 - y) %>% #invert y axis
  ggplot(aes(x = x, y = y)) +
  geom_point() +
  geom_text(aes(label = item), size = 3, vjust = -.2) +
  facet_wrap(~subject_id)

2.2 plot average raw responses

# pdf("my_plot_average.pdf", width = 7, height = 4)  # Width and height are in inches
d %>% 
  filter(trial_type == "free_sort") %>%
  mutate(y = 500 - y) %>% 
  group_by(item) %>% 
  mutate(mean_x = mean(x),
            mean_y = mean(y)) %>% 
  ungroup() %>% 
  distinct(item, .keep_all = TRUE) %>% 
  select(-subject_id) %>% 
  ggplot(aes(x = mean_x, y = mean_y)) +
  geom_point() +
  geom_text(aes(label = item), vjust = -1) +
  labs(title = "Average Raw Responses")

2.3 plot subject-specific Sorting RDMs

plot_rdm <- function(data, subject, desired_order) {
  df_matrix <- data %>%
    filter(subject_id == subject) %>%
    select(-subject_id) %>%
    select(itemA, itemB, freesort_distance) %>% 
    pivot_wider(names_from = itemB, values_from = freesort_distance) %>%
    column_to_rownames("itemA")
  
  df_matrix <- df_matrix[desired_order, desired_order]
  diag(df_matrix) <- 0 
  
  custom_palette <- colorRampPalette(c("red", "white", "blue"))(200)  # Generate 200 color levels

  corrplot(as.matrix(df_matrix), method = "color", is.corr = FALSE,
    col = custom_palette,  # Apply custom colors
    tl.col = "black", title = paste("Freesort RDM - ", subject),
    mar = c(0, 0, 2, 0), addgrid.col = "darkgray"
  )
}

subject_ids <- unique(df_freesort$subject_id)

#generate plots
# plots <- map(subject_ids, ~plot_rdm(df_freesort, .x, desired_order)) #run this to visualize plots in markdown  output

#save plots with subject name as file name
new_names <- paste0("subj_", str_pad(seq_along(subject_ids), width = 2, pad = "0")) 

#loop to save each plot as a separate PDF
for (i in seq_along(subject_ids)) {  
  pdf(file = paste0(base_path, "/study1-measure/figures/freesort/", new_names[i], ".pdf"), width = 8, height = 6)  #open a PDF device
  plot_rdm(df_freesort, subject_ids[i], desired_order)  
  dev.off()  #close the device to save the file
}

2.4 figure for paper: plot average freesort rdm

df_mean_freesort <- df_freesort %>% 
  unite("item_pairs", itemA, itemB, sep = "_", remove = FALSE) %>% 
  group_by(item_pairs) %>% 
  mutate(mean_freesort_distance = mean(freesort_distance)) %>% 
  ungroup()  
  
# matrixify and plot
df_mean_freesort_rdm <- df_mean_freesort %>% 
  distinct(item_pairs, .keep_all = TRUE) %>% 
  select(itemA, itemB, mean_freesort_distance) %>% 
  pivot_wider(names_from = itemB, values_from = mean_freesort_distance) %>% 
  column_to_rownames("itemA")  
  
df_mean_freesort_rdm <- df_mean_freesort_rdm[desired_order, desired_order]

diag(df_mean_freesort_rdm) <- NA #set diagonals 

#min max normalize
df_mean_freesort_matrix <- as.matrix(df_mean_freesort_rdm) 
min_val <- min(df_mean_freesort_matrix, na.rm = TRUE) 
max_val <- max(df_mean_freesort_matrix, na.rm = TRUE) 
df_mean_freesort_matrix <- (df_mean_freesort_matrix - min_val) / (max_val - min_val) 

#group_colors <- c("#2F2585", "#69ad44", "#f36b2d") 
group_colors <- c("#db9728", "#d0495a", "#00788a") 

#map each item to a color based on its group
item_colors <- setNames(c(rep(group_colors[1], 5), 
                          rep(group_colors[2], 5), 
                          rep(group_colors[3], 5)), 
                        desired_order)

#ensure colors align with the matrix order
label_colors <- item_colors[rownames(df_mean_freesort_matrix)]

custom_palette <- colorRampPalette(c("red", "white", "blue")) 

# pdf(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM.pdf"), width = 12, height = 12)

# pdf(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM.pdf"), width = 12, height = 12)
# 
# cairo_pdf(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM.pdf"),
#           width = 12, height = 12, bg = "transparent")

# png(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM.png"), width = 12, height = 12, bg = "transparent")  # For raster


png(
  file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM.png"),
  width = 12,
  height = 12,
  units = "in",
  res = 300,
  bg = "transparent"
)

#large image
corrplot(
  df_mean_freesort_matrix,
  method = "color",
  is.corr = FALSE,
  col = custom_palette(200),
  #tl.col = "black", #could be label_colors
  tl.col = label_colors, #could be label_colors
  tl.cex = 2,
  na.label = "square",       
  na.label.col = "darkgrey",   
  #title = "Average Freesort RDM",
  mar = c(0, 0, 2, 0),
  addgrid.col = "darkgray",
  cl.offset = 1.0,  # Push labels even further
  cl.cex = 2,     # Slightly smaller label text
  cl.ratio = 0.22,   # Widen the legend bar. Change this to push words far away
  cex.main = 3
)

#run this with the above to include legend for the domain colors
legend(
  "topleft",
  inset = c(0.1, 0.1), #move the legend outside the plot
  legend = c("Mind", "Action", "Body"),
  fill = group_colors,
  border = NA,
  bty = "n", #no box around the legend
  cex = 2, #text size
  xpd = TRUE, #allow drawing outside the plot area
  #ncol = 1,
  y.intersp = 1.4 #increase this for more space between rows
)

2.5 specific similarity relations reported in paper

compute_ci_freesort <- function(df, condition_label) {
  mean_distance <- mean(df$freesort_distance, na.rm = TRUE)
  sem <- sd(df$freesort_distance, na.rm = TRUE) / sqrt(nrow(df))
  ci_lower <- mean_distance - qt(0.975, df = nrow(df) - 1) * sem
  ci_upper <- mean_distance + qt(0.975, df = nrow(df) - 1) * sem

  cat(sprintf("%s: mean distance = %.3f [95%% CI: %.3f, %.3f]\n", 
              condition_label, mean_distance, ci_lower, ci_upper))
}

#compute and print within domain similarities
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "mind", domain_itemB == "mind"), "mind → mind")
## mind → mind: mean distance = 0.355 [95% CI: 0.338, 0.373]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "action", domain_itemB == "action"), "action → action")
## action → action: mean distance = 0.370 [95% CI: 0.351, 0.390]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "bio", domain_itemB == "bio"), "bio → bio")
## bio → bio: mean distance = 0.368 [95% CI: 0.350, 0.386]
#compute and print cross domain similarities
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "action", domain_itemB == "mind"), "mind → mind")
## mind → mind: mean distance = 0.565 [95% CI: 0.551, 0.580]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "bio", domain_itemB == "mind"), "action → action")
## action → action: mean distance = 0.574 [95% CI: 0.561, 0.587]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "bio", domain_itemB == "action"), "bio → bio")
## bio → bio: mean distance = 0.582 [95% CI: 0.568, 0.596]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "mind", domain_itemB == "action"), "bio → bio")
## bio → bio: mean distance = 0.565 [95% CI: 0.551, 0.580]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "mind", domain_itemB == "bio"), "bio → bio")
## bio → bio: mean distance = 0.574 [95% CI: 0.561, 0.587]
compute_ci_freesort(df_freesort %>% filter(domain_itemA == "action", domain_itemB == "bio"), "bio → bio")
## bio → bio: mean distance = 0.582 [95% CI: 0.568, 0.596]
#average cross-domain similarity
(.565 + .574 + .582 + .565 + .574 + .583) / 6
## [1] 0.5738333
#average within domain similarity
(.355 + .37 + .368) / 3
## [1] 0.3643333

2.5.1 plot only upper tri

This plot is not illustrated in the paper

pdf(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM-uppertri.pdf"), width = 12, height = 12)

#large image
corrplot(
  as.matrix(df_mean_freesort_matrix),
  method = "color",
  is.corr = FALSE,
  col = custom_palette(200),
  tl.col = label_colors, #could be label_colors
  tl.cex = 2,
  na.label = "square",       
  na.label.col = "darkgrey",   
  #title = "Average Freesort RDM",
  mar = c(0, 0, 2, 0),
  addgrid.col = "darkgray",
  cl.offset = 1.0,  #push labels even further
  cl.cex = 2,     #slightly smaller label text
  cl.ratio = 0.22,   #widen the legend bar. Change this to push words far away
  cex.main = 3,
  type = "upper"  #show only upper triangular matrix
)

2.5.2 lower tri

This plot is not illustrated in the paper. Upper and lower tri identical since the matrix is symmetrical.

pdf(file = paste0(rdms_outputs_path, "/freesort/averageFreesortRDM-lowertri.pdf"), width = 12, height = 12)

#large image
corrplot(
  as.matrix(df_mean_freesort_matrix),
  method = "color",
  is.corr = FALSE,
  col = custom_palette(200),
  #tl.col = "black", #could be label_colors
  # tl.col = label_colors, #could be label_colors
  # tl.cex = 2,
  tl.pos = "n",
  na.label = "square",       
  na.label.col = "darkgrey",   
  #title = "Freesort RDM lower tr",
  mar = c(0, 0, 2, 0),
  addgrid.col = "darkgray",
  cl.offset = 1.0,  #push labels even further
  cl.cex = 2,     #slightly smaller label text
  cl.ratio = 0.22,   #widen the legend bar. Change this to push words far away
  cex.main = 3,
  type = "lower"  #show only lower triangular matrix
)

3 Causal RDMs

3.1 plot subject-specific causal RDMs

plot_causal_rdm <- function(data, subject, desired_order) {
  
  df_matrix <- data %>%
    filter(subject_id == subject) %>%
    select(-subject_id) %>%
    select(itemA, itemB, causal_distance) %>% 
    pivot_wider(names_from = itemB, values_from = causal_distance) %>%
    column_to_rownames("itemA")
  
  df_matrix <- df_matrix[desired_order, desired_order]
  diag(df_matrix) <- 0
  
  custom_palette <- colorRampPalette(c("red", "white", "blue"))(200)

  
  corrplot(as.matrix(df_matrix), method = "color", is.corr = FALSE,
               col = custom_palette,  # Apply custom colors
           tl.col = "black", title = paste("Causal RDM - ", subject),
           mar = c(0, 0, 2, 0), addgrid.col = "darkgray")
}

subject_ids <- unique(df_causal$subject_id)

#generate plots
# plots <- map(subject_ids, ~plot_causal_rdm(df_causal, .x, desired_order)) #run this to generate plots

# save plots with subject name as file name
new_names <- paste0("subj_", str_pad(seq_along(subject_ids), width = 2, pad = "0"))

for (i in seq_along(subject_ids)) { 
    pdf(file = paste0(base_path, "/study1-measure/figures/causal/", new_names[i], ".pdf"), width = 8, height = 6)
  plot_causal_rdm(df_causal, subject_ids[i], desired_order)  
  dev.off()  
}

3.2 figure for paper: plot average causal RDM

df_mean_causal <- df_causal %>% 
  mutate(item_pairs = str_c(itemA, itemB, sep = "-"))  %>%
  group_by(item_pairs) %>%
  mutate(mean_causal_distance = mean(causal_distance, na.rm = TRUE)) %>%  #average across participants
  ungroup() 

df_mean_causal_rdm <- df_mean_causal %>%
  distinct(item_pairs, .keep_all = TRUE) %>% 
  select(itemA, itemB, mean_causal_distance) %>% 
  pivot_wider(names_from = itemB, values_from = mean_causal_distance) %>%
  column_to_rownames("itemA")  

df_mean_causal_rdm <- df_mean_causal_rdm[desired_order, desired_order]

diag(df_mean_causal_rdm) <- NA 

#min max normalize
df_mean_causal_matrix <- as.matrix(df_mean_causal_rdm) 
min_val <- min(df_mean_causal_matrix, na.rm = TRUE) 
max_val <- max(df_mean_causal_matrix, na.rm = TRUE) 
df_mean_causal_matrix <- (df_mean_causal_matrix - min_val) / (max_val - min_val) 

#colors
group_colors <- c("#db9728", "#d0495a", "#00788a") 


item_colors <- setNames(c(rep(group_colors[1], 5), #mind
                          rep(group_colors[2], 5), #action
                          rep(group_colors[3], 5)), #body
                        desired_order)

label_colors <- item_colors[rownames(df_mean_causal_rdm)]

custom_palette <- colorRampPalette(c("red", "white", "blue")) 

# pdf(file = paste0(rdms_outputs_path, "/causal/averageCausalRDM.pdf"), width = 12, height = 12)

png(
  file = paste0(rdms_outputs_path, "/causal/averageCausalRDM.png"),
  width = 12,
  height = 12,
  units = "in",
  res = 300,
  bg = "transparent"
)


#large image
corrplot(
  df_mean_causal_matrix,
  method = "color",
  is.corr = FALSE,
  col = custom_palette(200),
  #tl.col = "black", 
  tl.col = label_colors, 
  tl.cex = 2,
  mar = c(0, 0, 2, 0),
  addgrid.col = "darkgray",
  na.label = "square",       
  na.label.col = "darkgrey",
  cl.offset = 1.0, 
  cl.cex = 2, #label text
  cl.ratio = 0.22, #push words far away
  cex.main = 3
)

#run this together with above to include a domains color legend
legend(
  "topleft",
  inset = c(0.1, 0.1),
  legend = c("Mind", "Action", "Body"),
  fill = group_colors,
  border = NA,
  bty = "n",
  cex = 2, #text size
  xpd = TRUE, #draw outside the plot area
  #ncol = 1,
  y.intersp = 1.4 #more space between rows
)

3.3 specific causal relations reported in paper

compute_ci <- function(df, condition_label) {
  mean_distance <- mean(df$causal_distance, na.rm = TRUE)
  sem <- sd(df$causal_distance, na.rm = TRUE) / sqrt(nrow(df))
  ci_lower <- mean_distance - qt(0.975, df = nrow(df) - 1) * sem
  ci_upper <- mean_distance + qt(0.975, df = nrow(df) - 1) * sem

  cat(sprintf("%s: mean distance = %.3f [95%% CI: %.3f, %.3f]\n", 
              condition_label, mean_distance, ci_lower, ci_upper))
}


#defining and computing each condition
compute_ci(df_causal %>% filter(domain_itemA == "mind", domain_itemB == "action"), "mind → action")
## mind → action: mean distance = 0.323 [95% CI: 0.309, 0.337]
compute_ci(df_causal %>% filter(domain_itemA == "action", domain_itemB == "mind"), "action → mind")
## action → mind: mean distance = 0.415 [95% CI: 0.399, 0.431]
compute_ci(df_causal %>% filter(domain_itemA == "action", domain_itemB == "action"), "action → action")
## action → action: mean distance = 0.611 [95% CI: 0.592, 0.630]
compute_ci(df_causal %>% filter(domain_itemA == "bio", domain_itemB == "bio"), "body → body")
## body → body: mean distance = 0.381 [95% CI: 0.362, 0.400]
compute_ci(df_causal %>% filter(domain_itemA == "bio", domain_itemB == "action"), "body → action")
## body → action: mean distance = 0.442 [95% CI: 0.424, 0.460]
compute_ci(df_causal %>% filter(domain_itemA == "action", domain_itemB == "bio"), "action → body")
## action → body: mean distance = 0.427 [95% CI: 0.410, 0.445]
#Perception → Cognition
compute_ci(df_causal %>% 
             filter(domain_itemA == "mind", domain_itemB == "mind", 
                    itemA %in% c("hear something", "see something"),
                    itemB %in% c("choose what to do", "think about something", "remember something")),
           "perception → cognition")
## perception → cognition: mean distance = 0.149 [95% CI: 0.130, 0.167]
#Cognition → Perception
compute_ci(df_causal %>% 
             filter(domain_itemA == "mind", domain_itemB == "mind", 
                    itemB %in% c("hear something", "see something"),
                    itemA %in% c("choose what to do", "think about something", "remember something")),
           "cognition → perception")
## cognition → perception: mean distance = 0.477 [95% CI: 0.442, 0.512]
#Seeing → Sick
compute_ci(df_causal %>% filter(itemA == "see something", itemB == "get sick"), "seeing → sick")
## seeing → sick: mean distance = 0.390 [95% CI: 0.300, 0.481]
#Hearing → Sick
compute_ci(df_causal %>% filter(itemA == "hear something", itemB == "get sick"), "hearing → sick")
## hearing → sick: mean distance = 0.606 [95% CI: 0.514, 0.697]

4 Correlations: Sorting vs Causal Task

4.1 mean within-subject correlation

compute_subjectwise_rdm_correlation <- function(df_causal, df_freesort, desired_order, metric = "kendall") {
  subject_ids <- unique(df_causal$subject_id)
  subject_correlations <- numeric(length(subject_ids))
  
  for (i in seq_along(subject_ids)) {
    subject <- subject_ids[i]
    
    freesort_rdm <- df_freesort %>%
      filter(subject_id == subject) %>%
      select(itemA, itemB, freesort_distance) %>%
      pivot_wider(names_from = itemB, values_from = freesort_distance) %>%
      column_to_rownames("itemA") %>%
      as.matrix()
      diag(freesort_rdm) <- NA
    freesort_rdm <- freesort_rdm[desired_order, desired_order]
    
    causal_rdm <- df_causal %>%
      filter(subject_id == subject) %>%
      select(itemA, itemB, causal_distance) %>%
      pivot_wider(names_from = itemB, values_from = causal_distance) %>%
      column_to_rownames("itemA") %>%
      as.matrix()
        diag(causal_rdm) <- NA
    causal_rdm <- causal_rdm[desired_order, desired_order]
    
    freesort_flat <- as.numeric(freesort_rdm)
    causal_flat <- as.numeric(causal_rdm)
    
    subject_correlations[i] <- cor(freesort_flat, causal_flat, method = metric, use = "complete.obs")
  }
  
  return(subject_correlations)
}

subject_correlations <- compute_subjectwise_rdm_correlation(df_causal, df_freesort, desired_order)
hist(subject_correlations)

mean_correlation <- mean(subject_correlations, na.rm = TRUE)
mean_correlation
## [1] 0.03333567

4.2 t-test

t_test_result <- t.test(subject_correlations, mu = 0)
t_test_result
## 
##  One Sample t-test
## 
## data:  subject_correlations
## t = 3.7597, df = 49, p-value = 0.0004539
## alternative hypothesis: true mean is not equal to 0
## 95 percent confidence interval:
##  0.01551768 0.05115366
## sample estimates:
##  mean of x 
## 0.03333567

4.3 t-test assumptions

seem fine, so no need for perm test

shapiro_result <- shapiro.test(subject_correlations)
print(shapiro_result)
## 
##  Shapiro-Wilk normality test
## 
## data:  subject_correlations
## W = 0.9804, p-value = 0.5689
qqnorm(subject_correlations)
qqline(subject_correlations, col = "red")

hist(subject_correlations, main = "Histogram of Subject Correlations", 
     xlab = "Correlation", col = "lightblue", border = "black")

4.4 correlation between mean rdms

Note that this shows mean rdms are uncorrelated, while previously, we have seen a small positive correlation when we run subject-level correlations

#using Pearson's correlation (non-significant)
mean_rdm_correlations <- cor(as.numeric(df_mean_freesort_matrix), 
                             as.numeric(df_mean_causal_matrix), 
                             method = "pearson", use = "complete.obs")
cor.test(as.numeric(df_mean_freesort_matrix), 
         as.numeric(df_mean_causal_matrix), 
         method = "pearson", 
         use = "complete.obs")
## 
##  Pearson's product-moment correlation
## 
## data:  as.numeric(df_mean_freesort_matrix) and as.numeric(df_mean_causal_matrix)
## t = 0.78819, df = 208, p-value = 0.4315
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.08142201  0.18856724
## sample estimates:
##        cor 
## 0.05456999
#using Kendall's correlation (also non-significant)
n_permutations <- 10000
observed_tau <- cor(as.numeric(df_mean_freesort_matrix),
                    as.numeric(df_mean_causal_matrix),
                    method = "kendall",
                    use = "complete.obs")
observed_tau 
## [1] 0.01800242
null_dist <- replicate(n_permutations, {
  permuted <- sample(as.numeric(df_mean_causal_matrix))
  cor(as.numeric(df_mean_freesort_matrix), permuted, method = "kendall", use = "complete.obs")
})
p_value <- mean(null_dist >= observed_tau)  #one-sided
p_value
## [1] 0.3553

5 Correlations: Data vs Models

5.1 general functions for correlation plot

extract_upper_tri <- function(matrix) {
  matrix[upper.tri(matrix)]
}

compute_rdm_correlation <- function(data, emp_col, theor_col, desired_order, method = "kendall") {
  rdm1 <- data %>%
    select(itemA, itemB, {{ emp_col }}) %>%
    pivot_wider(names_from = itemB, values_from = {{ emp_col }}) %>%
    column_to_rownames("itemA") %>%
    as.matrix()
  
  rdm2 <- data %>%
    select(itemA, itemB, {{ theor_col }}) %>%
    pivot_wider(names_from = itemB, values_from = {{ theor_col }}) %>%
    column_to_rownames("itemA") %>%
    as.matrix()
  
  rdm1 <- rdm1[desired_order, desired_order]
  rdm2 <- rdm2[desired_order, desired_order]
  
  diag(rdm1) <- NA
  diag(rdm2) <- NA
  
  upper_rdm1 <- extract_upper_tri(rdm1)
  upper_rdm2 <- extract_upper_tri(rdm2)
  
  cor(upper_rdm1, upper_rdm2, method = method, use = "complete.obs")
}

5.2 noise ceiling functions

calculate_freesort_noise_ceiling <- function(data, subjects, desired_order, metric = "kendall") {
  # Initialize storage for correlations
  upper_bound <- numeric(length(subjects))
  lower_bound <- numeric(length(subjects))
  
  for (i in seq_along(subjects)) {
    subject <- subjects[i]
    
    #Subject's RDM
    subject_rdm <- data %>%
        filter(subject_id == subject) %>%
        select(subject_id, itemA, itemB, freesort_distance) %>% 
        pivot_wider(names_from = itemB, values_from = freesort_distance) %>%
        column_to_rownames("itemA") %>%
        as.matrix()
    subject_rdm <- subject_rdm[desired_order, desired_order]
    
    #Grand average RDM
    grand_avg_rdm <- data %>%
        select(subject_id, itemA, itemB, freesort_distance) %>% 
        unite("item_pairs", itemA, itemB, sep = "_", remove = FALSE) %>% 
        group_by(item_pairs) %>%
        mutate(freesort_distance = mean(freesort_distance)) %>%
        ungroup() %>% 
        distinct(item_pairs, .keep_all = TRUE) %>% 
        select(-c(item_pairs, subject_id)) %>% 
        pivot_wider(names_from = itemB, values_from = freesort_distance) %>%
        column_to_rownames("itemA") %>%
        as.matrix()
    grand_avg_rdm <- grand_avg_rdm[desired_order, desired_order]
    
    #Other subjects' average RDM
    other_avg_rdm <- data %>%
      select(subject_id, itemA, itemB, freesort_distance) %>% 
      filter(subject_id != subject) %>%
      unite("item_pairs", itemA, itemB, sep = "_", remove = FALSE) %>% 
      group_by(item_pairs) %>%
      mutate(freesort_distance = mean(freesort_distance)) %>%
      ungroup() %>% 
      distinct(item_pairs, .keep_all = TRUE) %>% 
      select(-c(item_pairs, subject_id)) %>% 
      pivot_wider(names_from = itemB, values_from = freesort_distance) %>%
      column_to_rownames("itemA") %>%
      as.matrix()
    other_avg_rdm <- other_avg_rdm[desired_order, desired_order]
    
    subject_upper <- as.numeric(subject_rdm[upper.tri(subject_rdm)])
    other_upper <- as.numeric(other_avg_rdm[upper.tri(other_avg_rdm)])
    grand_upper <- as.numeric(grand_avg_rdm[upper.tri(grand_avg_rdm)])
    
    lower_bound[i] <- cor(subject_upper, other_upper, method = metric) 
    upper_bound[i] <- cor(subject_upper, grand_upper, method = metric) 
  }
  
  # Return average bounds
  list(
    upper_bound = mean(upper_bound, na.rm = TRUE),
    lower_bound = mean(lower_bound, na.rm = TRUE)
  )
}

calculate_causal_noise_ceiling <- function(data, subjects, desired_order, metric = "kendall") {
  #Initialize storage for correlations
  upper_tri_upper_bound <- numeric(length(subjects))
  upper_tri_lower_bound <- numeric(length(subjects))
  
  lower_tri_upper_bound <- numeric(length(subjects))
  lower_tri_lower_bound <- numeric(length(subjects))
  
  #Loop through each subject
  for (i in seq_along(subjects)) {
    subject <- subjects[i]
    
    #subject's RDM
    subject_rdm <- data %>%
        filter(subject_id == subject) %>%
        select(subject_id, itemA, itemB, causal_distance) %>% 
        pivot_wider(names_from = itemB, values_from = causal_distance) %>%
        column_to_rownames("itemA") %>%
        as.matrix()
    subject_rdm <- subject_rdm[desired_order, desired_order]
    
    #grand average RDM
    grand_avg_rdm <- data %>%
        select(subject_id, itemA, itemB, causal_distance) %>% 
        unite("item_pairs", itemA, itemB, sep = "_", remove = FALSE) %>% 
        group_by(item_pairs) %>%
        mutate(causal_distance = mean(causal_distance)) %>%
        ungroup() %>% 
        distinct(item_pairs, .keep_all = TRUE) %>% 
        select(-c(item_pairs, subject_id)) %>% 
        pivot_wider(names_from = itemB, values_from = causal_distance) %>%
        column_to_rownames("itemA") %>%
        as.matrix()
    grand_avg_rdm <- grand_avg_rdm[desired_order, desired_order]
    
    #other subjects' average RDM
    other_avg_rdm <- data %>%
      select(subject_id, itemA, itemB, causal_distance) %>% 
      filter(subject_id != subject) %>%
      unite("item_pairs", itemA, itemB, sep = "_", remove = FALSE) %>% 
      group_by(item_pairs) %>%
      mutate(causal_distance = mean(causal_distance)) %>%
      ungroup() %>% 
      distinct(item_pairs, .keep_all = TRUE) %>% 
      select(-c(item_pairs, subject_id)) %>% 
      pivot_wider(names_from = itemB, values_from = causal_distance) %>%
      column_to_rownames("itemA") %>%
      as.matrix()
    other_avg_rdm <- other_avg_rdm[desired_order, desired_order]
    
    #upper tri noise ceiling
    subject_upper <- as.numeric(subject_rdm[upper.tri(subject_rdm)])
    other_upper <- as.numeric(other_avg_rdm[upper.tri(other_avg_rdm)])
    grand_upper <- as.numeric(grand_avg_rdm[upper.tri(grand_avg_rdm)])

    upper_tri_lower_bound[i] <- cor(subject_upper, other_upper, method = metric)
    upper_tri_upper_bound[i] <- cor(subject_upper, grand_upper, method = metric)

    #lower tri noise ceiling
    subject_lower <- as.numeric(subject_rdm[lower.tri(subject_rdm)])
    other_lower <- as.numeric(other_avg_rdm[lower.tri(other_avg_rdm)])
    grand_lower <- as.numeric(grand_avg_rdm[lower.tri(grand_avg_rdm)])
    
    lower_tri_lower_bound[i] <- cor(subject_lower, other_lower, method = metric)
    lower_tri_upper_bound[i] <- cor(subject_lower, grand_lower, method = metric)

  }
  
  #return average bounds
  list(
    upper_tri_upper_bound = mean(upper_tri_upper_bound, na.rm = TRUE),
    upper_tri_lower_bound = mean(upper_tri_lower_bound, na.rm = TRUE),
    
    lower_tri_upper_bound = mean(lower_tri_upper_bound, na.rm = TRUE),
    lower_tri_lower_bound = mean(lower_tri_lower_bound, na.rm = TRUE)
  )
}

5.3 calculate Sorting Task noise ceiling

#create correlations data frame
d_model_predictions <- read_csv(model_predictions_input_path)
## Rows: 210 Columns: 6
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (2): itemA, itemB
## dbl (4): theor1_Mind-Body-Action, theor2_PerCog-OdActSpAct-BodyStimBodynoSti...
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
d_freesort_correlations <- df_combined %>%
  left_join(d_model_predictions, by = c("itemA", "itemB")) %>%
  group_by(subject_id) %>%
  mutate(
    `1_Mind-Body-Action` = compute_rdm_correlation(pick(everything()), freesort_distance, `theor1_Mind-Body-Action`, desired_order),
    `2_PerCog-ObjDir-Stimulus` = compute_rdm_correlation(pick(everything()), freesort_distance, `theor2_PerCog-OdActSpAct-BodyStimBodynoStim`, desired_order),
    `3_Physical-Ethereal` = compute_rdm_correlation(pick(everything()), freesort_distance, `theor3_Physical-Ethereal`, desired_order),
    `4_cosine_similarity` = compute_rdm_correlation(pick(everything()), freesort_distance, `theor4_cosine_similarity`, desired_order),
    `5_causal_judgments` = compute_rdm_correlation(pick(everything()), freesort_distance, `causal_distance`, desired_order)
  ) %>%
  ungroup() %>%
  relocate(causal_distance, .after = "freesort_distance")

subjects_freesort <- unique(d_freesort_correlations$subject_id)
noise_ceiling_sorting <- calculate_freesort_noise_ceiling(d_freesort_correlations, subjects_freesort, desired_order)
noise_ceiling_sorting
## $upper_bound
## [1] 0.2587916
## 
## $lower_bound
## [1] 0.229508

5.4 calculate Causal Task noise ceiling

d_causal_correlations <-  df_combined %>%
  left_join(d_model_predictions, by = c("itemA", "itemB")) %>% 
  group_by(subject_id) %>%
  mutate(
    `1_Mind-Body-Action` = compute_rdm_correlation(pick(everything()), causal_distance, `theor1_Mind-Body-Action`, desired_order),
    `2_PerCog-ObjDir-Stimulus` = compute_rdm_correlation(pick(everything()), causal_distance, `theor2_PerCog-OdActSpAct-BodyStimBodynoStim`, desired_order),
    `3_Physical-Ethereal` = compute_rdm_correlation(pick(everything()), causal_distance, `theor3_Physical-Ethereal`, desired_order),
    `4_cosine_similarity` = compute_rdm_correlation(pick(everything()), causal_distance, `theor4_cosine_similarity`, desired_order),
    `5_causal_judgments` = compute_rdm_correlation(pick(everything()), causal_distance, `causal_distance`, desired_order)
  ) %>%
  ungroup() 

# Compute noise ceiling for all subjects
subjects_causal <- unique(d_causal_correlations$subject_id)
noise_ceiling_causal <- calculate_causal_noise_ceiling(d_causal_correlations,subjects_causal, desired_order)
noise_ceiling_causal
## $upper_tri_upper_bound
## [1] 0.4413385
## 
## $upper_tri_lower_bound
## [1] 0.4263549
## 
## $lower_tri_upper_bound
## [1] 0.4130283
## 
## $lower_tri_lower_bound
## [1] 0.3971654

5.5 figure: plot RSA corr between observed and model rdms

aggregated_freesort_correlations <- d_freesort_correlations %>% 
  rename("3 Category" = `1_Mind-Body-Action`,
         "6 Category" = `2_PerCog-ObjDir-Stimulus`,
         "2 Category" = `3_Physical-Ethereal`,
         "Cosine Similarity" = `4_cosine_similarity`) %>% 
  pivot_longer(
    cols = c("2 Category", "3 Category", "6 Category", "Cosine Similarity"),
    names_to = "correlation_type",
    values_to = "correlation_value"
  ) %>%
  group_by(subject_id, correlation_type) %>%
  summarize(mean_correlation = mean(correlation_value, na.rm = TRUE), .groups = "drop") %>% 
  mutate(method = "Sorting Task")

aggregated_causal_correlations <- d_causal_correlations %>% 
  rename("3 Category" = `1_Mind-Body-Action`,
         "6 Category" = `2_PerCog-ObjDir-Stimulus`,
         "2 Category" = `3_Physical-Ethereal`,
         "Cosine Similarity" = `4_cosine_similarity`) %>% 
  pivot_longer(
    cols = c("2 Category", "3 Category", "6 Category", "Cosine Similarity"),
    names_to = "correlation_type",
    values_to = "correlation_value"
  ) %>%
  group_by(subject_id, correlation_type) %>%
  summarize(mean_correlation = mean(correlation_value, na.rm = TRUE), .groups = "drop") %>% 
  mutate(method = "Causal Task")

#combine both corrs into one
d_combined_correlations <- bind_rows(aggregated_freesort_correlations, aggregated_causal_correlations) %>% 
  mutate(method = factor(method, levels = c("Sorting Task", "Causal Task"))) 

bar_plot <- ggplot(d_combined_correlations, 
                   aes(x = reorder(correlation_type, mean_correlation, FUN = function(x) -mean(x)), 
                        y = mean_correlation, fill = method)) +
  stat_summary(
    fun = "mean",
    geom = "bar",
    position = "dodge",  
    color = "black"
  ) +
  stat_summary(
    fun.data = mean_se,
    geom = "errorbar",
    width = 0.3,
    position = position_dodge(width = 0.9), 
    color = "black"
  ) +
  geom_point(aes(color = method), alpha = 0.3, size = 3, position = position_jitterdodge(jitter.width = 0.2)) +
  scale_fill_manual(values = c("Sorting Task" = "grey80", "Causal Task" = "#000000")) + #bar fill
  scale_color_manual(values = c("Sorting Task" = "darkgray", "Causal Task" = "#4A4A4A")) + #point colors
  
  geom_hline(yintercept = noise_ceiling_sorting$upper_bound, color = "darkgray", linetype = "dashed", linewidth = 0.5) +
  geom_hline(yintercept = noise_ceiling_sorting$lower_bound, color = "darkgray", linetype = "dashed", linewidth = 0.5) +
  geom_hline(yintercept = noise_ceiling_causal$upper_tri_upper_bound, color = "black", linetype = "dashed", linewidth = 0.5) +
  geom_hline(yintercept = noise_ceiling_causal$upper_tri_lower_bound, color = "black", linetype = "dashed", linewidth = 0.5) +
  theme_cowplot() +
  theme(
    panel.background = element_rect(fill = NA, color = NA),
    plot.background = element_rect(fill = NA, color = NA),
    axis.text.x = element_blank(),
    axis.text.y = element_text(size = 16),
    axis.ticks.length = unit(.25, "cm"),
    axis.title.x = element_blank(),
    axis.title.y = element_text(size = 18),
    legend.position = "top",
    legend.title = element_blank(),
    plot.title = element_text(size = 20, face = "bold", hjust = 0.5),
    plot.margin = margin(t = 10, r = 10, b = 20, l = 20)
  ) +
  labs(
    y = "Kendall's Tau Correlations",
    title = "Comparison of Theoretical and Empirical RDM Correlations"
  )

# Load and create RDM images
rdm_2_category <- ggdraw() + draw_image(paste0(rdms_input_path, "/theory_4.png"), scale = 1.2)
rdm_3_category <- ggdraw() + draw_image(paste0(rdms_input_path, "/theory_2.png"), scale = 1.2)
rdm_6_category <- ggdraw() + draw_image(paste0(rdms_input_path, "/theory_3.png"), scale = 1.2)
rdm_cosine <- ggdraw() + draw_image(paste0(rdms_input_path, "/theory_1.png"), scale = 1.2)

rdm_row <- rdm_3_category + rdm_6_category + rdm_2_category + rdm_cosine +
  plot_layout(ncol = 4)
# 
# plot_correlations <- bar_plot / rdm_row + 
#   plot_layout(heights = c(5, 2)) #adjusts height ratio


plot_correlations <- bar_plot / rdm_row +
  plot_layout(heights = c(5, 2)) &
  theme(plot.background = element_rect(fill = NA, color = NA))

plot_correlations

# ggsave(
#   filename = "figures/rdm_correlation_plot.pdf",
#   plot = plot_correlations,
#   device = "pdf",
#   #path = "path/to/directory",
#   width = 12,
#   height = 5,
#   units = "in"
# )

ggsave(
  filename = "figures/correlations/rdm_correlation_plot.png",  # ← use PNG instead of PDF
  plot = plot_correlations,
  device = "png",
  width = 12,
  height = 6,
  units = "in",
  dpi = 300,
  bg = "transparent"  # ← this is key
)

5.6 t-tests of correlations

This compared correlations to 0.

t_test_results <- d_combined_correlations %>%
  group_by(method, correlation_type) %>%
  dplyr::summarise(
    mean_mean_correlation = mean(mean_correlation, na.rm = TRUE),
    t_value = t.test(mean_correlation, mu = 0)$statistic,  
    p_value = t.test(mean_correlation, mu = 0)$p.value,
    .groups = "drop"
  )

print(t_test_results)
## # A tibble: 8 × 5
##   method       correlation_type  mean_mean_correlation t_value  p_value
##   <fct>        <chr>                             <dbl>   <dbl>    <dbl>
## 1 Sorting Task 2 Category                       0.144     6.00 2.35e- 7
## 2 Sorting Task 3 Category                       0.254     9.77 4.32e-13
## 3 Sorting Task 6 Category                       0.193    11.1  5.27e-15
## 4 Sorting Task Cosine Similarity                0.129     7.95 2.31e-10
## 5 Causal Task  2 Category                      -0.0823   -4.96 8.90e- 6
## 6 Causal Task  3 Category                      -0.0357   -2.93 5.10e- 3
## 7 Causal Task  6 Category                      -0.0569   -6.20 1.16e- 7
## 8 Causal Task  Cosine Similarity               -0.0797   -7.49 1.14e- 9

6 Linear models comparing the 4 hypotheses

6.1 freesort data

Reported in main paper

freesort_correlations <- d_freesort_correlations %>% 
  rename("mind_action_physiology_model" = `1_Mind-Body-Action`,
         "finegrained_mind_action_physiology_model" = `2_PerCog-ObjDir-Stimulus`,
         "physical_psychological_model" = `3_Physical-Ethereal`,
         "use_phrase_embeddings_model" = `4_cosine_similarity`) %>% 
  pivot_longer(
    cols = c("mind_action_physiology_model":"use_phrase_embeddings_model"),
    names_to = "RDM_type",
    values_to = "tau"
  ) %>%
  select(subject_id, itemA, itemB, domain_itemA, domain_itemB, trial_type_order, freesort_distance, causal_distance, RDM_type, tau) %>% 
  mutate(
    RDM_type = as.factor(RDM_type), 
    RDM_type = relevel(RDM_type, ref = "mind_action_physiology_model") 
  )

#ensure subject_id is a factor, for running simr later
freesort_correlations <- freesort_correlations %>%
  mutate(subject_id = as.factor(subject_id))

#summarize to 1 row per subject per model
subject_level_data <- freesort_correlations %>%
  distinct(subject_id, RDM_type, .keep_all = TRUE) %>%
  select(subject_id, RDM_type, tau)

#Note: mind action physiology is reference model. Estiamtes are mean kendall's taus
model_freesort <- lmerTest::lmer(tau ~ RDM_type + (1 | subject_id), 
              data = subject_level_data)

summary(model_freesort)
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: tau ~ RDM_type + (1 | subject_id)
##    Data: subject_level_data
## 
## REML criterion at convergence: -250.6
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.16101 -0.57970 -0.09903  0.53510  3.03194 
## 
## Random effects:
##  Groups     Name        Variance Std.Dev.
##  subject_id (Intercept) 0.013304 0.11534 
##  Residual               0.009363 0.09676 
## Number of obs: 200, groups:  subject_id, 50
## 
## Fixed effects:
##                                                   Estimate Std. Error        df
## (Intercept)                                        0.25383    0.02129  96.38760
## RDM_typefinegrained_mind_action_physiology_model  -0.06049    0.01935 147.00000
## RDM_typephysical_psychological_model              -0.11026    0.01935 147.00000
## RDM_typeuse_phrase_embeddings_model               -0.12482    0.01935 147.00000
##                                                  t value Pr(>|t|)    
## (Intercept)                                       11.922  < 2e-16 ***
## RDM_typefinegrained_mind_action_physiology_model  -3.126  0.00214 ** 
## RDM_typephysical_psychological_model              -5.698 6.42e-08 ***
## RDM_typeuse_phrase_embeddings_model               -6.450 1.52e-09 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##             (Intr) RDM_____ RDM_t__
## RDM_typ____ -0.454                 
## RDM_typph__ -0.454  0.500          
## RDM_typs___ -0.454  0.500    0.500

6.2 getting CIs

#for mind action and physiology model
d_target <- d_combined_correlations %>%
  filter(correlation_type == "3 Category", method == "Sorting Task")

t_test <- t.test(d_target$mean_correlation, mu = 0)

t_test$estimate #mean Kendall's tau
## mean of x 
## 0.2538288
t_test$conf.int 
## [1] 0.2016248 0.3060328
## attr(,"conf.level")
## [1] 0.95

7 Causal vs Freesort Models

7.1 wrangle data

Testing if the best model for the Freesort data (i.e. if Mind, action and physiology is good or bad at describing the Causal data)

taus_causal <- d_causal_correlations %>% 
  rename("mind_action_physiology_model" = `1_Mind-Body-Action`,
         "finegrained_mind_action_physiology_model" = `2_PerCog-ObjDir-Stimulus`,
         "physical_psychological_model" = `3_Physical-Ethereal`,
         "use_phrase_embeddings_model" = `4_cosine_similarity`) %>% 
  pivot_longer(
    cols = c("mind_action_physiology_model":"use_phrase_embeddings_model"),
    names_to = "RDM_type",
    values_to = "tau"
  ) %>%
  filter(RDM_type == "mind_action_physiology_model") %>% 
  rename("causal" = tau) %>% 
  distinct(subject_id, .keep_all = TRUE) %>% 
  select(subject_id, causal)

taus_freesort <- d_freesort_correlations %>% 
  rename("mind_action_physiology_model" = `1_Mind-Body-Action`,
         "finegrained_mind_action_physiology_model" = `2_PerCog-ObjDir-Stimulus`,
         "physical_psychological_model" = `3_Physical-Ethereal`,
         "use_phrase_embeddings_model" = `4_cosine_similarity`) %>% 
  pivot_longer(
    cols = c("mind_action_physiology_model":"use_phrase_embeddings_model"),
    names_to = "RDM_type",
    values_to = "tau"
  ) %>%
  filter(RDM_type == "mind_action_physiology_model") %>% 
  rename("freesort" = tau) %>% 
  distinct(subject_id, .keep_all = TRUE) %>% 
  select(subject_id, freesort)

combined_taus <- taus_causal %>% 
  left_join(taus_freesort,
            by = "subject_id") %>% 
  pivot_longer(cols = c(causal, freesort),
               names_to = "task",
               values_to = "tau") %>% 
  mutate(
    task = as.factor(task),  
    task = relevel(task, ref = "freesort") 
  ) %>% 
  mutate(subject_id = as.factor(subject_id))

model_prediction2 <- lmerTest::lmer(tau ~ task + (1|subject_id), data = combined_taus)
summary(model_prediction2)
## Linear mixed model fit by REML. t-tests use Satterthwaite's method [
## lmerModLmerTest]
## Formula: tau ~ task + (1 | subject_id)
##    Data: combined_taus
## 
## REML criterion at convergence: -95.4
## 
## Scaled residuals: 
##      Min       1Q   Median       3Q      Max 
## -2.16565 -0.55581  0.01977  0.50180  2.61279 
## 
## Random effects:
##  Groups     Name        Variance Std.Dev.
##  subject_id (Intercept) 0.002424 0.04924 
##  Residual               0.018143 0.13470 
## Number of obs: 100, groups:  subject_id, 50
## 
## Fixed effects:
##             Estimate Std. Error       df t value Pr(>|t|)    
## (Intercept)  0.25383    0.02028 96.65704   12.52  < 2e-16 ***
## taskcausal  -0.28949    0.02694 49.00000  -10.75 1.75e-14 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Correlation of Fixed Effects:
##            (Intr)
## taskcausal -0.664

7.2 plot: Map model explaining causal vs freesort

Compare Mind action physiology fit with freesort vs causal data

ggplot(combined_taus, aes(x = task, y = tau)) +
  geom_boxplot(outlier.shape = NA, fill = "lightgray") +
  geom_jitter(width = 0.2, alpha = 0.5) + 
  labs(title = "Comparison of Tau Values by Task",
       x = "Task",
       y = "Tau") +
  scale_fill_manual(values = c("Sorting Task" = "darkgray", "Causal Task" = "black")) + #bar fill
  scale_color_manual(values = c("Sorting Task" = "lightgray", "Causal Task" = "#4A4A4A")) +
  theme_bw()
## Warning: No shared levels found between `names(values)` of the manual scale and the
## data's fill values.
## Warning: No shared levels found between `names(values)` of the manual scale and the
## data's colour values.

ggplot(combined_taus, aes(x = task, y = tau, fill = task)) +
  geom_boxplot(outlier.shape = NA, alpha = 0.6) +
  geom_jitter(width = 0.2, alpha = 0.5, aes(color = task)) + 
  scale_fill_manual(values = c("freesort" = "darkgray", "causal" = "black")) +  
  scale_color_manual(values = c("freesort" = "darkgray", "causal" = "black")) + 
  labs(title = "Comparison of Tau Values by Task",
       x = "Task",
       y = "Tau") +
  theme_bw() +
  theme(legend.position = "none",
        axis.text.x = element_text(size = 14),
        axis.text.y = element_text(size = 14), 
        axis.title.x = element_text(size = 16), 
        axis.title.y = element_text(size = 16)) 

8 Power analysis using simr

For planning follow-up studies (i.e. replicating study 1 and 2).

8.1 simulate effect of all 4 models

# #suppose we want a study with 100 participants..
model_sim <- extend(model_freesort, along = "subject_id", n = 100)
length(unique(simr::getData(model_sim)$subject_id)) #confirm that N is now 100
## [1] 100
# Leave this code commented out
# #and we want to replicate the lower limit effect size, assuming pilot is inflated..
# ## fixef(model_sim)["(Intercept)"] <- 0.202

# #what would the power be?
#sim_inercept <- powerSim(model_sim, fixed("(Intercept)", "t"), nsim = 100)
#sim_inercept


#========OR 

#Uncomment this and run (if commented out)
 
# #what would the power curve be for different sample sizes for the intercept?
# #intercept
pbapply::pboptions(type = "none")   # disables the progress bars

curve_intercept <- powerCurve(model_sim, fixed("(Intercept)", "t"),
                    along = "subject_id",
                    breaks = c(2, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
                    nsim = 100)
## Simulating: |                                                                  |Simulating: |=                                                                 |Simulating: |==                                                                |Simulating: |===                                                               |Simulating: |====                                                              |Simulating: |=====                                                             |Simulating: |======                                                            |Simulating: |=======                                                           |Simulating: |========                                                          |Simulating: |=========                                                         |Simulating: |==========                                                        |Simulating: |===========                                                       |Simulating: |============                                                      |Simulating: |=============                                                     |Simulating: |==============                                                    |Simulating: |===============                                                   |Simulating: |================                                                  |Simulating: |=================                                                 |Simulating: |==================                                                |Simulating: |===================                                               |Simulating: |====================                                              |Simulating: |=====================                                             |Simulating: |======================                                            |Simulating: |=======================                                           |Simulating: |========================                                          |Simulating: |=========================                                         |Simulating: |==========================                                        |Simulating: |===========================                                       |Simulating: |============================                                      |Simulating: |=============================                                     |Simulating: |==============================                                    |Simulating: |===============================                                   |Simulating: |================================                                  |Simulating: |=================================                                 |Simulating: |==================================                                |Simulating: |===================================                               |Simulating: |====================================                              |Simulating: |=====================================                             |Simulating: |======================================                            |Simulating: |=======================================                           |Simulating: |========================================                          |Simulating: |=========================================                         |Simulating: |==========================================                        |Simulating: |===========================================                       |Simulating: |============================================                      |Simulating: |=============================================                     |Simulating: |==============================================                    |Simulating: |===============================================                   |Simulating: |================================================                  |Simulating: |=================================================                 |Simulating: |==================================================                |Simulating: |===================================================               |Simulating: |====================================================              |Simulating: |=====================================================             |Simulating: |======================================================            |Simulating: |=======================================================           |Simulating: |========================================================          |Simulating: |=========================================================         |Simulating: |==========================================================        |Simulating: |===========================================================       |Simulating: |============================================================      |Simulating: |=============================================================     |Simulating: |==============================================================    |Simulating: |===============================================================   |Simulating: |================================================================  |Simulating: |================================================================= |Simulating: |==================================================================|( 1/12) ( 1/12) Simulating: |                                                          |( 1/12) Simulating: |=                                                         |( 1/12) Simulating: |==                                                        |( 1/12) Simulating: |===                                                       |( 1/12) Simulating: |====                                                      |( 1/12) Simulating: |=====                                                     |( 1/12) Simulating: |======                                                    |( 1/12) Simulating: |=======                                                   |( 1/12) Simulating: |========                                                  |( 1/12) Simulating: |=========                                                 |( 1/12) Simulating: |==========                                                |( 1/12) Simulating: |===========                                               |( 1/12) Simulating: |============                                              |( 1/12) Simulating: |=============                                             |( 1/12) Simulating: |==============                                            |( 1/12) Simulating: |===============                                           |( 1/12) Simulating: |================                                          |( 1/12) Simulating: |=================                                         |( 1/12) Simulating: |==================                                        |( 1/12) Simulating: |===================                                       |( 1/12) Simulating: |====================                                      |( 1/12) Simulating: |=====================                                     |( 1/12) Simulating: |======================                                    |( 1/12) Simulating: |=======================                                   |( 1/12) Simulating: |========================                                  |( 1/12) Simulating: |=========================                                 |( 1/12) Simulating: |==========================                                |( 1/12) Simulating: |===========================                               |( 1/12) Simulating: |============================                              |( 1/12) Simulating: |=============================                             |( 1/12) Simulating: |==============================                            |( 1/12) Simulating: |===============================                           |( 1/12) Simulating: |================================                          |( 1/12) Simulating: |=================================                         |( 1/12) Simulating: |==================================                        |( 1/12) Simulating: |===================================                       |( 1/12) Simulating: |====================================                      |( 1/12) Simulating: |=====================================                     |( 1/12) Simulating: |======================================                    |( 1/12) Simulating: |=======================================                   |( 1/12) Simulating: |========================================                  |( 1/12) Simulating: |=========================================                 |( 1/12) Simulating: |==========================================                |( 1/12) Simulating: |===========================================               |( 1/12) Simulating: |============================================              |( 1/12) Simulating: |=============================================             |( 1/12) Simulating: |==============================================            |( 1/12) Simulating: |===============================================           |( 1/12) Simulating: |================================================          |( 1/12) Simulating: |=================================================         |( 1/12) Simulating: |==================================================        |( 1/12) Simulating: |===================================================       |( 1/12) Simulating: |====================================================      |( 1/12) Simulating: |=====================================================     |( 1/12) Simulating: |======================================================    |( 1/12) Simulating: |=======================================================   |( 1/12) Simulating: |========================================================  |( 1/12) Simulating: |========================================================= |( 1/12) Simulating: |==========================================================|( 1/12) ( 2/12) ( 2/12) Simulating: |                                                          |( 2/12) Simulating: |=                                                         |( 2/12) Simulating: |==                                                        |( 2/12) Simulating: |===                                                       |( 2/12) Simulating: |====                                                      |( 2/12) Simulating: |=====                                                     |( 2/12) Simulating: |======                                                    |( 2/12) Simulating: |=======                                                   |( 2/12) Simulating: |========                                                  |( 2/12) Simulating: |=========                                                 |( 2/12) Simulating: |==========                                                |( 2/12) Simulating: |===========                                               |( 2/12) Simulating: |============                                              |( 2/12) Simulating: |=============                                             |( 2/12) Simulating: |==============                                            |( 2/12) Simulating: |===============                                           |( 2/12) Simulating: |================                                          |( 2/12) Simulating: |=================                                         |( 2/12) Simulating: |==================                                        |( 2/12) Simulating: |===================                                       |( 2/12) Simulating: |====================                                      |( 2/12) Simulating: |=====================                                     |( 2/12) Simulating: |======================                                    |( 2/12) Simulating: |=======================                                   |( 2/12) Simulating: |========================                                  |( 2/12) Simulating: |=========================                                 |( 2/12) Simulating: |==========================                                |( 2/12) Simulating: |===========================                               |( 2/12) Simulating: |============================                              |( 2/12) Simulating: |=============================                             |( 2/12) Simulating: |==============================                            |( 2/12) Simulating: |===============================                           |( 2/12) Simulating: |================================                          |( 2/12) Simulating: |=================================                         |( 2/12) Simulating: |==================================                        |( 2/12) Simulating: |===================================                       |( 2/12) Simulating: |====================================                      |( 2/12) Simulating: |=====================================                     |( 2/12) Simulating: |======================================                    |( 2/12) Simulating: |=======================================                   |( 2/12) Simulating: |========================================                  |( 2/12) Simulating: |=========================================                 |( 2/12) Simulating: |==========================================                |( 2/12) Simulating: |===========================================               |( 2/12) Simulating: |============================================              |( 2/12) Simulating: |=============================================             |( 2/12) Simulating: |==============================================            |( 2/12) Simulating: |===============================================           |( 2/12) Simulating: |================================================          |( 2/12) Simulating: |=================================================         |( 2/12) Simulating: |==================================================        |( 2/12) Simulating: |===================================================       |( 2/12) Simulating: |====================================================      |( 2/12) Simulating: |=====================================================     |( 2/12) Simulating: |======================================================    |( 2/12) Simulating: |=======================================================   |( 2/12) Simulating: |========================================================  |( 2/12) Simulating: |========================================================= |( 2/12) Simulating: |==========================================================|( 2/12) ( 3/12) ( 3/12) Simulating: |                                                          |( 3/12) Simulating: |=                                                         |( 3/12) Simulating: |==                                                        |( 3/12) Simulating: |===                                                       |( 3/12) Simulating: |====                                                      |( 3/12) Simulating: |=====                                                     |( 3/12) Simulating: |======                                                    |( 3/12) Simulating: |=======                                                   |( 3/12) Simulating: |========                                                  |( 3/12) Simulating: |=========                                                 |( 3/12) Simulating: |==========                                                |( 3/12) Simulating: |===========                                               |( 3/12) Simulating: |============                                              |( 3/12) Simulating: |=============                                             |( 3/12) Simulating: |==============                                            |( 3/12) Simulating: |===============                                           |( 3/12) Simulating: |================                                          |( 3/12) Simulating: |=================                                         |( 3/12) Simulating: |==================                                        |( 3/12) Simulating: |===================                                       |( 3/12) Simulating: |====================                                      |( 3/12) Simulating: |=====================                                     |( 3/12) Simulating: |======================                                    |( 3/12) Simulating: |=======================                                   |( 3/12) Simulating: |========================                                  |( 3/12) Simulating: |=========================                                 |( 3/12) Simulating: |==========================                                |( 3/12) Simulating: |===========================                               |( 3/12) Simulating: |============================                              |( 3/12) Simulating: |=============================                             |( 3/12) Simulating: |==============================                            |( 3/12) Simulating: |===============================                           |( 3/12) Simulating: |================================                          |( 3/12) Simulating: |=================================                         |( 3/12) Simulating: |==================================                        |( 3/12) Simulating: |===================================                       |( 3/12) Simulating: |====================================                      |( 3/12) Simulating: |=====================================                     |( 3/12) Simulating: |======================================                    |( 3/12) Simulating: |=======================================                   |( 3/12) Simulating: |========================================                  |( 3/12) Simulating: |=========================================                 |( 3/12) Simulating: |==========================================                |( 3/12) Simulating: |===========================================               |( 3/12) Simulating: |============================================              |( 3/12) Simulating: |=============================================             |( 3/12) Simulating: |==============================================            |( 3/12) Simulating: |===============================================           |( 3/12) Simulating: |================================================          |( 3/12) Simulating: |=================================================         |( 3/12) Simulating: |==================================================        |( 3/12) Simulating: |===================================================       |( 3/12) Simulating: |====================================================      |( 3/12) Simulating: |=====================================================     |( 3/12) Simulating: |======================================================    |( 3/12) Simulating: |=======================================================   |( 3/12) Simulating: |========================================================  |( 3/12) Simulating: |========================================================= |( 3/12) Simulating: |==========================================================|( 3/12) ( 4/12) ( 4/12) Simulating: |                                                          |( 4/12) Simulating: |=                                                         |( 4/12) Simulating: |==                                                        |( 4/12) Simulating: |===                                                       |( 4/12) Simulating: |====                                                      |( 4/12) Simulating: |=====                                                     |( 4/12) Simulating: |======                                                    |( 4/12) Simulating: |=======                                                   |( 4/12) Simulating: |========                                                  |( 4/12) Simulating: |=========                                                 |( 4/12) Simulating: |==========                                                |( 4/12) Simulating: |===========                                               |( 4/12) Simulating: |============                                              |( 4/12) Simulating: |=============                                             |( 4/12) Simulating: |==============                                            |( 4/12) Simulating: |===============                                           |( 4/12) Simulating: |================                                          |( 4/12) Simulating: |=================                                         |( 4/12) Simulating: |==================                                        |( 4/12) Simulating: |===================                                       |( 4/12) Simulating: |====================                                      |( 4/12) Simulating: |=====================                                     |( 4/12) Simulating: |======================                                    |( 4/12) Simulating: |=======================                                   |( 4/12) Simulating: |========================                                  |( 4/12) Simulating: |=========================                                 |( 4/12) Simulating: |==========================                                |( 4/12) Simulating: |===========================                               |( 4/12) Simulating: |============================                              |( 4/12) Simulating: |=============================                             |( 4/12) Simulating: |==============================                            |( 4/12) Simulating: |===============================                           |( 4/12) Simulating: |================================                          |( 4/12) Simulating: |=================================                         |( 4/12) Simulating: |==================================                        |( 4/12) Simulating: |===================================                       |( 4/12) Simulating: |====================================                      |( 4/12) Simulating: |=====================================                     |( 4/12) Simulating: |======================================                    |( 4/12) Simulating: |=======================================                   |( 4/12) Simulating: |========================================                  |( 4/12) Simulating: |=========================================                 |( 4/12) Simulating: |==========================================                |( 4/12) Simulating: |===========================================               |( 4/12) Simulating: |============================================              |( 4/12) Simulating: |=============================================             |( 4/12) Simulating: |==============================================            |( 4/12) Simulating: |===============================================           |( 4/12) Simulating: |================================================          |( 4/12) Simulating: |=================================================         |( 4/12) Simulating: |==================================================        |( 4/12) Simulating: |===================================================       |( 4/12) Simulating: |====================================================      |( 4/12) Simulating: |=====================================================     |( 4/12) Simulating: |======================================================    |( 4/12) Simulating: |=======================================================   |( 4/12) Simulating: |========================================================  |( 4/12) Simulating: |========================================================= |( 4/12) Simulating: |==========================================================|( 4/12) ( 5/12) ( 5/12) Simulating: |                                                          |( 5/12) Simulating: |=                                                         |( 5/12) Simulating: |==                                                        |( 5/12) Simulating: |===                                                       |( 5/12) Simulating: |====                                                      |( 5/12) Simulating: |=====                                                     |( 5/12) Simulating: |======                                                    |( 5/12) Simulating: |=======                                                   |( 5/12) Simulating: |========                                                  |( 5/12) Simulating: |=========                                                 |( 5/12) Simulating: |==========                                                |( 5/12) Simulating: |===========                                               |( 5/12) Simulating: |============                                              |( 5/12) Simulating: |=============                                             |( 5/12) Simulating: |==============                                            |( 5/12) Simulating: |===============                                           |( 5/12) Simulating: |================                                          |( 5/12) Simulating: |=================                                         |( 5/12) Simulating: |==================                                        |( 5/12) Simulating: |===================                                       |( 5/12) Simulating: |====================                                      |( 5/12) Simulating: |=====================                                     |( 5/12) Simulating: |======================                                    |( 5/12) Simulating: |=======================                                   |( 5/12) Simulating: |========================                                  |( 5/12) Simulating: |=========================                                 |( 5/12) Simulating: |==========================                                |( 5/12) Simulating: |===========================                               |( 5/12) Simulating: |============================                              |( 5/12) Simulating: |=============================                             |( 5/12) Simulating: |==============================                            |( 5/12) Simulating: |===============================                           |( 5/12) Simulating: |================================                          |( 5/12) Simulating: |=================================                         |( 5/12) Simulating: |==================================                        |( 5/12) Simulating: |===================================                       |( 5/12) Simulating: |====================================                      |( 5/12) Simulating: |=====================================                     |( 5/12) Simulating: |======================================                    |( 5/12) Simulating: |=======================================                   |( 5/12) Simulating: |========================================                  |( 5/12) Simulating: |=========================================                 |( 5/12) Simulating: |==========================================                |( 5/12) Simulating: |===========================================               |( 5/12) Simulating: |============================================              |( 5/12) Simulating: |=============================================             |( 5/12) Simulating: |==============================================            |( 5/12) Simulating: |===============================================           |( 5/12) Simulating: |================================================          |( 5/12) Simulating: |=================================================         |( 5/12) Simulating: |==================================================        |( 5/12) Simulating: |===================================================       |( 5/12) Simulating: |====================================================      |( 5/12) Simulating: |=====================================================     |( 5/12) Simulating: |======================================================    |( 5/12) Simulating: |=======================================================   |( 5/12) Simulating: |========================================================  |( 5/12) Simulating: |========================================================= |( 5/12) Simulating: |==========================================================|( 5/12) ( 6/12) ( 6/12) Simulating: |                                                          |( 6/12) Simulating: |=                                                         |( 6/12) Simulating: |==                                                        |( 6/12) Simulating: |===                                                       |( 6/12) Simulating: |====                                                      |( 6/12) Simulating: |=====                                                     |( 6/12) Simulating: |======                                                    |( 6/12) Simulating: |=======                                                   |( 6/12) Simulating: |========                                                  |( 6/12) Simulating: |=========                                                 |( 6/12) Simulating: |==========                                                |( 6/12) Simulating: |===========                                               |( 6/12) Simulating: |============                                              |( 6/12) Simulating: |=============                                             |( 6/12) Simulating: |==============                                            |( 6/12) Simulating: |===============                                           |( 6/12) Simulating: |================                                          |( 6/12) Simulating: |=================                                         |( 6/12) Simulating: |==================                                        |( 6/12) Simulating: |===================                                       |( 6/12) Simulating: |====================                                      |( 6/12) Simulating: |=====================                                     |( 6/12) Simulating: |======================                                    |( 6/12) Simulating: |=======================                                   |( 6/12) Simulating: |========================                                  |( 6/12) Simulating: |=========================                                 |( 6/12) Simulating: |==========================                                |( 6/12) Simulating: |===========================                               |( 6/12) Simulating: |============================                              |( 6/12) Simulating: |=============================                             |( 6/12) Simulating: |==============================                            |( 6/12) Simulating: |===============================                           |( 6/12) Simulating: |================================                          |( 6/12) Simulating: |=================================                         |( 6/12) Simulating: |==================================                        |( 6/12) Simulating: |===================================                       |( 6/12) Simulating: |====================================                      |( 6/12) Simulating: |=====================================                     |( 6/12) Simulating: |======================================                    |( 6/12) Simulating: |=======================================                   |( 6/12) Simulating: |========================================                  |( 6/12) Simulating: |=========================================                 |( 6/12) Simulating: |==========================================                |( 6/12) Simulating: |===========================================               |( 6/12) Simulating: |============================================              |( 6/12) Simulating: |=============================================             |( 6/12) Simulating: |==============================================            |( 6/12) Simulating: |===============================================           |( 6/12) Simulating: |================================================          |( 6/12) Simulating: |=================================================         |( 6/12) Simulating: |==================================================        |( 6/12) Simulating: |===================================================       |( 6/12) Simulating: |====================================================      |( 6/12) Simulating: |=====================================================     |( 6/12) Simulating: |======================================================    |( 6/12) Simulating: |=======================================================   |( 6/12) Simulating: |========================================================  |( 6/12) Simulating: |========================================================= |( 6/12) Simulating: |==========================================================|( 6/12) ( 7/12) ( 7/12) Simulating: |                                                          |( 7/12) Simulating: |=                                                         |( 7/12) Simulating: |==                                                        |( 7/12) Simulating: |===                                                       |( 7/12) Simulating: |====                                                      |( 7/12) Simulating: |=====                                                     |( 7/12) Simulating: |======                                                    |( 7/12) Simulating: |=======                                                   |( 7/12) Simulating: |========                                                  |( 7/12) Simulating: |=========                                                 |( 7/12) Simulating: |==========                                                |( 7/12) Simulating: |===========                                               |( 7/12) Simulating: |============                                              |( 7/12) Simulating: |=============                                             |( 7/12) Simulating: |==============                                            |( 7/12) Simulating: |===============                                           |( 7/12) Simulating: |================                                          |( 7/12) Simulating: |=================                                         |( 7/12) Simulating: |==================                                        |( 7/12) Simulating: |===================                                       |( 7/12) Simulating: |====================                                      |( 7/12) Simulating: |=====================                                     |( 7/12) Simulating: |======================                                    |( 7/12) Simulating: |=======================                                   |( 7/12) Simulating: |========================                                  |( 7/12) Simulating: |=========================                                 |( 7/12) Simulating: |==========================                                |( 7/12) Simulating: |===========================                               |( 7/12) Simulating: |============================                              |( 7/12) Simulating: |=============================                             |( 7/12) Simulating: |==============================                            |( 7/12) Simulating: |===============================                           |( 7/12) Simulating: |================================                          |( 7/12) Simulating: |=================================                         |( 7/12) Simulating: |==================================                        |( 7/12) Simulating: |===================================                       |( 7/12) Simulating: |====================================                      |( 7/12) Simulating: |=====================================                     |( 7/12) Simulating: |======================================                    |( 7/12) Simulating: |=======================================                   |( 7/12) Simulating: |========================================                  |( 7/12) Simulating: |=========================================                 |( 7/12) Simulating: |==========================================                |( 7/12) Simulating: |===========================================               |( 7/12) Simulating: |============================================              |( 7/12) Simulating: |=============================================             |( 7/12) Simulating: |==============================================            |( 7/12) Simulating: |===============================================           |( 7/12) Simulating: |================================================          |( 7/12) Simulating: |=================================================         |( 7/12) Simulating: |==================================================        |( 7/12) Simulating: |===================================================       |( 7/12) Simulating: |====================================================      |( 7/12) Simulating: |=====================================================     |( 7/12) Simulating: |======================================================    |( 7/12) Simulating: |=======================================================   |( 7/12) Simulating: |========================================================  |( 7/12) Simulating: |========================================================= |( 7/12) Simulating: |==========================================================|( 7/12) ( 8/12) ( 8/12) Simulating: |                                                          |( 8/12) Simulating: |=                                                         |( 8/12) Simulating: |==                                                        |( 8/12) Simulating: |===                                                       |( 8/12) Simulating: |====                                                      |( 8/12) Simulating: |=====                                                     |( 8/12) Simulating: |======                                                    |( 8/12) Simulating: |=======                                                   |( 8/12) Simulating: |========                                                  |( 8/12) Simulating: |=========                                                 |( 8/12) Simulating: |==========                                                |( 8/12) Simulating: |===========                                               |( 8/12) Simulating: |============                                              |( 8/12) Simulating: |=============                                             |( 8/12) Simulating: |==============                                            |( 8/12) Simulating: |===============                                           |( 8/12) Simulating: |================                                          |( 8/12) Simulating: |=================                                         |( 8/12) Simulating: |==================                                        |( 8/12) Simulating: |===================                                       |( 8/12) Simulating: |====================                                      |( 8/12) Simulating: |=====================                                     |( 8/12) Simulating: |======================                                    |( 8/12) Simulating: |=======================                                   |( 8/12) Simulating: |========================                                  |( 8/12) Simulating: |=========================                                 |( 8/12) Simulating: |==========================                                |( 8/12) Simulating: |===========================                               |( 8/12) Simulating: |============================                              |( 8/12) Simulating: |=============================                             |( 8/12) Simulating: |==============================                            |( 8/12) Simulating: |===============================                           |( 8/12) Simulating: |================================                          |( 8/12) Simulating: |=================================                         |( 8/12) Simulating: |==================================                        |( 8/12) Simulating: |===================================                       |( 8/12) Simulating: |====================================                      |( 8/12) Simulating: |=====================================                     |( 8/12) Simulating: |======================================                    |( 8/12) Simulating: |=======================================                   |( 8/12) Simulating: |========================================                  |( 8/12) Simulating: |=========================================                 |( 8/12) Simulating: |==========================================                |( 8/12) Simulating: |===========================================               |( 8/12) Simulating: |============================================              |( 8/12) Simulating: |=============================================             |( 8/12) Simulating: |==============================================            |( 8/12) Simulating: |===============================================           |( 8/12) Simulating: |================================================          |( 8/12) Simulating: |=================================================         |( 8/12) Simulating: |==================================================        |( 8/12) Simulating: |===================================================       |( 8/12) Simulating: |====================================================      |( 8/12) Simulating: |=====================================================     |( 8/12) Simulating: |======================================================    |( 8/12) Simulating: |=======================================================   |( 8/12) Simulating: |========================================================  |( 8/12) Simulating: |========================================================= |( 8/12) Simulating: |==========================================================|( 8/12) ( 9/12) ( 9/12) Simulating: |                                                          |( 9/12) Simulating: |=                                                         |( 9/12) Simulating: |==                                                        |( 9/12) Simulating: |===                                                       |( 9/12) Simulating: |====                                                      |( 9/12) Simulating: |=====                                                     |( 9/12) Simulating: |======                                                    |( 9/12) Simulating: |=======                                                   |( 9/12) Simulating: |========                                                  |( 9/12) Simulating: |=========                                                 |( 9/12) Simulating: |==========                                                |( 9/12) Simulating: |===========                                               |( 9/12) Simulating: |============                                              |( 9/12) Simulating: |=============                                             |( 9/12) Simulating: |==============                                            |( 9/12) Simulating: |===============                                           |( 9/12) Simulating: |================                                          |( 9/12) Simulating: |=================                                         |( 9/12) Simulating: |==================                                        |( 9/12) Simulating: |===================                                       |( 9/12) Simulating: |====================                                      |( 9/12) Simulating: |=====================                                     |( 9/12) Simulating: |======================                                    |( 9/12) Simulating: |=======================                                   |( 9/12) Simulating: |========================                                  |( 9/12) Simulating: |=========================                                 |( 9/12) Simulating: |==========================                                |( 9/12) Simulating: |===========================                               |( 9/12) Simulating: |============================                              |( 9/12) Simulating: |=============================                             |( 9/12) Simulating: |==============================                            |( 9/12) Simulating: |===============================                           |( 9/12) Simulating: |================================                          |( 9/12) Simulating: |=================================                         |( 9/12) Simulating: |==================================                        |( 9/12) Simulating: |===================================                       |( 9/12) Simulating: |====================================                      |( 9/12) Simulating: |=====================================                     |( 9/12) Simulating: |======================================                    |( 9/12) Simulating: |=======================================                   |( 9/12) Simulating: |========================================                  |( 9/12) Simulating: |=========================================                 |( 9/12) Simulating: |==========================================                |( 9/12) Simulating: |===========================================               |( 9/12) Simulating: |============================================              |( 9/12) Simulating: |=============================================             |( 9/12) Simulating: |==============================================            |( 9/12) Simulating: |===============================================           |( 9/12) Simulating: |================================================          |( 9/12) Simulating: |=================================================         |( 9/12) Simulating: |==================================================        |( 9/12) Simulating: |===================================================       |( 9/12) Simulating: |====================================================      |( 9/12) Simulating: |=====================================================     |( 9/12) Simulating: |======================================================    |( 9/12) Simulating: |=======================================================   |( 9/12) Simulating: |========================================================  |( 9/12) Simulating: |========================================================= |( 9/12) Simulating: |==========================================================|( 9/12) (10/12) (10/12) Simulating: |                                                          |(10/12) Simulating: |=                                                         |(10/12) Simulating: |==                                                        |(10/12) Simulating: |===                                                       |(10/12) Simulating: |====                                                      |(10/12) Simulating: |=====                                                     |(10/12) Simulating: |======                                                    |(10/12) Simulating: |=======                                                   |(10/12) Simulating: |========                                                  |(10/12) Simulating: |=========                                                 |(10/12) Simulating: |==========                                                |(10/12) Simulating: |===========                                               |(10/12) Simulating: |============                                              |(10/12) Simulating: |=============                                             |(10/12) Simulating: |==============                                            |(10/12) Simulating: |===============                                           |(10/12) Simulating: |================                                          |(10/12) Simulating: |=================                                         |(10/12) Simulating: |==================                                        |(10/12) Simulating: |===================                                       |(10/12) Simulating: |====================                                      |(10/12) Simulating: |=====================                                     |(10/12) Simulating: |======================                                    |(10/12) Simulating: |=======================                                   |(10/12) Simulating: |========================                                  |(10/12) Simulating: |=========================                                 |(10/12) Simulating: |==========================                                |(10/12) Simulating: |===========================                               |(10/12) Simulating: |============================                              |(10/12) Simulating: |=============================                             |(10/12) Simulating: |==============================                            |(10/12) Simulating: |===============================                           |(10/12) Simulating: |================================                          |(10/12) Simulating: |=================================                         |(10/12) Simulating: |==================================                        |(10/12) Simulating: |===================================                       |(10/12) Simulating: |====================================                      |(10/12) Simulating: |=====================================                     |(10/12) Simulating: |======================================                    |(10/12) Simulating: |=======================================                   |(10/12) Simulating: |========================================                  |(10/12) Simulating: |=========================================                 |(10/12) Simulating: |==========================================                |(10/12) Simulating: |===========================================               |(10/12) Simulating: |============================================              |(10/12) Simulating: |=============================================             |(10/12) Simulating: |==============================================            |(10/12) Simulating: |===============================================           |(10/12) Simulating: |================================================          |(10/12) Simulating: |=================================================         |(10/12) Simulating: |==================================================        |(10/12) Simulating: |===================================================       |(10/12) Simulating: |====================================================      |(10/12) Simulating: |=====================================================     |(10/12) Simulating: |======================================================    |(10/12) Simulating: |=======================================================   |(10/12) Simulating: |========================================================  |(10/12) Simulating: |========================================================= |(10/12) Simulating: |==========================================================|(10/12) (11/12) (11/12) Simulating: |                                                          |(11/12) Simulating: |=                                                         |(11/12) Simulating: |==                                                        |(11/12) Simulating: |===                                                       |(11/12) Simulating: |====                                                      |(11/12) Simulating: |=====                                                     |(11/12) Simulating: |======                                                    |(11/12) Simulating: |=======                                                   |(11/12) Simulating: |========                                                  |(11/12) Simulating: |=========                                                 |(11/12) Simulating: |==========                                                |(11/12) Simulating: |===========                                               |(11/12) Simulating: |============                                              |(11/12) Simulating: |=============                                             |(11/12) Simulating: |==============                                            |(11/12) Simulating: |===============                                           |(11/12) Simulating: |================                                          |(11/12) Simulating: |=================                                         |(11/12) Simulating: |==================                                        |(11/12) Simulating: |===================                                       |(11/12) Simulating: |====================                                      |(11/12) Simulating: |=====================                                     |(11/12) Simulating: |======================                                    |(11/12) Simulating: |=======================                                   |(11/12) Simulating: |========================                                  |(11/12) Simulating: |=========================                                 |(11/12) Simulating: |==========================                                |(11/12) Simulating: |===========================                               |(11/12) Simulating: |============================                              |(11/12) Simulating: |=============================                             |(11/12) Simulating: |==============================                            |(11/12) Simulating: |===============================                           |(11/12) Simulating: |================================                          |(11/12) Simulating: |=================================                         |(11/12) Simulating: |==================================                        |(11/12) Simulating: |===================================                       |(11/12) Simulating: |====================================                      |(11/12) Simulating: |=====================================                     |(11/12) Simulating: |======================================                    |(11/12) Simulating: |=======================================                   |(11/12) Simulating: |========================================                  |(11/12) Simulating: |=========================================                 |(11/12) Simulating: |==========================================                |(11/12) Simulating: |===========================================               |(11/12) Simulating: |============================================              |(11/12) Simulating: |=============================================             |(11/12) Simulating: |==============================================            |(11/12) Simulating: |===============================================           |(11/12) Simulating: |================================================          |(11/12) Simulating: |=================================================         |(11/12) Simulating: |==================================================        |(11/12) Simulating: |===================================================       |(11/12) Simulating: |====================================================      |(11/12) Simulating: |=====================================================     |(11/12) Simulating: |======================================================    |(11/12) Simulating: |=======================================================   |(11/12) Simulating: |========================================================  |(11/12) Simulating: |========================================================= |(11/12) Simulating: |==========================================================|(11/12) (12/12) (12/12) Simulating: |                                                          |(12/12) Simulating: |=                                                         |(12/12) Simulating: |==                                                        |(12/12) Simulating: |===                                                       |(12/12) Simulating: |====                                                      |(12/12) Simulating: |=====                                                     |(12/12) Simulating: |======                                                    |(12/12) Simulating: |=======                                                   |(12/12) Simulating: |========                                                  |(12/12) Simulating: |=========                                                 |(12/12) Simulating: |==========                                                |(12/12) Simulating: |===========                                               |(12/12) Simulating: |============                                              |(12/12) Simulating: |=============                                             |(12/12) Simulating: |==============                                            |(12/12) Simulating: |===============                                           |(12/12) Simulating: |================                                          |(12/12) Simulating: |=================                                         |(12/12) Simulating: |==================                                        |(12/12) Simulating: |===================                                       |(12/12) Simulating: |====================                                      |(12/12) Simulating: |=====================                                     |(12/12) Simulating: |======================                                    |(12/12) Simulating: |=======================                                   |(12/12) Simulating: |========================                                  |(12/12) Simulating: |=========================                                 |(12/12) Simulating: |==========================                                |(12/12) Simulating: |===========================                               |(12/12) Simulating: |============================                              |(12/12) Simulating: |=============================                             |(12/12) Simulating: |==============================                            |(12/12) Simulating: |===============================                           |(12/12) Simulating: |================================                          |(12/12) Simulating: |=================================                         |(12/12) Simulating: |==================================                        |(12/12) Simulating: |===================================                       |(12/12) Simulating: |====================================                      |(12/12) Simulating: |=====================================                     |(12/12) Simulating: |======================================                    |(12/12) Simulating: |=======================================                   |(12/12) Simulating: |========================================                  |(12/12) Simulating: |=========================================                 |(12/12) Simulating: |==========================================                |(12/12) Simulating: |===========================================               |(12/12) Simulating: |============================================              |(12/12) Simulating: |=============================================             |(12/12) Simulating: |==============================================            |(12/12) Simulating: |===============================================           |(12/12) Simulating: |================================================          |(12/12) Simulating: |=================================================         |(12/12) Simulating: |==================================================        |(12/12) Simulating: |===================================================       |(12/12) Simulating: |====================================================      |(12/12) Simulating: |=====================================================     |(12/12) Simulating: |======================================================    |(12/12) Simulating: |=======================================================   |(12/12) Simulating: |========================================================  |(12/12) Simulating: |========================================================= |(12/12) Simulating: |==========================================================|(12/12) 
curve_fine_grained <- powerCurve(model_sim,
                    fixed("RDM_typefinegrained_mind_action_physiology_model", "t"),
                    along = "subject_id",
                    breaks = c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
                    nsim = 100)
## Simulating: |                                                                  |Simulating: |=                                                                 |Simulating: |==                                                                |Simulating: |===                                                               |Simulating: |====                                                              |Simulating: |=====                                                             |Simulating: |======                                                            |Simulating: |=======                                                           |Simulating: |========                                                          |Simulating: |=========                                                         |Simulating: |==========                                                        |Simulating: |===========                                                       |Simulating: |============                                                      |Simulating: |=============                                                     |Simulating: |==============                                                    |Simulating: |===============                                                   |Simulating: |================                                                  |Simulating: |=================                                                 |Simulating: |==================                                                |Simulating: |===================                                               |Simulating: |====================                                              |Simulating: |=====================                                             |Simulating: |======================                                            |Simulating: |=======================                                           |Simulating: |========================                                          |Simulating: |=========================                                         |Simulating: |==========================                                        |Simulating: |===========================                                       |Simulating: |============================                                      |Simulating: |=============================                                     |Simulating: |==============================                                    |Simulating: |===============================                                   |Simulating: |================================                                  |Simulating: |=================================                                 |Simulating: |==================================                                |Simulating: |===================================                               |Simulating: |====================================                              |Simulating: |=====================================                             |Simulating: |======================================                            |Simulating: |=======================================                           |Simulating: |========================================                          |Simulating: |=========================================                         |Simulating: |==========================================                        |Simulating: |===========================================                       |Simulating: |============================================                      |Simulating: |=============================================                     |Simulating: |==============================================                    |Simulating: |===============================================                   |Simulating: |================================================                  |Simulating: |=================================================                 |Simulating: |==================================================                |Simulating: |===================================================               |Simulating: |====================================================              |Simulating: |=====================================================             |Simulating: |======================================================            |Simulating: |=======================================================           |Simulating: |========================================================          |Simulating: |=========================================================         |Simulating: |==========================================================        |Simulating: |===========================================================       |Simulating: |============================================================      |Simulating: |=============================================================     |Simulating: |==============================================================    |Simulating: |===============================================================   |Simulating: |================================================================  |Simulating: |================================================================= |Simulating: |==================================================================|( 1/10) ( 1/10) Simulating: |                                                          |( 1/10) Simulating: |=                                                         |( 1/10) Simulating: |==                                                        |( 1/10) Simulating: |===                                                       |( 1/10) Simulating: |====                                                      |( 1/10) Simulating: |=====                                                     |( 1/10) Simulating: |======                                                    |( 1/10) Simulating: |=======                                                   |( 1/10) Simulating: |========                                                  |( 1/10) Simulating: |=========                                                 |( 1/10) Simulating: |==========                                                |( 1/10) Simulating: |===========                                               |( 1/10) Simulating: |============                                              |( 1/10) Simulating: |=============                                             |( 1/10) Simulating: |==============                                            |( 1/10) Simulating: |===============                                           |( 1/10) Simulating: |================                                          |( 1/10) Simulating: |=================                                         |( 1/10) Simulating: |==================                                        |( 1/10) Simulating: |===================                                       |( 1/10) Simulating: |====================                                      |( 1/10) Simulating: |=====================                                     |( 1/10) Simulating: |======================                                    |( 1/10) Simulating: |=======================                                   |( 1/10) Simulating: |========================                                  |( 1/10) Simulating: |=========================                                 |( 1/10) Simulating: |==========================                                |( 1/10) Simulating: |===========================                               |( 1/10) Simulating: |============================                              |( 1/10) Simulating: |=============================                             |( 1/10) Simulating: |==============================                            |( 1/10) Simulating: |===============================                           |( 1/10) Simulating: |================================                          |( 1/10) Simulating: |=================================                         |( 1/10) Simulating: |==================================                        |( 1/10) Simulating: |===================================                       |( 1/10) Simulating: |====================================                      |( 1/10) Simulating: |=====================================                     |( 1/10) Simulating: |======================================                    |( 1/10) Simulating: |=======================================                   |( 1/10) Simulating: |========================================                  |( 1/10) Simulating: |=========================================                 |( 1/10) Simulating: |==========================================                |( 1/10) Simulating: |===========================================               |( 1/10) Simulating: |============================================              |( 1/10) Simulating: |=============================================             |( 1/10) Simulating: |==============================================            |( 1/10) Simulating: |===============================================           |( 1/10) Simulating: |================================================          |( 1/10) Simulating: |=================================================         |( 1/10) Simulating: |==================================================        |( 1/10) Simulating: |===================================================       |( 1/10) Simulating: |====================================================      |( 1/10) Simulating: |=====================================================     |( 1/10) Simulating: |======================================================    |( 1/10) Simulating: |=======================================================   |( 1/10) Simulating: |========================================================  |( 1/10) Simulating: |========================================================= |( 1/10) Simulating: |==========================================================|( 1/10) ( 2/10) ( 2/10) Simulating: |                                                          |( 2/10) Simulating: |=                                                         |( 2/10) Simulating: |==                                                        |( 2/10) Simulating: |===                                                       |( 2/10) Simulating: |====                                                      |( 2/10) Simulating: |=====                                                     |( 2/10) Simulating: |======                                                    |( 2/10) Simulating: |=======                                                   |( 2/10) Simulating: |========                                                  |( 2/10) Simulating: |=========                                                 |( 2/10) Simulating: |==========                                                |( 2/10) Simulating: |===========                                               |( 2/10) Simulating: |============                                              |( 2/10) Simulating: |=============                                             |( 2/10) Simulating: |==============                                            |( 2/10) Simulating: |===============                                           |( 2/10) Simulating: |================                                          |( 2/10) Simulating: |=================                                         |( 2/10) Simulating: |==================                                        |( 2/10) Simulating: |===================                                       |( 2/10) Simulating: |====================                                      |( 2/10) Simulating: |=====================                                     |( 2/10) Simulating: |======================                                    |( 2/10) Simulating: |=======================                                   |( 2/10) Simulating: |========================                                  |( 2/10) Simulating: |=========================                                 |( 2/10) Simulating: |==========================                                |( 2/10) Simulating: |===========================                               |( 2/10) Simulating: |============================                              |( 2/10) Simulating: |=============================                             |( 2/10) Simulating: |==============================                            |( 2/10) Simulating: |===============================                           |( 2/10) Simulating: |================================                          |( 2/10) Simulating: |=================================                         |( 2/10) Simulating: |==================================                        |( 2/10) Simulating: |===================================                       |( 2/10) Simulating: |====================================                      |( 2/10) Simulating: |=====================================                     |( 2/10) Simulating: |======================================                    |( 2/10) Simulating: |=======================================                   |( 2/10) Simulating: |========================================                  |( 2/10) Simulating: |=========================================                 |( 2/10) Simulating: |==========================================                |( 2/10) Simulating: |===========================================               |( 2/10) Simulating: |============================================              |( 2/10) Simulating: |=============================================             |( 2/10) Simulating: |==============================================            |( 2/10) Simulating: |===============================================           |( 2/10) Simulating: |================================================          |( 2/10) Simulating: |=================================================         |( 2/10) Simulating: |==================================================        |( 2/10) Simulating: |===================================================       |( 2/10) Simulating: |====================================================      |( 2/10) Simulating: |=====================================================     |( 2/10) Simulating: |======================================================    |( 2/10) Simulating: |=======================================================   |( 2/10) Simulating: |========================================================  |( 2/10) Simulating: |========================================================= |( 2/10) Simulating: |==========================================================|( 2/10) ( 3/10) ( 3/10) Simulating: |                                                          |( 3/10) Simulating: |=                                                         |( 3/10) Simulating: |==                                                        |( 3/10) Simulating: |===                                                       |( 3/10) Simulating: |====                                                      |( 3/10) Simulating: |=====                                                     |( 3/10) Simulating: |======                                                    |( 3/10) Simulating: |=======                                                   |( 3/10) Simulating: |========                                                  |( 3/10) Simulating: |=========                                                 |( 3/10) Simulating: |==========                                                |( 3/10) Simulating: |===========                                               |( 3/10) Simulating: |============                                              |( 3/10) Simulating: |=============                                             |( 3/10) Simulating: |==============                                            |( 3/10) Simulating: |===============                                           |( 3/10) Simulating: |================                                          |( 3/10) Simulating: |=================                                         |( 3/10) Simulating: |==================                                        |( 3/10) Simulating: |===================                                       |( 3/10) Simulating: |====================                                      |( 3/10) Simulating: |=====================                                     |( 3/10) Simulating: |======================                                    |( 3/10) Simulating: |=======================                                   |( 3/10) Simulating: |========================                                  |( 3/10) Simulating: |=========================                                 |( 3/10) Simulating: |==========================                                |( 3/10) Simulating: |===========================                               |( 3/10) Simulating: |============================                              |( 3/10) Simulating: |=============================                             |( 3/10) Simulating: |==============================                            |( 3/10) Simulating: |===============================                           |( 3/10) Simulating: |================================                          |( 3/10) Simulating: |=================================                         |( 3/10) Simulating: |==================================                        |( 3/10) Simulating: |===================================                       |( 3/10) Simulating: |====================================                      |( 3/10) Simulating: |=====================================                     |( 3/10) Simulating: |======================================                    |( 3/10) Simulating: |=======================================                   |( 3/10) Simulating: |========================================                  |( 3/10) Simulating: |=========================================                 |( 3/10) Simulating: |==========================================                |( 3/10) Simulating: |===========================================               |( 3/10) Simulating: |============================================              |( 3/10) Simulating: |=============================================             |( 3/10) Simulating: |==============================================            |( 3/10) Simulating: |===============================================           |( 3/10) Simulating: |================================================          |( 3/10) Simulating: |=================================================         |( 3/10) Simulating: |==================================================        |( 3/10) Simulating: |===================================================       |( 3/10) Simulating: |====================================================      |( 3/10) Simulating: |=====================================================     |( 3/10) Simulating: |======================================================    |( 3/10) Simulating: |=======================================================   |( 3/10) Simulating: |========================================================  |( 3/10) Simulating: |========================================================= |( 3/10) Simulating: |==========================================================|( 3/10) ( 4/10) ( 4/10) Simulating: |                                                          |( 4/10) Simulating: |=                                                         |( 4/10) Simulating: |==                                                        |( 4/10) Simulating: |===                                                       |( 4/10) Simulating: |====                                                      |( 4/10) Simulating: |=====                                                     |( 4/10) Simulating: |======                                                    |( 4/10) Simulating: |=======                                                   |( 4/10) Simulating: |========                                                  |( 4/10) Simulating: |=========                                                 |( 4/10) Simulating: |==========                                                |( 4/10) Simulating: |===========                                               |( 4/10) Simulating: |============                                              |( 4/10) Simulating: |=============                                             |( 4/10) Simulating: |==============                                            |( 4/10) Simulating: |===============                                           |( 4/10) Simulating: |================                                          |( 4/10) Simulating: |=================                                         |( 4/10) Simulating: |==================                                        |( 4/10) Simulating: |===================                                       |( 4/10) Simulating: |====================                                      |( 4/10) Simulating: |=====================                                     |( 4/10) Simulating: |======================                                    |( 4/10) Simulating: |=======================                                   |( 4/10) Simulating: |========================                                  |( 4/10) Simulating: |=========================                                 |( 4/10) Simulating: |==========================                                |( 4/10) Simulating: |===========================                               |( 4/10) Simulating: |============================                              |( 4/10) Simulating: |=============================                             |( 4/10) Simulating: |==============================                            |( 4/10) Simulating: |===============================                           |( 4/10) Simulating: |================================                          |( 4/10) Simulating: |=================================                         |( 4/10) Simulating: |==================================                        |( 4/10) Simulating: |===================================                       |( 4/10) Simulating: |====================================                      |( 4/10) Simulating: |=====================================                     |( 4/10) Simulating: |======================================                    |( 4/10) Simulating: |=======================================                   |( 4/10) Simulating: |========================================                  |( 4/10) Simulating: |=========================================                 |( 4/10) Simulating: |==========================================                |( 4/10) Simulating: |===========================================               |( 4/10) Simulating: |============================================              |( 4/10) Simulating: |=============================================             |( 4/10) Simulating: |==============================================            |( 4/10) Simulating: |===============================================           |( 4/10) Simulating: |================================================          |( 4/10) Simulating: |=================================================         |( 4/10) Simulating: |==================================================        |( 4/10) Simulating: |===================================================       |( 4/10) Simulating: |====================================================      |( 4/10) Simulating: |=====================================================     |( 4/10) Simulating: |======================================================    |( 4/10) Simulating: |=======================================================   |( 4/10) Simulating: |========================================================  |( 4/10) Simulating: |========================================================= |( 4/10) Simulating: |==========================================================|( 4/10) ( 5/10) ( 5/10) Simulating: |                                                          |( 5/10) Simulating: |=                                                         |( 5/10) Simulating: |==                                                        |( 5/10) Simulating: |===                                                       |( 5/10) Simulating: |====                                                      |( 5/10) Simulating: |=====                                                     |( 5/10) Simulating: |======                                                    |( 5/10) Simulating: |=======                                                   |( 5/10) Simulating: |========                                                  |( 5/10) Simulating: |=========                                                 |( 5/10) Simulating: |==========                                                |( 5/10) Simulating: |===========                                               |( 5/10) Simulating: |============                                              |( 5/10) Simulating: |=============                                             |( 5/10) Simulating: |==============                                            |( 5/10) Simulating: |===============                                           |( 5/10) Simulating: |================                                          |( 5/10) Simulating: |=================                                         |( 5/10) Simulating: |==================                                        |( 5/10) Simulating: |===================                                       |( 5/10) Simulating: |====================                                      |( 5/10) Simulating: |=====================                                     |( 5/10) Simulating: |======================                                    |( 5/10) Simulating: |=======================                                   |( 5/10) Simulating: |========================                                  |( 5/10) Simulating: |=========================                                 |( 5/10) Simulating: |==========================                                |( 5/10) Simulating: |===========================                               |( 5/10) Simulating: |============================                              |( 5/10) Simulating: |=============================                             |( 5/10) Simulating: |==============================                            |( 5/10) Simulating: |===============================                           |( 5/10) Simulating: |================================                          |( 5/10) Simulating: |=================================                         |( 5/10) Simulating: |==================================                        |( 5/10) Simulating: |===================================                       |( 5/10) Simulating: |====================================                      |( 5/10) Simulating: |=====================================                     |( 5/10) Simulating: |======================================                    |( 5/10) Simulating: |=======================================                   |( 5/10) Simulating: |========================================                  |( 5/10) Simulating: |=========================================                 |( 5/10) Simulating: |==========================================                |( 5/10) Simulating: |===========================================               |( 5/10) Simulating: |============================================              |( 5/10) Simulating: |=============================================             |( 5/10) Simulating: |==============================================            |( 5/10) Simulating: |===============================================           |( 5/10) Simulating: |================================================          |( 5/10) Simulating: |=================================================         |( 5/10) Simulating: |==================================================        |( 5/10) Simulating: |===================================================       |( 5/10) Simulating: |====================================================      |( 5/10) Simulating: |=====================================================     |( 5/10) Simulating: |======================================================    |( 5/10) Simulating: |=======================================================   |( 5/10) Simulating: |========================================================  |( 5/10) Simulating: |========================================================= |( 5/10) Simulating: |==========================================================|( 5/10) ( 6/10) ( 6/10) Simulating: |                                                          |( 6/10) Simulating: |=                                                         |( 6/10) Simulating: |==                                                        |( 6/10) Simulating: |===                                                       |( 6/10) Simulating: |====                                                      |( 6/10) Simulating: |=====                                                     |( 6/10) Simulating: |======                                                    |( 6/10) Simulating: |=======                                                   |( 6/10) Simulating: |========                                                  |( 6/10) Simulating: |=========                                                 |( 6/10) Simulating: |==========                                                |( 6/10) Simulating: |===========                                               |( 6/10) Simulating: |============                                              |( 6/10) Simulating: |=============                                             |( 6/10) Simulating: |==============                                            |( 6/10) Simulating: |===============                                           |( 6/10) Simulating: |================                                          |( 6/10) Simulating: |=================                                         |( 6/10) Simulating: |==================                                        |( 6/10) Simulating: |===================                                       |( 6/10) Simulating: |====================                                      |( 6/10) Simulating: |=====================                                     |( 6/10) Simulating: |======================                                    |( 6/10) Simulating: |=======================                                   |( 6/10) Simulating: |========================                                  |( 6/10) Simulating: |=========================                                 |( 6/10) Simulating: |==========================                                |( 6/10) Simulating: |===========================                               |( 6/10) Simulating: |============================                              |( 6/10) Simulating: |=============================                             |( 6/10) Simulating: |==============================                            |( 6/10) Simulating: |===============================                           |( 6/10) Simulating: |================================                          |( 6/10) Simulating: |=================================                         |( 6/10) Simulating: |==================================                        |( 6/10) Simulating: |===================================                       |( 6/10) Simulating: |====================================                      |( 6/10) Simulating: |=====================================                     |( 6/10) Simulating: |======================================                    |( 6/10) Simulating: |=======================================                   |( 6/10) Simulating: |========================================                  |( 6/10) Simulating: |=========================================                 |( 6/10) Simulating: |==========================================                |( 6/10) Simulating: |===========================================               |( 6/10) Simulating: |============================================              |( 6/10) Simulating: |=============================================             |( 6/10) Simulating: |==============================================            |( 6/10) Simulating: |===============================================           |( 6/10) Simulating: |================================================          |( 6/10) Simulating: |=================================================         |( 6/10) Simulating: |==================================================        |( 6/10) Simulating: |===================================================       |( 6/10) Simulating: |====================================================      |( 6/10) Simulating: |=====================================================     |( 6/10) Simulating: |======================================================    |( 6/10) Simulating: |=======================================================   |( 6/10) Simulating: |========================================================  |( 6/10) Simulating: |========================================================= |( 6/10) Simulating: |==========================================================|( 6/10) ( 7/10) ( 7/10) Simulating: |                                                          |( 7/10) Simulating: |=                                                         |( 7/10) Simulating: |==                                                        |( 7/10) Simulating: |===                                                       |( 7/10) Simulating: |====                                                      |( 7/10) Simulating: |=====                                                     |( 7/10) Simulating: |======                                                    |( 7/10) Simulating: |=======                                                   |( 7/10) Simulating: |========                                                  |( 7/10) Simulating: |=========                                                 |( 7/10) Simulating: |==========                                                |( 7/10) Simulating: |===========                                               |( 7/10) Simulating: |============                                              |( 7/10) Simulating: |=============                                             |( 7/10) Simulating: |==============                                            |( 7/10) Simulating: |===============                                           |( 7/10) Simulating: |================                                          |( 7/10) Simulating: |=================                                         |( 7/10) Simulating: |==================                                        |( 7/10) Simulating: |===================                                       |( 7/10) Simulating: |====================                                      |( 7/10) Simulating: |=====================                                     |( 7/10) Simulating: |======================                                    |( 7/10) Simulating: |=======================                                   |( 7/10) Simulating: |========================                                  |( 7/10) Simulating: |=========================                                 |( 7/10) Simulating: |==========================                                |( 7/10) Simulating: |===========================                               |( 7/10) Simulating: |============================                              |( 7/10) Simulating: |=============================                             |( 7/10) Simulating: |==============================                            |( 7/10) Simulating: |===============================                           |( 7/10) Simulating: |================================                          |( 7/10) Simulating: |=================================                         |( 7/10) Simulating: |==================================                        |( 7/10) Simulating: |===================================                       |( 7/10) Simulating: |====================================                      |( 7/10) Simulating: |=====================================                     |( 7/10) Simulating: |======================================                    |( 7/10) Simulating: |=======================================                   |( 7/10) Simulating: |========================================                  |( 7/10) Simulating: |=========================================                 |( 7/10) Simulating: |==========================================                |( 7/10) Simulating: |===========================================               |( 7/10) Simulating: |============================================              |( 7/10) Simulating: |=============================================             |( 7/10) Simulating: |==============================================            |( 7/10) Simulating: |===============================================           |( 7/10) Simulating: |================================================          |( 7/10) Simulating: |=================================================         |( 7/10) Simulating: |==================================================        |( 7/10) Simulating: |===================================================       |( 7/10) Simulating: |====================================================      |( 7/10) Simulating: |=====================================================     |( 7/10) Simulating: |======================================================    |( 7/10) Simulating: |=======================================================   |( 7/10) Simulating: |========================================================  |( 7/10) Simulating: |========================================================= |( 7/10) Simulating: |==========================================================|( 7/10) ( 8/10) ( 8/10) Simulating: |                                                          |( 8/10) Simulating: |=                                                         |( 8/10) Simulating: |==                                                        |( 8/10) Simulating: |===                                                       |( 8/10) Simulating: |====                                                      |( 8/10) Simulating: |=====                                                     |( 8/10) Simulating: |======                                                    |( 8/10) Simulating: |=======                                                   |( 8/10) Simulating: |========                                                  |( 8/10) Simulating: |=========                                                 |( 8/10) Simulating: |==========                                                |( 8/10) Simulating: |===========                                               |( 8/10) Simulating: |============                                              |( 8/10) Simulating: |=============                                             |( 8/10) Simulating: |==============                                            |( 8/10) Simulating: |===============                                           |( 8/10) Simulating: |================                                          |( 8/10) Simulating: |=================                                         |( 8/10) Simulating: |==================                                        |( 8/10) Simulating: |===================                                       |( 8/10) Simulating: |====================                                      |( 8/10) Simulating: |=====================                                     |( 8/10) Simulating: |======================                                    |( 8/10) Simulating: |=======================                                   |( 8/10) Simulating: |========================                                  |( 8/10) Simulating: |=========================                                 |( 8/10) Simulating: |==========================                                |( 8/10) Simulating: |===========================                               |( 8/10) Simulating: |============================                              |( 8/10) Simulating: |=============================                             |( 8/10) Simulating: |==============================                            |( 8/10) Simulating: |===============================                           |( 8/10) Simulating: |================================                          |( 8/10) Simulating: |=================================                         |( 8/10) Simulating: |==================================                        |( 8/10) Simulating: |===================================                       |( 8/10) Simulating: |====================================                      |( 8/10) Simulating: |=====================================                     |( 8/10) Simulating: |======================================                    |( 8/10) Simulating: |=======================================                   |( 8/10) Simulating: |========================================                  |( 8/10) Simulating: |=========================================                 |( 8/10) Simulating: |==========================================                |( 8/10) Simulating: |===========================================               |( 8/10) Simulating: |============================================              |( 8/10) Simulating: |=============================================             |( 8/10) Simulating: |==============================================            |( 8/10) Simulating: |===============================================           |( 8/10) Simulating: |================================================          |( 8/10) Simulating: |=================================================         |( 8/10) Simulating: |==================================================        |( 8/10) Simulating: |===================================================       |( 8/10) Simulating: |====================================================      |( 8/10) Simulating: |=====================================================     |( 8/10) Simulating: |======================================================    |( 8/10) Simulating: |=======================================================   |( 8/10) Simulating: |========================================================  |( 8/10) Simulating: |========================================================= |( 8/10) Simulating: |==========================================================|( 8/10) ( 9/10) ( 9/10) Simulating: |                                                          |( 9/10) Simulating: |=                                                         |( 9/10) Simulating: |==                                                        |( 9/10) Simulating: |===                                                       |( 9/10) Simulating: |====                                                      |( 9/10) Simulating: |=====                                                     |( 9/10) Simulating: |======                                                    |( 9/10) Simulating: |=======                                                   |( 9/10) Simulating: |========                                                  |( 9/10) Simulating: |=========                                                 |( 9/10) Simulating: |==========                                                |( 9/10) Simulating: |===========                                               |( 9/10) Simulating: |============                                              |( 9/10) Simulating: |=============                                             |( 9/10) Simulating: |==============                                            |( 9/10) Simulating: |===============                                           |( 9/10) Simulating: |================                                          |( 9/10) Simulating: |=================                                         |( 9/10) Simulating: |==================                                        |( 9/10) Simulating: |===================                                       |( 9/10) Simulating: |====================                                      |( 9/10) Simulating: |=====================                                     |( 9/10) Simulating: |======================                                    |( 9/10) Simulating: |=======================                                   |( 9/10) Simulating: |========================                                  |( 9/10) Simulating: |=========================                                 |( 9/10) Simulating: |==========================                                |( 9/10) Simulating: |===========================                               |( 9/10) Simulating: |============================                              |( 9/10) Simulating: |=============================                             |( 9/10) Simulating: |==============================                            |( 9/10) Simulating: |===============================                           |( 9/10) Simulating: |================================                          |( 9/10) Simulating: |=================================                         |( 9/10) Simulating: |==================================                        |( 9/10) Simulating: |===================================                       |( 9/10) Simulating: |====================================                      |( 9/10) Simulating: |=====================================                     |( 9/10) Simulating: |======================================                    |( 9/10) Simulating: |=======================================                   |( 9/10) Simulating: |========================================                  |( 9/10) Simulating: |=========================================                 |( 9/10) Simulating: |==========================================                |( 9/10) Simulating: |===========================================               |( 9/10) Simulating: |============================================              |( 9/10) Simulating: |=============================================             |( 9/10) Simulating: |==============================================            |( 9/10) Simulating: |===============================================           |( 9/10) Simulating: |================================================          |( 9/10) Simulating: |=================================================         |( 9/10) Simulating: |==================================================        |( 9/10) Simulating: |===================================================       |( 9/10) Simulating: |====================================================      |( 9/10) Simulating: |=====================================================     |( 9/10) Simulating: |======================================================    |( 9/10) Simulating: |=======================================================   |( 9/10) Simulating: |========================================================  |( 9/10) Simulating: |========================================================= |( 9/10) Simulating: |==========================================================|( 9/10) (10/10) (10/10) Simulating: |                                                          |(10/10) Simulating: |=                                                         |(10/10) Simulating: |==                                                        |(10/10) Simulating: |===                                                       |(10/10) Simulating: |====                                                      |(10/10) Simulating: |=====                                                     |(10/10) Simulating: |======                                                    |(10/10) Simulating: |=======                                                   |(10/10) Simulating: |========                                                  |(10/10) Simulating: |=========                                                 |(10/10) Simulating: |==========                                                |(10/10) Simulating: |===========                                               |(10/10) Simulating: |============                                              |(10/10) Simulating: |=============                                             |(10/10) Simulating: |==============                                            |(10/10) Simulating: |===============                                           |(10/10) Simulating: |================                                          |(10/10) Simulating: |=================                                         |(10/10) Simulating: |==================                                        |(10/10) Simulating: |===================                                       |(10/10) Simulating: |====================                                      |(10/10) Simulating: |=====================                                     |(10/10) Simulating: |======================                                    |(10/10) Simulating: |=======================                                   |(10/10) Simulating: |========================                                  |(10/10) Simulating: |=========================                                 |(10/10) Simulating: |==========================                                |(10/10) Simulating: |===========================                               |(10/10) Simulating: |============================                              |(10/10) Simulating: |=============================                             |(10/10) Simulating: |==============================                            |(10/10) Simulating: |===============================                           |(10/10) Simulating: |================================                          |(10/10) Simulating: |=================================                         |(10/10) Simulating: |==================================                        |(10/10) Simulating: |===================================                       |(10/10) Simulating: |====================================                      |(10/10) Simulating: |=====================================                     |(10/10) Simulating: |======================================                    |(10/10) Simulating: |=======================================                   |(10/10) Simulating: |========================================                  |(10/10) Simulating: |=========================================                 |(10/10) Simulating: |==========================================                |(10/10) Simulating: |===========================================               |(10/10) Simulating: |============================================              |(10/10) Simulating: |=============================================             |(10/10) Simulating: |==============================================            |(10/10) Simulating: |===============================================           |(10/10) Simulating: |================================================          |(10/10) Simulating: |=================================================         |(10/10) Simulating: |==================================================        |(10/10) Simulating: |===================================================       |(10/10) Simulating: |====================================================      |(10/10) Simulating: |=====================================================     |(10/10) Simulating: |======================================================    |(10/10) Simulating: |=======================================================   |(10/10) Simulating: |========================================================  |(10/10) Simulating: |========================================================= |(10/10) Simulating: |==========================================================|(10/10) 
curve_dualism <- powerCurve(model_sim,
                    fixed("RDM_typephysical_psychological_model", "t"),
                    along = "subject_id",
                    breaks = c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
                    nsim = 100)
## Simulating: |                                                                  |Simulating: |=                                                                 |Simulating: |==                                                                |Simulating: |===                                                               |Simulating: |====                                                              |Simulating: |=====                                                             |Simulating: |======                                                            |Simulating: |=======                                                           |Simulating: |========                                                          |Simulating: |=========                                                         |Simulating: |==========                                                        |Simulating: |===========                                                       |Simulating: |============                                                      |Simulating: |=============                                                     |Simulating: |==============                                                    |Simulating: |===============                                                   |Simulating: |================                                                  |Simulating: |=================                                                 |Simulating: |==================                                                |Simulating: |===================                                               |Simulating: |====================                                              |Simulating: |=====================                                             |Simulating: |======================                                            |Simulating: |=======================                                           |Simulating: |========================                                          |Simulating: |=========================                                         |Simulating: |==========================                                        |Simulating: |===========================                                       |Simulating: |============================                                      |Simulating: |=============================                                     |Simulating: |==============================                                    |Simulating: |===============================                                   |Simulating: |================================                                  |Simulating: |=================================                                 |Simulating: |==================================                                |Simulating: |===================================                               |Simulating: |====================================                              |Simulating: |=====================================                             |Simulating: |======================================                            |Simulating: |=======================================                           |Simulating: |========================================                          |Simulating: |=========================================                         |Simulating: |==========================================                        |Simulating: |===========================================                       |Simulating: |============================================                      |Simulating: |=============================================                     |Simulating: |==============================================                    |Simulating: |===============================================                   |Simulating: |================================================                  |Simulating: |=================================================                 |Simulating: |==================================================                |Simulating: |===================================================               |Simulating: |====================================================              |Simulating: |=====================================================             |Simulating: |======================================================            |Simulating: |=======================================================           |Simulating: |========================================================          |Simulating: |=========================================================         |Simulating: |==========================================================        |Simulating: |===========================================================       |Simulating: |============================================================      |Simulating: |=============================================================     |Simulating: |==============================================================    |Simulating: |===============================================================   |Simulating: |================================================================  |Simulating: |================================================================= |Simulating: |==================================================================|( 1/10) ( 1/10) Simulating: |                                                          |( 1/10) Simulating: |=                                                         |( 1/10) Simulating: |==                                                        |( 1/10) Simulating: |===                                                       |( 1/10) Simulating: |====                                                      |( 1/10) Simulating: |=====                                                     |( 1/10) Simulating: |======                                                    |( 1/10) Simulating: |=======                                                   |( 1/10) Simulating: |========                                                  |( 1/10) Simulating: |=========                                                 |( 1/10) Simulating: |==========                                                |( 1/10) Simulating: |===========                                               |( 1/10) Simulating: |============                                              |( 1/10) Simulating: |=============                                             |( 1/10) Simulating: |==============                                            |( 1/10) Simulating: |===============                                           |( 1/10) Simulating: |================                                          |( 1/10) Simulating: |=================                                         |( 1/10) Simulating: |==================                                        |( 1/10) Simulating: |===================                                       |( 1/10) Simulating: |====================                                      |( 1/10) Simulating: |=====================                                     |( 1/10) Simulating: |======================                                    |( 1/10) Simulating: |=======================                                   |( 1/10) Simulating: |========================                                  |( 1/10) Simulating: |=========================                                 |( 1/10) Simulating: |==========================                                |( 1/10) Simulating: |===========================                               |( 1/10) Simulating: |============================                              |( 1/10) Simulating: |=============================                             |( 1/10) Simulating: |==============================                            |( 1/10) Simulating: |===============================                           |( 1/10) Simulating: |================================                          |( 1/10) Simulating: |=================================                         |( 1/10) Simulating: |==================================                        |( 1/10) Simulating: |===================================                       |( 1/10) Simulating: |====================================                      |( 1/10) Simulating: |=====================================                     |( 1/10) Simulating: |======================================                    |( 1/10) Simulating: |=======================================                   |( 1/10) Simulating: |========================================                  |( 1/10) Simulating: |=========================================                 |( 1/10) Simulating: |==========================================                |( 1/10) Simulating: |===========================================               |( 1/10) Simulating: |============================================              |( 1/10) Simulating: |=============================================             |( 1/10) Simulating: |==============================================            |( 1/10) Simulating: |===============================================           |( 1/10) Simulating: |================================================          |( 1/10) Simulating: |=================================================         |( 1/10) Simulating: |==================================================        |( 1/10) Simulating: |===================================================       |( 1/10) Simulating: |====================================================      |( 1/10) Simulating: |=====================================================     |( 1/10) Simulating: |======================================================    |( 1/10) Simulating: |=======================================================   |( 1/10) Simulating: |========================================================  |( 1/10) Simulating: |========================================================= |( 1/10) Simulating: |==========================================================|( 1/10) ( 2/10) ( 2/10) Simulating: |                                                          |( 2/10) Simulating: |=                                                         |( 2/10) Simulating: |==                                                        |( 2/10) Simulating: |===                                                       |( 2/10) Simulating: |====                                                      |( 2/10) Simulating: |=====                                                     |( 2/10) Simulating: |======                                                    |( 2/10) Simulating: |=======                                                   |( 2/10) Simulating: |========                                                  |( 2/10) Simulating: |=========                                                 |( 2/10) Simulating: |==========                                                |( 2/10) Simulating: |===========                                               |( 2/10) Simulating: |============                                              |( 2/10) Simulating: |=============                                             |( 2/10) Simulating: |==============                                            |( 2/10) Simulating: |===============                                           |( 2/10) Simulating: |================                                          |( 2/10) Simulating: |=================                                         |( 2/10) Simulating: |==================                                        |( 2/10) Simulating: |===================                                       |( 2/10) Simulating: |====================                                      |( 2/10) Simulating: |=====================                                     |( 2/10) Simulating: |======================                                    |( 2/10) Simulating: |=======================                                   |( 2/10) Simulating: |========================                                  |( 2/10) Simulating: |=========================                                 |( 2/10) Simulating: |==========================                                |( 2/10) Simulating: |===========================                               |( 2/10) Simulating: |============================                              |( 2/10) Simulating: |=============================                             |( 2/10) Simulating: |==============================                            |( 2/10) Simulating: |===============================                           |( 2/10) Simulating: |================================                          |( 2/10) Simulating: |=================================                         |( 2/10) Simulating: |==================================                        |( 2/10) Simulating: |===================================                       |( 2/10) Simulating: |====================================                      |( 2/10) Simulating: |=====================================                     |( 2/10) Simulating: |======================================                    |( 2/10) Simulating: |=======================================                   |( 2/10) Simulating: |========================================                  |( 2/10) Simulating: |=========================================                 |( 2/10) Simulating: |==========================================                |( 2/10) Simulating: |===========================================               |( 2/10) Simulating: |============================================              |( 2/10) Simulating: |=============================================             |( 2/10) Simulating: |==============================================            |( 2/10) Simulating: |===============================================           |( 2/10) Simulating: |================================================          |( 2/10) Simulating: |=================================================         |( 2/10) Simulating: |==================================================        |( 2/10) Simulating: |===================================================       |( 2/10) Simulating: |====================================================      |( 2/10) Simulating: |=====================================================     |( 2/10) Simulating: |======================================================    |( 2/10) Simulating: |=======================================================   |( 2/10) Simulating: |========================================================  |( 2/10) Simulating: |========================================================= |( 2/10) Simulating: |==========================================================|( 2/10) ( 3/10) ( 3/10) Simulating: |                                                          |( 3/10) Simulating: |=                                                         |( 3/10) Simulating: |==                                                        |( 3/10) Simulating: |===                                                       |( 3/10) Simulating: |====                                                      |( 3/10) Simulating: |=====                                                     |( 3/10) Simulating: |======                                                    |( 3/10) Simulating: |=======                                                   |( 3/10) Simulating: |========                                                  |( 3/10) Simulating: |=========                                                 |( 3/10) Simulating: |==========                                                |( 3/10) Simulating: |===========                                               |( 3/10) Simulating: |============                                              |( 3/10) Simulating: |=============                                             |( 3/10) Simulating: |==============                                            |( 3/10) Simulating: |===============                                           |( 3/10) Simulating: |================                                          |( 3/10) Simulating: |=================                                         |( 3/10) Simulating: |==================                                        |( 3/10) Simulating: |===================                                       |( 3/10) Simulating: |====================                                      |( 3/10) Simulating: |=====================                                     |( 3/10) Simulating: |======================                                    |( 3/10) Simulating: |=======================                                   |( 3/10) Simulating: |========================                                  |( 3/10) Simulating: |=========================                                 |( 3/10) Simulating: |==========================                                |( 3/10) Simulating: |===========================                               |( 3/10) Simulating: |============================                              |( 3/10) Simulating: |=============================                             |( 3/10) Simulating: |==============================                            |( 3/10) Simulating: |===============================                           |( 3/10) Simulating: |================================                          |( 3/10) Simulating: |=================================                         |( 3/10) Simulating: |==================================                        |( 3/10) Simulating: |===================================                       |( 3/10) Simulating: |====================================                      |( 3/10) Simulating: |=====================================                     |( 3/10) Simulating: |======================================                    |( 3/10) Simulating: |=======================================                   |( 3/10) Simulating: |========================================                  |( 3/10) Simulating: |=========================================                 |( 3/10) Simulating: |==========================================                |( 3/10) Simulating: |===========================================               |( 3/10) Simulating: |============================================              |( 3/10) Simulating: |=============================================             |( 3/10) Simulating: |==============================================            |( 3/10) Simulating: |===============================================           |( 3/10) Simulating: |================================================          |( 3/10) Simulating: |=================================================         |( 3/10) Simulating: |==================================================        |( 3/10) Simulating: |===================================================       |( 3/10) Simulating: |====================================================      |( 3/10) Simulating: |=====================================================     |( 3/10) Simulating: |======================================================    |( 3/10) Simulating: |=======================================================   |( 3/10) Simulating: |========================================================  |( 3/10) Simulating: |========================================================= |( 3/10) Simulating: |==========================================================|( 3/10) ( 4/10) ( 4/10) Simulating: |                                                          |( 4/10) Simulating: |=                                                         |( 4/10) Simulating: |==                                                        |( 4/10) Simulating: |===                                                       |( 4/10) Simulating: |====                                                      |( 4/10) Simulating: |=====                                                     |( 4/10) Simulating: |======                                                    |( 4/10) Simulating: |=======                                                   |( 4/10) Simulating: |========                                                  |( 4/10) Simulating: |=========                                                 |( 4/10) Simulating: |==========                                                |( 4/10) Simulating: |===========                                               |( 4/10) Simulating: |============                                              |( 4/10) Simulating: |=============                                             |( 4/10) Simulating: |==============                                            |( 4/10) Simulating: |===============                                           |( 4/10) Simulating: |================                                          |( 4/10) Simulating: |=================                                         |( 4/10) Simulating: |==================                                        |( 4/10) Simulating: |===================                                       |( 4/10) Simulating: |====================                                      |( 4/10) Simulating: |=====================                                     |( 4/10) Simulating: |======================                                    |( 4/10) Simulating: |=======================                                   |( 4/10) Simulating: |========================                                  |( 4/10) Simulating: |=========================                                 |( 4/10) Simulating: |==========================                                |( 4/10) Simulating: |===========================                               |( 4/10) Simulating: |============================                              |( 4/10) Simulating: |=============================                             |( 4/10) Simulating: |==============================                            |( 4/10) Simulating: |===============================                           |( 4/10) Simulating: |================================                          |( 4/10) Simulating: |=================================                         |( 4/10) Simulating: |==================================                        |( 4/10) Simulating: |===================================                       |( 4/10) Simulating: |====================================                      |( 4/10) Simulating: |=====================================                     |( 4/10) Simulating: |======================================                    |( 4/10) Simulating: |=======================================                   |( 4/10) Simulating: |========================================                  |( 4/10) Simulating: |=========================================                 |( 4/10) Simulating: |==========================================                |( 4/10) Simulating: |===========================================               |( 4/10) Simulating: |============================================              |( 4/10) Simulating: |=============================================             |( 4/10) Simulating: |==============================================            |( 4/10) Simulating: |===============================================           |( 4/10) Simulating: |================================================          |( 4/10) Simulating: |=================================================         |( 4/10) Simulating: |==================================================        |( 4/10) Simulating: |===================================================       |( 4/10) Simulating: |====================================================      |( 4/10) Simulating: |=====================================================     |( 4/10) Simulating: |======================================================    |( 4/10) Simulating: |=======================================================   |( 4/10) Simulating: |========================================================  |( 4/10) Simulating: |========================================================= |( 4/10) Simulating: |==========================================================|( 4/10) ( 5/10) ( 5/10) Simulating: |                                                          |( 5/10) Simulating: |=                                                         |( 5/10) Simulating: |==                                                        |( 5/10) Simulating: |===                                                       |( 5/10) Simulating: |====                                                      |( 5/10) Simulating: |=====                                                     |( 5/10) Simulating: |======                                                    |( 5/10) Simulating: |=======                                                   |( 5/10) Simulating: |========                                                  |( 5/10) Simulating: |=========                                                 |( 5/10) Simulating: |==========                                                |( 5/10) Simulating: |===========                                               |( 5/10) Simulating: |============                                              |( 5/10) Simulating: |=============                                             |( 5/10) Simulating: |==============                                            |( 5/10) Simulating: |===============                                           |( 5/10) Simulating: |================                                          |( 5/10) Simulating: |=================                                         |( 5/10) Simulating: |==================                                        |( 5/10) Simulating: |===================                                       |( 5/10) Simulating: |====================                                      |( 5/10) Simulating: |=====================                                     |( 5/10) Simulating: |======================                                    |( 5/10) Simulating: |=======================                                   |( 5/10) Simulating: |========================                                  |( 5/10) Simulating: |=========================                                 |( 5/10) Simulating: |==========================                                |( 5/10) Simulating: |===========================                               |( 5/10) Simulating: |============================                              |( 5/10) Simulating: |=============================                             |( 5/10) Simulating: |==============================                            |( 5/10) Simulating: |===============================                           |( 5/10) Simulating: |================================                          |( 5/10) Simulating: |=================================                         |( 5/10) Simulating: |==================================                        |( 5/10) Simulating: |===================================                       |( 5/10) Simulating: |====================================                      |( 5/10) Simulating: |=====================================                     |( 5/10) Simulating: |======================================                    |( 5/10) Simulating: |=======================================                   |( 5/10) Simulating: |========================================                  |( 5/10) Simulating: |=========================================                 |( 5/10) Simulating: |==========================================                |( 5/10) Simulating: |===========================================               |( 5/10) Simulating: |============================================              |( 5/10) Simulating: |=============================================             |( 5/10) Simulating: |==============================================            |( 5/10) Simulating: |===============================================           |( 5/10) Simulating: |================================================          |( 5/10) Simulating: |=================================================         |( 5/10) Simulating: |==================================================        |( 5/10) Simulating: |===================================================       |( 5/10) Simulating: |====================================================      |( 5/10) Simulating: |=====================================================     |( 5/10) Simulating: |======================================================    |( 5/10) Simulating: |=======================================================   |( 5/10) Simulating: |========================================================  |( 5/10) Simulating: |========================================================= |( 5/10) Simulating: |==========================================================|( 5/10) ( 6/10) ( 6/10) Simulating: |                                                          |( 6/10) Simulating: |=                                                         |( 6/10) Simulating: |==                                                        |( 6/10) Simulating: |===                                                       |( 6/10) Simulating: |====                                                      |( 6/10) Simulating: |=====                                                     |( 6/10) Simulating: |======                                                    |( 6/10) Simulating: |=======                                                   |( 6/10) Simulating: |========                                                  |( 6/10) Simulating: |=========                                                 |( 6/10) Simulating: |==========                                                |( 6/10) Simulating: |===========                                               |( 6/10) Simulating: |============                                              |( 6/10) Simulating: |=============                                             |( 6/10) Simulating: |==============                                            |( 6/10) Simulating: |===============                                           |( 6/10) Simulating: |================                                          |( 6/10) Simulating: |=================                                         |( 6/10) Simulating: |==================                                        |( 6/10) Simulating: |===================                                       |( 6/10) Simulating: |====================                                      |( 6/10) Simulating: |=====================                                     |( 6/10) Simulating: |======================                                    |( 6/10) Simulating: |=======================                                   |( 6/10) Simulating: |========================                                  |( 6/10) Simulating: |=========================                                 |( 6/10) Simulating: |==========================                                |( 6/10) Simulating: |===========================                               |( 6/10) Simulating: |============================                              |( 6/10) Simulating: |=============================                             |( 6/10) Simulating: |==============================                            |( 6/10) Simulating: |===============================                           |( 6/10) Simulating: |================================                          |( 6/10) Simulating: |=================================                         |( 6/10) Simulating: |==================================                        |( 6/10) Simulating: |===================================                       |( 6/10) Simulating: |====================================                      |( 6/10) Simulating: |=====================================                     |( 6/10) Simulating: |======================================                    |( 6/10) Simulating: |=======================================                   |( 6/10) Simulating: |========================================                  |( 6/10) Simulating: |=========================================                 |( 6/10) Simulating: |==========================================                |( 6/10) Simulating: |===========================================               |( 6/10) Simulating: |============================================              |( 6/10) Simulating: |=============================================             |( 6/10) Simulating: |==============================================            |( 6/10) Simulating: |===============================================           |( 6/10) Simulating: |================================================          |( 6/10) Simulating: |=================================================         |( 6/10) Simulating: |==================================================        |( 6/10) Simulating: |===================================================       |( 6/10) Simulating: |====================================================      |( 6/10) Simulating: |=====================================================     |( 6/10) Simulating: |======================================================    |( 6/10) Simulating: |=======================================================   |( 6/10) Simulating: |========================================================  |( 6/10) Simulating: |========================================================= |( 6/10) Simulating: |==========================================================|( 6/10) ( 7/10) ( 7/10) Simulating: |                                                          |( 7/10) Simulating: |=                                                         |( 7/10) Simulating: |==                                                        |( 7/10) Simulating: |===                                                       |( 7/10) Simulating: |====                                                      |( 7/10) Simulating: |=====                                                     |( 7/10) Simulating: |======                                                    |( 7/10) Simulating: |=======                                                   |( 7/10) Simulating: |========                                                  |( 7/10) Simulating: |=========                                                 |( 7/10) Simulating: |==========                                                |( 7/10) Simulating: |===========                                               |( 7/10) Simulating: |============                                              |( 7/10) Simulating: |=============                                             |( 7/10) Simulating: |==============                                            |( 7/10) Simulating: |===============                                           |( 7/10) Simulating: |================                                          |( 7/10) Simulating: |=================                                         |( 7/10) Simulating: |==================                                        |( 7/10) Simulating: |===================                                       |( 7/10) Simulating: |====================                                      |( 7/10) Simulating: |=====================                                     |( 7/10) Simulating: |======================                                    |( 7/10) Simulating: |=======================                                   |( 7/10) Simulating: |========================                                  |( 7/10) Simulating: |=========================                                 |( 7/10) Simulating: |==========================                                |( 7/10) Simulating: |===========================                               |( 7/10) Simulating: |============================                              |( 7/10) Simulating: |=============================                             |( 7/10) Simulating: |==============================                            |( 7/10) Simulating: |===============================                           |( 7/10) Simulating: |================================                          |( 7/10) Simulating: |=================================                         |( 7/10) Simulating: |==================================                        |( 7/10) Simulating: |===================================                       |( 7/10) Simulating: |====================================                      |( 7/10) Simulating: |=====================================                     |( 7/10) Simulating: |======================================                    |( 7/10) Simulating: |=======================================                   |( 7/10) Simulating: |========================================                  |( 7/10) Simulating: |=========================================                 |( 7/10) Simulating: |==========================================                |( 7/10) Simulating: |===========================================               |( 7/10) Simulating: |============================================              |( 7/10) Simulating: |=============================================             |( 7/10) Simulating: |==============================================            |( 7/10) Simulating: |===============================================           |( 7/10) Simulating: |================================================          |( 7/10) Simulating: |=================================================         |( 7/10) Simulating: |==================================================        |( 7/10) Simulating: |===================================================       |( 7/10) Simulating: |====================================================      |( 7/10) Simulating: |=====================================================     |( 7/10) Simulating: |======================================================    |( 7/10) Simulating: |=======================================================   |( 7/10) Simulating: |========================================================  |( 7/10) Simulating: |========================================================= |( 7/10) Simulating: |==========================================================|( 7/10) ( 8/10) ( 8/10) Simulating: |                                                          |( 8/10) Simulating: |=                                                         |( 8/10) Simulating: |==                                                        |( 8/10) Simulating: |===                                                       |( 8/10) Simulating: |====                                                      |( 8/10) Simulating: |=====                                                     |( 8/10) Simulating: |======                                                    |( 8/10) Simulating: |=======                                                   |( 8/10) Simulating: |========                                                  |( 8/10) Simulating: |=========                                                 |( 8/10) Simulating: |==========                                                |( 8/10) Simulating: |===========                                               |( 8/10) Simulating: |============                                              |( 8/10) Simulating: |=============                                             |( 8/10) Simulating: |==============                                            |( 8/10) Simulating: |===============                                           |( 8/10) Simulating: |================                                          |( 8/10) Simulating: |=================                                         |( 8/10) Simulating: |==================                                        |( 8/10) Simulating: |===================                                       |( 8/10) Simulating: |====================                                      |( 8/10) Simulating: |=====================                                     |( 8/10) Simulating: |======================                                    |( 8/10) Simulating: |=======================                                   |( 8/10) Simulating: |========================                                  |( 8/10) Simulating: |=========================                                 |( 8/10) Simulating: |==========================                                |( 8/10) Simulating: |===========================                               |( 8/10) Simulating: |============================                              |( 8/10) Simulating: |=============================                             |( 8/10) Simulating: |==============================                            |( 8/10) Simulating: |===============================                           |( 8/10) Simulating: |================================                          |( 8/10) Simulating: |=================================                         |( 8/10) Simulating: |==================================                        |( 8/10) Simulating: |===================================                       |( 8/10) Simulating: |====================================                      |( 8/10) Simulating: |=====================================                     |( 8/10) Simulating: |======================================                    |( 8/10) Simulating: |=======================================                   |( 8/10) Simulating: |========================================                  |( 8/10) Simulating: |=========================================                 |( 8/10) Simulating: |==========================================                |( 8/10) Simulating: |===========================================               |( 8/10) Simulating: |============================================              |( 8/10) Simulating: |=============================================             |( 8/10) Simulating: |==============================================            |( 8/10) Simulating: |===============================================           |( 8/10) Simulating: |================================================          |( 8/10) Simulating: |=================================================         |( 8/10) Simulating: |==================================================        |( 8/10) Simulating: |===================================================       |( 8/10) Simulating: |====================================================      |( 8/10) Simulating: |=====================================================     |( 8/10) Simulating: |======================================================    |( 8/10) Simulating: |=======================================================   |( 8/10) Simulating: |========================================================  |( 8/10) Simulating: |========================================================= |( 8/10) Simulating: |==========================================================|( 8/10) ( 9/10) ( 9/10) Simulating: |                                                          |( 9/10) Simulating: |=                                                         |( 9/10) Simulating: |==                                                        |( 9/10) Simulating: |===                                                       |( 9/10) Simulating: |====                                                      |( 9/10) Simulating: |=====                                                     |( 9/10) Simulating: |======                                                    |( 9/10) Simulating: |=======                                                   |( 9/10) Simulating: |========                                                  |( 9/10) Simulating: |=========                                                 |( 9/10) Simulating: |==========                                                |( 9/10) Simulating: |===========                                               |( 9/10) Simulating: |============                                              |( 9/10) Simulating: |=============                                             |( 9/10) Simulating: |==============                                            |( 9/10) Simulating: |===============                                           |( 9/10) Simulating: |================                                          |( 9/10) Simulating: |=================                                         |( 9/10) Simulating: |==================                                        |( 9/10) Simulating: |===================                                       |( 9/10) Simulating: |====================                                      |( 9/10) Simulating: |=====================                                     |( 9/10) Simulating: |======================                                    |( 9/10) Simulating: |=======================                                   |( 9/10) Simulating: |========================                                  |( 9/10) Simulating: |=========================                                 |( 9/10) Simulating: |==========================                                |( 9/10) Simulating: |===========================                               |( 9/10) Simulating: |============================                              |( 9/10) Simulating: |=============================                             |( 9/10) Simulating: |==============================                            |( 9/10) Simulating: |===============================                           |( 9/10) Simulating: |================================                          |( 9/10) Simulating: |=================================                         |( 9/10) Simulating: |==================================                        |( 9/10) Simulating: |===================================                       |( 9/10) Simulating: |====================================                      |( 9/10) Simulating: |=====================================                     |( 9/10) Simulating: |======================================                    |( 9/10) Simulating: |=======================================                   |( 9/10) Simulating: |========================================                  |( 9/10) Simulating: |=========================================                 |( 9/10) Simulating: |==========================================                |( 9/10) Simulating: |===========================================               |( 9/10) Simulating: |============================================              |( 9/10) Simulating: |=============================================             |( 9/10) Simulating: |==============================================            |( 9/10) Simulating: |===============================================           |( 9/10) Simulating: |================================================          |( 9/10) Simulating: |=================================================         |( 9/10) Simulating: |==================================================        |( 9/10) Simulating: |===================================================       |( 9/10) Simulating: |====================================================      |( 9/10) Simulating: |=====================================================     |( 9/10) Simulating: |======================================================    |( 9/10) Simulating: |=======================================================   |( 9/10) Simulating: |========================================================  |( 9/10) Simulating: |========================================================= |( 9/10) Simulating: |==========================================================|( 9/10) (10/10) (10/10) Simulating: |                                                          |(10/10) Simulating: |=                                                         |(10/10) Simulating: |==                                                        |(10/10) Simulating: |===                                                       |(10/10) Simulating: |====                                                      |(10/10) Simulating: |=====                                                     |(10/10) Simulating: |======                                                    |(10/10) Simulating: |=======                                                   |(10/10) Simulating: |========                                                  |(10/10) Simulating: |=========                                                 |(10/10) Simulating: |==========                                                |(10/10) Simulating: |===========                                               |(10/10) Simulating: |============                                              |(10/10) Simulating: |=============                                             |(10/10) Simulating: |==============                                            |(10/10) Simulating: |===============                                           |(10/10) Simulating: |================                                          |(10/10) Simulating: |=================                                         |(10/10) Simulating: |==================                                        |(10/10) Simulating: |===================                                       |(10/10) Simulating: |====================                                      |(10/10) Simulating: |=====================                                     |(10/10) Simulating: |======================                                    |(10/10) Simulating: |=======================                                   |(10/10) Simulating: |========================                                  |(10/10) Simulating: |=========================                                 |(10/10) Simulating: |==========================                                |(10/10) Simulating: |===========================                               |(10/10) Simulating: |============================                              |(10/10) Simulating: |=============================                             |(10/10) Simulating: |==============================                            |(10/10) Simulating: |===============================                           |(10/10) Simulating: |================================                          |(10/10) Simulating: |=================================                         |(10/10) Simulating: |==================================                        |(10/10) Simulating: |===================================                       |(10/10) Simulating: |====================================                      |(10/10) Simulating: |=====================================                     |(10/10) Simulating: |======================================                    |(10/10) Simulating: |=======================================                   |(10/10) Simulating: |========================================                  |(10/10) Simulating: |=========================================                 |(10/10) Simulating: |==========================================                |(10/10) Simulating: |===========================================               |(10/10) Simulating: |============================================              |(10/10) Simulating: |=============================================             |(10/10) Simulating: |==============================================            |(10/10) Simulating: |===============================================           |(10/10) Simulating: |================================================          |(10/10) Simulating: |=================================================         |(10/10) Simulating: |==================================================        |(10/10) Simulating: |===================================================       |(10/10) Simulating: |====================================================      |(10/10) Simulating: |=====================================================     |(10/10) Simulating: |======================================================    |(10/10) Simulating: |=======================================================   |(10/10) Simulating: |========================================================  |(10/10) Simulating: |========================================================= |(10/10) Simulating: |==========================================================|(10/10) 
curve_embeddings <- powerCurve(model_sim,
                    fixed("RDM_typeuse_phrase_embeddings_model", "t"),
                    along = "subject_id",
                    breaks = c(10, 20, 30, 40, 50, 60, 70, 80, 90, 100),
                    nsim = 100)
## Simulating: |                                                                  |Simulating: |=                                                                 |Simulating: |==                                                                |Simulating: |===                                                               |Simulating: |====                                                              |Simulating: |=====                                                             |Simulating: |======                                                            |Simulating: |=======                                                           |Simulating: |========                                                          |Simulating: |=========                                                         |Simulating: |==========                                                        |Simulating: |===========                                                       |Simulating: |============                                                      |Simulating: |=============                                                     |Simulating: |==============                                                    |Simulating: |===============                                                   |Simulating: |================                                                  |Simulating: |=================                                                 |Simulating: |==================                                                |Simulating: |===================                                               |Simulating: |====================                                              |Simulating: |=====================                                             |Simulating: |======================                                            |Simulating: |=======================                                           |Simulating: |========================                                          |Simulating: |=========================                                         |Simulating: |==========================                                        |Simulating: |===========================                                       |Simulating: |============================                                      |Simulating: |=============================                                     |Simulating: |==============================                                    |Simulating: |===============================                                   |Simulating: |================================                                  |Simulating: |=================================                                 |Simulating: |==================================                                |Simulating: |===================================                               |Simulating: |====================================                              |Simulating: |=====================================                             |Simulating: |======================================                            |Simulating: |=======================================                           |Simulating: |========================================                          |Simulating: |=========================================                         |Simulating: |==========================================                        |Simulating: |===========================================                       |Simulating: |============================================                      |Simulating: |=============================================                     |Simulating: |==============================================                    |Simulating: |===============================================                   |Simulating: |================================================                  |Simulating: |=================================================                 |Simulating: |==================================================                |Simulating: |===================================================               |Simulating: |====================================================              |Simulating: |=====================================================             |Simulating: |======================================================            |Simulating: |=======================================================           |Simulating: |========================================================          |Simulating: |=========================================================         |Simulating: |==========================================================        |Simulating: |===========================================================       |Simulating: |============================================================      |Simulating: |=============================================================     |Simulating: |==============================================================    |Simulating: |===============================================================   |Simulating: |================================================================  |Simulating: |================================================================= |Simulating: |==================================================================|( 1/10) ( 1/10) Simulating: |                                                          |( 1/10) Simulating: |=                                                         |( 1/10) Simulating: |==                                                        |( 1/10) Simulating: |===                                                       |( 1/10) Simulating: |====                                                      |( 1/10) Simulating: |=====                                                     |( 1/10) Simulating: |======                                                    |( 1/10) Simulating: |=======                                                   |( 1/10) Simulating: |========                                                  |( 1/10) Simulating: |=========                                                 |( 1/10) Simulating: |==========                                                |( 1/10) Simulating: |===========                                               |( 1/10) Simulating: |============                                              |( 1/10) Simulating: |=============                                             |( 1/10) Simulating: |==============                                            |( 1/10) Simulating: |===============                                           |( 1/10) Simulating: |================                                          |( 1/10) Simulating: |=================                                         |( 1/10) Simulating: |==================                                        |( 1/10) Simulating: |===================                                       |( 1/10) Simulating: |====================                                      |( 1/10) Simulating: |=====================                                     |( 1/10) Simulating: |======================                                    |( 1/10) Simulating: |=======================                                   |( 1/10) Simulating: |========================                                  |( 1/10) Simulating: |=========================                                 |( 1/10) Simulating: |==========================                                |( 1/10) Simulating: |===========================                               |( 1/10) Simulating: |============================                              |( 1/10) Simulating: |=============================                             |( 1/10) Simulating: |==============================                            |( 1/10) Simulating: |===============================                           |( 1/10) Simulating: |================================                          |( 1/10) Simulating: |=================================                         |( 1/10) Simulating: |==================================                        |( 1/10) Simulating: |===================================                       |( 1/10) Simulating: |====================================                      |( 1/10) Simulating: |=====================================                     |( 1/10) Simulating: |======================================                    |( 1/10) Simulating: |=======================================                   |( 1/10) Simulating: |========================================                  |( 1/10) Simulating: |=========================================                 |( 1/10) Simulating: |==========================================                |( 1/10) Simulating: |===========================================               |( 1/10) Simulating: |============================================              |( 1/10) Simulating: |=============================================             |( 1/10) Simulating: |==============================================            |( 1/10) Simulating: |===============================================           |( 1/10) Simulating: |================================================          |( 1/10) Simulating: |=================================================         |( 1/10) Simulating: |==================================================        |( 1/10) Simulating: |===================================================       |( 1/10) Simulating: |====================================================      |( 1/10) Simulating: |=====================================================     |( 1/10) Simulating: |======================================================    |( 1/10) Simulating: |=======================================================   |( 1/10) Simulating: |========================================================  |( 1/10) Simulating: |========================================================= |( 1/10) Simulating: |==========================================================|( 1/10) ( 2/10) ( 2/10) Simulating: |                                                          |( 2/10) Simulating: |=                                                         |( 2/10) Simulating: |==                                                        |( 2/10) Simulating: |===                                                       |( 2/10) Simulating: |====                                                      |( 2/10) Simulating: |=====                                                     |( 2/10) Simulating: |======                                                    |( 2/10) Simulating: |=======                                                   |( 2/10) Simulating: |========                                                  |( 2/10) Simulating: |=========                                                 |( 2/10) Simulating: |==========                                                |( 2/10) Simulating: |===========                                               |( 2/10) Simulating: |============                                              |( 2/10) Simulating: |=============                                             |( 2/10) Simulating: |==============                                            |( 2/10) Simulating: |===============                                           |( 2/10) Simulating: |================                                          |( 2/10) Simulating: |=================                                         |( 2/10) Simulating: |==================                                        |( 2/10) Simulating: |===================                                       |( 2/10) Simulating: |====================                                      |( 2/10) Simulating: |=====================                                     |( 2/10) Simulating: |======================                                    |( 2/10) Simulating: |=======================                                   |( 2/10) Simulating: |========================                                  |( 2/10) Simulating: |=========================                                 |( 2/10) Simulating: |==========================                                |( 2/10) Simulating: |===========================                               |( 2/10) Simulating: |============================                              |( 2/10) Simulating: |=============================                             |( 2/10) Simulating: |==============================                            |( 2/10) Simulating: |===============================                           |( 2/10) Simulating: |================================                          |( 2/10) Simulating: |=================================                         |( 2/10) Simulating: |==================================                        |( 2/10) Simulating: |===================================                       |( 2/10) Simulating: |====================================                      |( 2/10) Simulating: |=====================================                     |( 2/10) Simulating: |======================================                    |( 2/10) Simulating: |=======================================                   |( 2/10) Simulating: |========================================                  |( 2/10) Simulating: |=========================================                 |( 2/10) Simulating: |==========================================                |( 2/10) Simulating: |===========================================               |( 2/10) Simulating: |============================================              |( 2/10) Simulating: |=============================================             |( 2/10) Simulating: |==============================================            |( 2/10) Simulating: |===============================================           |( 2/10) Simulating: |================================================          |( 2/10) Simulating: |=================================================         |( 2/10) Simulating: |==================================================        |( 2/10) Simulating: |===================================================       |( 2/10) Simulating: |====================================================      |( 2/10) Simulating: |=====================================================     |( 2/10) Simulating: |======================================================    |( 2/10) Simulating: |=======================================================   |( 2/10) Simulating: |========================================================  |( 2/10) Simulating: |========================================================= |( 2/10) Simulating: |==========================================================|( 2/10) ( 3/10) ( 3/10) Simulating: |                                                          |( 3/10) Simulating: |=                                                         |( 3/10) Simulating: |==                                                        |( 3/10) Simulating: |===                                                       |( 3/10) Simulating: |====                                                      |( 3/10) Simulating: |=====                                                     |( 3/10) Simulating: |======                                                    |( 3/10) Simulating: |=======                                                   |( 3/10) Simulating: |========                                                  |( 3/10) Simulating: |=========                                                 |( 3/10) Simulating: |==========                                                |( 3/10) Simulating: |===========                                               |( 3/10) Simulating: |============                                              |( 3/10) Simulating: |=============                                             |( 3/10) Simulating: |==============                                            |( 3/10) Simulating: |===============                                           |( 3/10) Simulating: |================                                          |( 3/10) Simulating: |=================                                         |( 3/10) Simulating: |==================                                        |( 3/10) Simulating: |===================                                       |( 3/10) Simulating: |====================                                      |( 3/10) Simulating: |=====================                                     |( 3/10) Simulating: |======================                                    |( 3/10) Simulating: |=======================                                   |( 3/10) Simulating: |========================                                  |( 3/10) Simulating: |=========================                                 |( 3/10) Simulating: |==========================                                |( 3/10) Simulating: |===========================                               |( 3/10) Simulating: |============================                              |( 3/10) Simulating: |=============================                             |( 3/10) Simulating: |==============================                            |( 3/10) Simulating: |===============================                           |( 3/10) Simulating: |================================                          |( 3/10) Simulating: |=================================                         |( 3/10) Simulating: |==================================                        |( 3/10) Simulating: |===================================                       |( 3/10) Simulating: |====================================                      |( 3/10) Simulating: |=====================================                     |( 3/10) Simulating: |======================================                    |( 3/10) Simulating: |=======================================                   |( 3/10) Simulating: |========================================                  |( 3/10) Simulating: |=========================================                 |( 3/10) Simulating: |==========================================                |( 3/10) Simulating: |===========================================               |( 3/10) Simulating: |============================================              |( 3/10) Simulating: |=============================================             |( 3/10) Simulating: |==============================================            |( 3/10) Simulating: |===============================================           |( 3/10) Simulating: |================================================          |( 3/10) Simulating: |=================================================         |( 3/10) Simulating: |==================================================        |( 3/10) Simulating: |===================================================       |( 3/10) Simulating: |====================================================      |( 3/10) Simulating: |=====================================================     |( 3/10) Simulating: |======================================================    |( 3/10) Simulating: |=======================================================   |( 3/10) Simulating: |========================================================  |( 3/10) Simulating: |========================================================= |( 3/10) Simulating: |==========================================================|( 3/10) ( 4/10) ( 4/10) Simulating: |                                                          |( 4/10) Simulating: |=                                                         |( 4/10) Simulating: |==                                                        |( 4/10) Simulating: |===                                                       |( 4/10) Simulating: |====                                                      |( 4/10) Simulating: |=====                                                     |( 4/10) Simulating: |======                                                    |( 4/10) Simulating: |=======                                                   |( 4/10) Simulating: |========                                                  |( 4/10) Simulating: |=========                                                 |( 4/10) Simulating: |==========                                                |( 4/10) Simulating: |===========                                               |( 4/10) Simulating: |============                                              |( 4/10) Simulating: |=============                                             |( 4/10) Simulating: |==============                                            |( 4/10) Simulating: |===============                                           |( 4/10) Simulating: |================                                          |( 4/10) Simulating: |=================                                         |( 4/10) Simulating: |==================                                        |( 4/10) Simulating: |===================                                       |( 4/10) Simulating: |====================                                      |( 4/10) Simulating: |=====================                                     |( 4/10) Simulating: |======================                                    |( 4/10) Simulating: |=======================                                   |( 4/10) Simulating: |========================                                  |( 4/10) Simulating: |=========================                                 |( 4/10) Simulating: |==========================                                |( 4/10) Simulating: |===========================                               |( 4/10) Simulating: |============================                              |( 4/10) Simulating: |=============================                             |( 4/10) Simulating: |==============================                            |( 4/10) Simulating: |===============================                           |( 4/10) Simulating: |================================                          |( 4/10) Simulating: |=================================                         |( 4/10) Simulating: |==================================                        |( 4/10) Simulating: |===================================                       |( 4/10) Simulating: |====================================                      |( 4/10) Simulating: |=====================================                     |( 4/10) Simulating: |======================================                    |( 4/10) Simulating: |=======================================                   |( 4/10) Simulating: |========================================                  |( 4/10) Simulating: |=========================================                 |( 4/10) Simulating: |==========================================                |( 4/10) Simulating: |===========================================               |( 4/10) Simulating: |============================================              |( 4/10) Simulating: |=============================================             |( 4/10) Simulating: |==============================================            |( 4/10) Simulating: |===============================================           |( 4/10) Simulating: |================================================          |( 4/10) Simulating: |=================================================         |( 4/10) Simulating: |==================================================        |( 4/10) Simulating: |===================================================       |( 4/10) Simulating: |====================================================      |( 4/10) Simulating: |=====================================================     |( 4/10) Simulating: |======================================================    |( 4/10) Simulating: |=======================================================   |( 4/10) Simulating: |========================================================  |( 4/10) Simulating: |========================================================= |( 4/10) Simulating: |==========================================================|( 4/10) ( 5/10) ( 5/10) Simulating: |                                                          |( 5/10) Simulating: |=                                                         |( 5/10) Simulating: |==                                                        |( 5/10) Simulating: |===                                                       |( 5/10) Simulating: |====                                                      |( 5/10) Simulating: |=====                                                     |( 5/10) Simulating: |======                                                    |( 5/10) Simulating: |=======                                                   |( 5/10) Simulating: |========                                                  |( 5/10) Simulating: |=========                                                 |( 5/10) Simulating: |==========                                                |( 5/10) Simulating: |===========                                               |( 5/10) Simulating: |============                                              |( 5/10) Simulating: |=============                                             |( 5/10) Simulating: |==============                                            |( 5/10) Simulating: |===============                                           |( 5/10) Simulating: |================                                          |( 5/10) Simulating: |=================                                         |( 5/10) Simulating: |==================                                        |( 5/10) Simulating: |===================                                       |( 5/10) Simulating: |====================                                      |( 5/10) Simulating: |=====================                                     |( 5/10) Simulating: |======================                                    |( 5/10) Simulating: |=======================                                   |( 5/10) Simulating: |========================                                  |( 5/10) Simulating: |=========================                                 |( 5/10) Simulating: |==========================                                |( 5/10) Simulating: |===========================                               |( 5/10) Simulating: |============================                              |( 5/10) Simulating: |=============================                             |( 5/10) Simulating: |==============================                            |( 5/10) Simulating: |===============================                           |( 5/10) Simulating: |================================                          |( 5/10) Simulating: |=================================                         |( 5/10) Simulating: |==================================                        |( 5/10) Simulating: |===================================                       |( 5/10) Simulating: |====================================                      |( 5/10) Simulating: |=====================================                     |( 5/10) Simulating: |======================================                    |( 5/10) Simulating: |=======================================                   |( 5/10) Simulating: |========================================                  |( 5/10) Simulating: |=========================================                 |( 5/10) Simulating: |==========================================                |( 5/10) Simulating: |===========================================               |( 5/10) Simulating: |============================================              |( 5/10) Simulating: |=============================================             |( 5/10) Simulating: |==============================================            |( 5/10) Simulating: |===============================================           |( 5/10) Simulating: |================================================          |( 5/10) Simulating: |=================================================         |( 5/10) Simulating: |==================================================        |( 5/10) Simulating: |===================================================       |( 5/10) Simulating: |====================================================      |( 5/10) Simulating: |=====================================================     |( 5/10) Simulating: |======================================================    |( 5/10) Simulating: |=======================================================   |( 5/10) Simulating: |========================================================  |( 5/10) Simulating: |========================================================= |( 5/10) Simulating: |==========================================================|( 5/10) ( 6/10) ( 6/10) Simulating: |                                                          |( 6/10) Simulating: |=                                                         |( 6/10) Simulating: |==                                                        |( 6/10) Simulating: |===                                                       |( 6/10) Simulating: |====                                                      |( 6/10) Simulating: |=====                                                     |( 6/10) Simulating: |======                                                    |( 6/10) Simulating: |=======                                                   |( 6/10) Simulating: |========                                                  |( 6/10) Simulating: |=========                                                 |( 6/10) Simulating: |==========                                                |( 6/10) Simulating: |===========                                               |( 6/10) Simulating: |============                                              |( 6/10) Simulating: |=============                                             |( 6/10) Simulating: |==============                                            |( 6/10) Simulating: |===============                                           |( 6/10) Simulating: |================                                          |( 6/10) Simulating: |=================                                         |( 6/10) Simulating: |==================                                        |( 6/10) Simulating: |===================                                       |( 6/10) Simulating: |====================                                      |( 6/10) Simulating: |=====================                                     |( 6/10) Simulating: |======================                                    |( 6/10) Simulating: |=======================                                   |( 6/10) Simulating: |========================                                  |( 6/10) Simulating: |=========================                                 |( 6/10) Simulating: |==========================                                |( 6/10) Simulating: |===========================                               |( 6/10) Simulating: |============================                              |( 6/10) Simulating: |=============================                             |( 6/10) Simulating: |==============================                            |( 6/10) Simulating: |===============================                           |( 6/10) Simulating: |================================                          |( 6/10) Simulating: |=================================                         |( 6/10) Simulating: |==================================                        |( 6/10) Simulating: |===================================                       |( 6/10) Simulating: |====================================                      |( 6/10) Simulating: |=====================================                     |( 6/10) Simulating: |======================================                    |( 6/10) Simulating: |=======================================                   |( 6/10) Simulating: |========================================                  |( 6/10) Simulating: |=========================================                 |( 6/10) Simulating: |==========================================                |( 6/10) Simulating: |===========================================               |( 6/10) Simulating: |============================================              |( 6/10) Simulating: |=============================================             |( 6/10) Simulating: |==============================================            |( 6/10) Simulating: |===============================================           |( 6/10) Simulating: |================================================          |( 6/10) Simulating: |=================================================         |( 6/10) Simulating: |==================================================        |( 6/10) Simulating: |===================================================       |( 6/10) Simulating: |====================================================      |( 6/10) Simulating: |=====================================================     |( 6/10) Simulating: |======================================================    |( 6/10) Simulating: |=======================================================   |( 6/10) Simulating: |========================================================  |( 6/10) Simulating: |========================================================= |( 6/10) Simulating: |==========================================================|( 6/10) ( 7/10) ( 7/10) Simulating: |                                                          |( 7/10) Simulating: |=                                                         |( 7/10) Simulating: |==                                                        |( 7/10) Simulating: |===                                                       |( 7/10) Simulating: |====                                                      |( 7/10) Simulating: |=====                                                     |( 7/10) Simulating: |======                                                    |( 7/10) Simulating: |=======                                                   |( 7/10) Simulating: |========                                                  |( 7/10) Simulating: |=========                                                 |( 7/10) Simulating: |==========                                                |( 7/10) Simulating: |===========                                               |( 7/10) Simulating: |============                                              |( 7/10) Simulating: |=============                                             |( 7/10) Simulating: |==============                                            |( 7/10) Simulating: |===============                                           |( 7/10) Simulating: |================                                          |( 7/10) Simulating: |=================                                         |( 7/10) Simulating: |==================                                        |( 7/10) Simulating: |===================                                       |( 7/10) Simulating: |====================                                      |( 7/10) Simulating: |=====================                                     |( 7/10) Simulating: |======================                                    |( 7/10) Simulating: |=======================                                   |( 7/10) Simulating: |========================                                  |( 7/10) Simulating: |=========================                                 |( 7/10) Simulating: |==========================                                |( 7/10) Simulating: |===========================                               |( 7/10) Simulating: |============================                              |( 7/10) Simulating: |=============================                             |( 7/10) Simulating: |==============================                            |( 7/10) Simulating: |===============================                           |( 7/10) Simulating: |================================                          |( 7/10) Simulating: |=================================                         |( 7/10) Simulating: |==================================                        |( 7/10) Simulating: |===================================                       |( 7/10) Simulating: |====================================                      |( 7/10) Simulating: |=====================================                     |( 7/10) Simulating: |======================================                    |( 7/10) Simulating: |=======================================                   |( 7/10) Simulating: |========================================                  |( 7/10) Simulating: |=========================================                 |( 7/10) Simulating: |==========================================                |( 7/10) Simulating: |===========================================               |( 7/10) Simulating: |============================================              |( 7/10) Simulating: |=============================================             |( 7/10) Simulating: |==============================================            |( 7/10) Simulating: |===============================================           |( 7/10) Simulating: |================================================          |( 7/10) Simulating: |=================================================         |( 7/10) Simulating: |==================================================        |( 7/10) Simulating: |===================================================       |( 7/10) Simulating: |====================================================      |( 7/10) Simulating: |=====================================================     |( 7/10) Simulating: |======================================================    |( 7/10) Simulating: |=======================================================   |( 7/10) Simulating: |========================================================  |( 7/10) Simulating: |========================================================= |( 7/10) Simulating: |==========================================================|( 7/10) ( 8/10) ( 8/10) Simulating: |                                                          |( 8/10) Simulating: |=                                                         |( 8/10) Simulating: |==                                                        |( 8/10) Simulating: |===                                                       |( 8/10) Simulating: |====                                                      |( 8/10) Simulating: |=====                                                     |( 8/10) Simulating: |======                                                    |( 8/10) Simulating: |=======                                                   |( 8/10) Simulating: |========                                                  |( 8/10) Simulating: |=========                                                 |( 8/10) Simulating: |==========                                                |( 8/10) Simulating: |===========                                               |( 8/10) Simulating: |============                                              |( 8/10) Simulating: |=============                                             |( 8/10) Simulating: |==============                                            |( 8/10) Simulating: |===============                                           |( 8/10) Simulating: |================                                          |( 8/10) Simulating: |=================                                         |( 8/10) Simulating: |==================                                        |( 8/10) Simulating: |===================                                       |( 8/10) Simulating: |====================                                      |( 8/10) Simulating: |=====================                                     |( 8/10) Simulating: |======================                                    |( 8/10) Simulating: |=======================                                   |( 8/10) Simulating: |========================                                  |( 8/10) Simulating: |=========================                                 |( 8/10) Simulating: |==========================                                |( 8/10) Simulating: |===========================                               |( 8/10) Simulating: |============================                              |( 8/10) Simulating: |=============================                             |( 8/10) Simulating: |==============================                            |( 8/10) Simulating: |===============================                           |( 8/10) Simulating: |================================                          |( 8/10) Simulating: |=================================                         |( 8/10) Simulating: |==================================                        |( 8/10) Simulating: |===================================                       |( 8/10) Simulating: |====================================                      |( 8/10) Simulating: |=====================================                     |( 8/10) Simulating: |======================================                    |( 8/10) Simulating: |=======================================                   |( 8/10) Simulating: |========================================                  |( 8/10) Simulating: |=========================================                 |( 8/10) Simulating: |==========================================                |( 8/10) Simulating: |===========================================               |( 8/10) Simulating: |============================================              |( 8/10) Simulating: |=============================================             |( 8/10) Simulating: |==============================================            |( 8/10) Simulating: |===============================================           |( 8/10) Simulating: |================================================          |( 8/10) Simulating: |=================================================         |( 8/10) Simulating: |==================================================        |( 8/10) Simulating: |===================================================       |( 8/10) Simulating: |====================================================      |( 8/10) Simulating: |=====================================================     |( 8/10) Simulating: |======================================================    |( 8/10) Simulating: |=======================================================   |( 8/10) Simulating: |========================================================  |( 8/10) Simulating: |========================================================= |( 8/10) Simulating: |==========================================================|( 8/10) ( 9/10) ( 9/10) Simulating: |                                                          |( 9/10) Simulating: |=                                                         |( 9/10) Simulating: |==                                                        |( 9/10) Simulating: |===                                                       |( 9/10) Simulating: |====                                                      |( 9/10) Simulating: |=====                                                     |( 9/10) Simulating: |======                                                    |( 9/10) Simulating: |=======                                                   |( 9/10) Simulating: |========                                                  |( 9/10) Simulating: |=========                                                 |( 9/10) Simulating: |==========                                                |( 9/10) Simulating: |===========                                               |( 9/10) Simulating: |============                                              |( 9/10) Simulating: |=============                                             |( 9/10) Simulating: |==============                                            |( 9/10) Simulating: |===============                                           |( 9/10) Simulating: |================                                          |( 9/10) Simulating: |=================                                         |( 9/10) Simulating: |==================                                        |( 9/10) Simulating: |===================                                       |( 9/10) Simulating: |====================                                      |( 9/10) Simulating: |=====================                                     |( 9/10) Simulating: |======================                                    |( 9/10) Simulating: |=======================                                   |( 9/10) Simulating: |========================                                  |( 9/10) Simulating: |=========================                                 |( 9/10) Simulating: |==========================                                |( 9/10) Simulating: |===========================                               |( 9/10) Simulating: |============================                              |( 9/10) Simulating: |=============================                             |( 9/10) Simulating: |==============================                            |( 9/10) Simulating: |===============================                           |( 9/10) Simulating: |================================                          |( 9/10) Simulating: |=================================                         |( 9/10) Simulating: |==================================                        |( 9/10) Simulating: |===================================                       |( 9/10) Simulating: |====================================                      |( 9/10) Simulating: |=====================================                     |( 9/10) Simulating: |======================================                    |( 9/10) Simulating: |=======================================                   |( 9/10) Simulating: |========================================                  |( 9/10) Simulating: |=========================================                 |( 9/10) Simulating: |==========================================                |( 9/10) Simulating: |===========================================               |( 9/10) Simulating: |============================================              |( 9/10) Simulating: |=============================================             |( 9/10) Simulating: |==============================================            |( 9/10) Simulating: |===============================================           |( 9/10) Simulating: |================================================          |( 9/10) Simulating: |=================================================         |( 9/10) Simulating: |==================================================        |( 9/10) Simulating: |===================================================       |( 9/10) Simulating: |====================================================      |( 9/10) Simulating: |=====================================================     |( 9/10) Simulating: |======================================================    |( 9/10) Simulating: |=======================================================   |( 9/10) Simulating: |========================================================  |( 9/10) Simulating: |========================================================= |( 9/10) Simulating: |==========================================================|( 9/10) (10/10) (10/10) Simulating: |                                                          |(10/10) Simulating: |=                                                         |(10/10) Simulating: |==                                                        |(10/10) Simulating: |===                                                       |(10/10) Simulating: |====                                                      |(10/10) Simulating: |=====                                                     |(10/10) Simulating: |======                                                    |(10/10) Simulating: |=======                                                   |(10/10) Simulating: |========                                                  |(10/10) Simulating: |=========                                                 |(10/10) Simulating: |==========                                                |(10/10) Simulating: |===========                                               |(10/10) Simulating: |============                                              |(10/10) Simulating: |=============                                             |(10/10) Simulating: |==============                                            |(10/10) Simulating: |===============                                           |(10/10) Simulating: |================                                          |(10/10) Simulating: |=================                                         |(10/10) Simulating: |==================                                        |(10/10) Simulating: |===================                                       |(10/10) Simulating: |====================                                      |(10/10) Simulating: |=====================                                     |(10/10) Simulating: |======================                                    |(10/10) Simulating: |=======================                                   |(10/10) Simulating: |========================                                  |(10/10) Simulating: |=========================                                 |(10/10) Simulating: |==========================                                |(10/10) Simulating: |===========================                               |(10/10) Simulating: |============================                              |(10/10) Simulating: |=============================                             |(10/10) Simulating: |==============================                            |(10/10) Simulating: |===============================                           |(10/10) Simulating: |================================                          |(10/10) Simulating: |=================================                         |(10/10) Simulating: |==================================                        |(10/10) Simulating: |===================================                       |(10/10) Simulating: |====================================                      |(10/10) Simulating: |=====================================                     |(10/10) Simulating: |======================================                    |(10/10) Simulating: |=======================================                   |(10/10) Simulating: |========================================                  |(10/10) Simulating: |=========================================                 |(10/10) Simulating: |==========================================                |(10/10) Simulating: |===========================================               |(10/10) Simulating: |============================================              |(10/10) Simulating: |=============================================             |(10/10) Simulating: |==============================================            |(10/10) Simulating: |===============================================           |(10/10) Simulating: |================================================          |(10/10) Simulating: |=================================================         |(10/10) Simulating: |==================================================        |(10/10) Simulating: |===================================================       |(10/10) Simulating: |====================================================      |(10/10) Simulating: |=====================================================     |(10/10) Simulating: |======================================================    |(10/10) Simulating: |=======================================================   |(10/10) Simulating: |========================================================  |(10/10) Simulating: |========================================================= |(10/10) Simulating: |==========================================================|(10/10) 

8.2 plot power curves

# ## all curves
plot(curve_intercept)

print(curve_intercept)
## Power for predictor '(Intercept)', (95% confidence interval),
## by number of levels in subject_id:
##       2: 31.00% (22.13, 41.03) - 8 rows
##       5: 87.00% (78.80, 92.89) - 20 rows
##      10: 99.00% (94.55, 99.97) - 40 rows
##      20: 100.0% (96.38, 100.0) - 80 rows
##      30: 100.0% (96.38, 100.0) - 120 rows
##      40: 100.0% (96.38, 100.0) - 160 rows
##      50: 100.0% (96.38, 100.0) - 200 rows
##      60: 100.0% (96.38, 100.0) - 240 rows
##      70: 100.0% (96.38, 100.0) - 280 rows
##      80: 100.0% (96.38, 100.0) - 320 rows
##      90: 100.0% (96.38, 100.0) - 360 rows
##     100: 100.0% (96.38, 100.0) - 400 rows
## 
## Time elapsed: 0 h 0 m 30 s
plot(curve_fine_grained)

print(curve_fine_grained)
## Power for predictor 'RDM_typefinegrained_mind_action_physiology_model', (95% confidence interval),
## by number of levels in subject_id:
##      10: 30.00% (21.24, 39.98) - 40 rows
##      20: 49.00% (38.86, 59.20) - 80 rows
##      30: 63.00% (52.76, 72.44) - 120 rows
##      40: 74.00% (64.27, 82.26) - 160 rows
##      50: 82.00% (73.05, 88.97) - 200 rows
##      60: 90.00% (82.38, 95.10) - 240 rows
##      70: 92.00% (84.84, 96.48) - 280 rows
##      80: 99.00% (94.55, 99.97) - 320 rows
##      90: 100.0% (96.38, 100.0) - 360 rows
##     100: 99.00% (94.55, 99.97) - 400 rows
## 
## Time elapsed: 0 h 0 m 25 s
plot(curve_dualism)

print(curve_dualism)
## Power for predictor 'RDM_typephysical_psychological_model', (95% confidence interval),
## by number of levels in subject_id:
##      10: 71.00% (61.07, 79.64) - 40 rows
##      20: 92.00% (84.84, 96.48) - 80 rows
##      30: 98.00% (92.96, 99.76) - 120 rows
##      40: 100.0% (96.38, 100.0) - 160 rows
##      50: 100.0% (96.38, 100.0) - 200 rows
##      60: 100.0% (96.38, 100.0) - 240 rows
##      70: 100.0% (96.38, 100.0) - 280 rows
##      80: 100.0% (96.38, 100.0) - 320 rows
##      90: 100.0% (96.38, 100.0) - 360 rows
##     100: 100.0% (96.38, 100.0) - 400 rows
## 
## Time elapsed: 0 h 0 m 25 s
plot(curve_embeddings)

print(curve_embeddings)
## Power for predictor 'RDM_typeuse_phrase_embeddings_model', (95% confidence interval),
## by number of levels in subject_id:
##      10: 78.00% (68.61, 85.67) - 40 rows
##      20: 98.00% (92.96, 99.76) - 80 rows
##      30: 99.00% (94.55, 99.97) - 120 rows
##      40: 100.0% (96.38, 100.0) - 160 rows
##      50: 100.0% (96.38, 100.0) - 200 rows
##      60: 100.0% (96.38, 100.0) - 240 rows
##      70: 100.0% (96.38, 100.0) - 280 rows
##      80: 100.0% (96.38, 100.0) - 320 rows
##      90: 100.0% (96.38, 100.0) - 360 rows
##     100: 100.0% (96.38, 100.0) - 400 rows
## 
## Time elapsed: 0 h 0 m 25 s

8.3 simulate effect of MAP model on Causal v Freesort

# summary(model_prediction2)

# #suppose we want a study with 100 participants..
model_sim2 <- extend(model_prediction2, along = "subject_id", n = 100)
length(unique(simr::getData(model_sim2)$subject_id)) #confirm that N is now 100
## [1] 100
# #and we want to replicate the lower limit effect size, assuming pilot is inflated..
# #fixef(model_sim2)["(Intercept)"] <- -0.202

# #4
curve_difference <- powerCurve(model_sim2, fixed("taskcausal",
                                               "t"),
                    along = "subject_id", breaks = c(2, 5, 10, 20, 30, 40, 50,
                                                     60, 70, 80, 90, 100, 
                                                     110, 150), nsim = 100)
## Calculating power at 14 sample sizes along subject_id
## Simulating: |                                                                  |Simulating: |=                                                                 |Simulating: |==                                                                |Simulating: |===                                                               |Simulating: |====                                                              |Simulating: |=====                                                             |Simulating: |======                                                            |Simulating: |=======                                                           |Simulating: |========                                                          |Simulating: |=========                                                         |Simulating: |==========                                                        |Simulating: |===========                                                       |Simulating: |============                                                      |Simulating: |=============                                                     |Simulating: |==============                                                    |Simulating: |===============                                                   |Simulating: |================                                                  |Simulating: |=================                                                 |Simulating: |==================                                                |Simulating: |===================                                               |Simulating: |====================                                              |Simulating: |=====================                                             |Simulating: |======================                                            |Simulating: |=======================                                           |Simulating: |========================                                          |Simulating: |=========================                                         |Simulating: |==========================                                        |Simulating: |===========================                                       |Simulating: |============================                                      |Simulating: |=============================                                     |Simulating: |==============================                                    |Simulating: |===============================                                   |Simulating: |================================                                  |Simulating: |=================================                                 |Simulating: |==================================                                |Simulating: |===================================                               |Simulating: |====================================                              |Simulating: |=====================================                             |Simulating: |======================================                            |Simulating: |=======================================                           |Simulating: |========================================                          |Simulating: |=========================================                         |Simulating: |==========================================                        |Simulating: |===========================================                       |Simulating: |============================================                      |Simulating: |=============================================                     |Simulating: |==============================================                    |Simulating: |===============================================                   |Simulating: |================================================                  |Simulating: |=================================================                 |Simulating: |==================================================                |Simulating: |===================================================               |Simulating: |====================================================              |Simulating: |=====================================================             |Simulating: |======================================================            |Simulating: |=======================================================           |Simulating: |========================================================          |Simulating: |=========================================================         |Simulating: |==========================================================        |Simulating: |===========================================================       |Simulating: |============================================================      |Simulating: |=============================================================     |Simulating: |==============================================================    |Simulating: |===============================================================   |Simulating: |================================================================  |Simulating: |================================================================= |Simulating: |==================================================================|( 1/14) ( 1/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=                                                         |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |===                                                       |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |====                                                      |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=====                                                     |( 1/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=======                                                   |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |========                                                  |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=========                                                 |( 1/14) Simulating: |==========                                                |( 1/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=============                                             |( 1/14) Simulating: |==============                                            |( 1/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |================                                          |( 1/14) Simulating: |=================                                         |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==================                                        |( 1/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=====================                                     |( 1/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=======================                                   |( 1/14) Simulating: |========================                                  |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=========================                                 |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |===========================                               |( 1/14) Simulating: |============================                              |( 1/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==============================                            |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |================================                          |( 1/14) Simulating: |=================================                         |( 1/14) Simulating: |==================================                        |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |===================================                       |( 1/14) Simulating: |====================================                      |( 1/14) Simulating: |=====================================                     |( 1/14) Simulating: |======================================                    |( 1/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |========================================                  |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=========================================                 |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==========================================                |( 1/14) Simulating: |===========================================               |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |============================================              |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=============================================             |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==============================================            |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |================================================          |( 1/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==================================================        |( 1/14) Simulating: |===================================================       |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |====================================================      |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |=====================================================     |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |======================================================    |( 1/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |========================================================  |( 1/14) Simulating: |========================================================= |
## boundary (singular) fit: see help('isSingular')
## ( 1/14) Simulating: |==========================================================|( 1/14) ( 2/14) ( 2/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=                                                         |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===                                                       |( 2/14) Simulating: |====                                                      |( 2/14) Simulating: |=====                                                     |( 2/14) Simulating: |======                                                    |( 2/14) Simulating: |=======                                                   |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |========                                                  |( 2/14) Simulating: |=========                                                 |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==========                                                |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=============                                             |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==============                                            |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |================                                          |( 2/14) Simulating: |=================                                         |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==================                                        |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=====================                                     |( 2/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=======================                                   |( 2/14) Simulating: |========================                                  |( 2/14) Simulating: |=========================                                 |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===========================                               |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |============================                              |( 2/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==============================                            |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===============================                           |( 2/14) Simulating: |================================                          |( 2/14) Simulating: |=================================                         |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==================================                        |( 2/14) Simulating: |===================================                       |( 2/14) Simulating: |====================================                      |( 2/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |======================================                    |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |========================================                  |( 2/14) Simulating: |=========================================                 |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==========================================                |( 2/14) Simulating: |===========================================               |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |============================================              |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=============================================             |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==============================================            |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |================================================          |( 2/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==================================================        |( 2/14) Simulating: |===================================================       |( 2/14) Simulating: |====================================================      |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=====================================================     |( 2/14) Simulating: |======================================================    |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |========================================================  |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |========================================================= |
## boundary (singular) fit: see help('isSingular')
## ( 2/14) Simulating: |==========================================================|( 2/14) ( 3/14) ( 3/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=                                                         |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==                                                        |( 3/14) Simulating: |===                                                       |( 3/14) Simulating: |====                                                      |( 3/14) Simulating: |=====                                                     |( 3/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=======                                                   |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |========                                                  |( 3/14) Simulating: |=========                                                 |( 3/14) Simulating: |==========                                                |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=============                                             |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==============                                            |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |================                                          |( 3/14) Simulating: |=================                                         |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==================                                        |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=====================                                     |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=======================                                   |( 3/14) Simulating: |========================                                  |( 3/14) Simulating: |=========================                                 |( 3/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |===========================                               |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |============================                              |( 3/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==============================                            |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |================================                          |( 3/14) Simulating: |=================================                         |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==================================                        |( 3/14) Simulating: |===================================                       |( 3/14) Simulating: |====================================                      |( 3/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |======================================                    |( 3/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |========================================                  |( 3/14) Simulating: |=========================================                 |( 3/14) Simulating: |==========================================                |( 3/14) Simulating: |===========================================               |( 3/14) Simulating: |============================================              |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |=============================================             |( 3/14) Simulating: |==============================================            |( 3/14) Simulating: |===============================================           |( 3/14) Simulating: |================================================          |( 3/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |==================================================        |( 3/14) Simulating: |===================================================       |( 3/14) Simulating: |====================================================      |( 3/14) Simulating: |=====================================================     |
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |======================================================    |( 3/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 3/14) Simulating: |========================================================  |( 3/14) Simulating: |========================================================= |( 3/14) Simulating: |==========================================================|( 3/14) ( 4/14) ( 4/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=                                                         |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |==                                                        |( 4/14) Simulating: |===                                                       |( 4/14) Simulating: |====                                                      |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=====                                                     |( 4/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=======                                                   |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |========                                                  |( 4/14) Simulating: |=========                                                 |( 4/14) Simulating: |==========                                                |( 4/14) Simulating: |===========                                               |( 4/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=============                                             |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |==============                                            |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |================                                          |( 4/14) Simulating: |=================                                         |( 4/14) Simulating: |==================                                        |( 4/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |====================                                      |( 4/14) Simulating: |=====================                                     |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |========================                                  |( 4/14) Simulating: |=========================                                 |( 4/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |===========================                               |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |============================                              |( 4/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |==============================                            |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |===============================                           |( 4/14) Simulating: |================================                          |( 4/14) Simulating: |=================================                         |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |==================================                        |( 4/14) Simulating: |===================================                       |( 4/14) Simulating: |====================================                      |( 4/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |======================================                    |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |========================================                  |( 4/14) Simulating: |=========================================                 |( 4/14) Simulating: |==========================================                |( 4/14) Simulating: |===========================================               |( 4/14) Simulating: |============================================              |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |=============================================             |( 4/14) Simulating: |==============================================            |( 4/14) Simulating: |===============================================           |( 4/14) Simulating: |================================================          |( 4/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |==================================================        |( 4/14) Simulating: |===================================================       |( 4/14) Simulating: |====================================================      |( 4/14) Simulating: |=====================================================     |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |======================================================    |( 4/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## ( 4/14) Simulating: |========================================================  |( 4/14) Simulating: |========================================================= |( 4/14) Simulating: |==========================================================|( 4/14) ( 5/14) ( 5/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=                                                         |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==                                                        |( 5/14) Simulating: |===                                                       |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |====                                                      |( 5/14) Simulating: |=====                                                     |( 5/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=======                                                   |( 5/14) Simulating: |========                                                  |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=========                                                 |( 5/14) Simulating: |==========                                                |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |===========                                               |( 5/14) Simulating: |============                                              |( 5/14) Simulating: |=============                                             |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==============                                            |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |================                                          |( 5/14) Simulating: |=================                                         |( 5/14) Simulating: |==================                                        |( 5/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=====================                                     |( 5/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |========================                                  |( 5/14) Simulating: |=========================                                 |( 5/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |===========================                               |( 5/14) Simulating: |============================                              |( 5/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==============================                            |( 5/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |================================                          |( 5/14) Simulating: |=================================                         |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==================================                        |( 5/14) Simulating: |===================================                       |( 5/14) Simulating: |====================================                      |( 5/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |======================================                    |( 5/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |========================================                  |( 5/14) Simulating: |=========================================                 |( 5/14) Simulating: |==========================================                |( 5/14) Simulating: |===========================================               |( 5/14) Simulating: |============================================              |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |=============================================             |( 5/14) Simulating: |==============================================            |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |================================================          |( 5/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==================================================        |( 5/14) Simulating: |===================================================       |( 5/14) Simulating: |====================================================      |( 5/14) Simulating: |=====================================================     |( 5/14) Simulating: |======================================================    |( 5/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |========================================================  |( 5/14) Simulating: |========================================================= |
## boundary (singular) fit: see help('isSingular')
## ( 5/14) Simulating: |==========================================================|( 5/14) ( 6/14) ( 6/14) Simulating: |                                                          |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |=                                                         |( 6/14) Simulating: |==                                                        |( 6/14) Simulating: |===                                                       |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |====                                                      |( 6/14) Simulating: |=====                                                     |( 6/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |=======                                                   |( 6/14) Simulating: |========                                                  |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |=========                                                 |( 6/14) Simulating: |==========                                                |( 6/14) Simulating: |===========                                               |( 6/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |=============                                             |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |==============                                            |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |================                                          |( 6/14) Simulating: |=================                                         |( 6/14) Simulating: |==================                                        |( 6/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |=====================                                     |( 6/14) Simulating: |======================                                    |( 6/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |========================                                  |( 6/14) Simulating: |=========================                                 |( 6/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |===========================                               |( 6/14) Simulating: |============================                              |( 6/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |==============================                            |( 6/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |================================                          |( 6/14) Simulating: |=================================                         |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |==================================                        |( 6/14) Simulating: |===================================                       |( 6/14) Simulating: |====================================                      |( 6/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |======================================                    |( 6/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |========================================                  |( 6/14) Simulating: |=========================================                 |( 6/14) Simulating: |==========================================                |( 6/14) Simulating: |===========================================               |( 6/14) Simulating: |============================================              |( 6/14) Simulating: |=============================================             |( 6/14) Simulating: |==============================================            |( 6/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |================================================          |( 6/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |==================================================        |( 6/14) Simulating: |===================================================       |( 6/14) Simulating: |====================================================      |( 6/14) Simulating: |=====================================================     |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |======================================================    |( 6/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |========================================================  |( 6/14) Simulating: |========================================================= |
## boundary (singular) fit: see help('isSingular')
## ( 6/14) Simulating: |==========================================================|( 6/14) ( 7/14) ( 7/14) Simulating: |                                                          |( 7/14) Simulating: |=                                                         |( 7/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |===                                                       |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |====                                                      |( 7/14) Simulating: |=====                                                     |( 7/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |=======                                                   |( 7/14) Simulating: |========                                                  |( 7/14) Simulating: |=========                                                 |( 7/14) Simulating: |==========                                                |( 7/14) Simulating: |===========                                               |( 7/14) Simulating: |============                                              |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |=============                                             |( 7/14) Simulating: |==============                                            |( 7/14) Simulating: |===============                                           |( 7/14) Simulating: |================                                          |( 7/14) Simulating: |=================                                         |( 7/14) Simulating: |==================                                        |( 7/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |=====================                                     |( 7/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |========================                                  |( 7/14) Simulating: |=========================                                 |( 7/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |===========================                               |( 7/14) Simulating: |============================                              |( 7/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |==============================                            |( 7/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |================================                          |( 7/14) Simulating: |=================================                         |( 7/14) Simulating: |==================================                        |( 7/14) Simulating: |===================================                       |( 7/14) Simulating: |====================================                      |( 7/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |======================================                    |( 7/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |========================================                  |( 7/14) Simulating: |=========================================                 |( 7/14) Simulating: |==========================================                |( 7/14) Simulating: |===========================================               |( 7/14) Simulating: |============================================              |( 7/14) Simulating: |=============================================             |( 7/14) Simulating: |==============================================            |( 7/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |================================================          |( 7/14) Simulating: |=================================================         |( 7/14) Simulating: |==================================================        |( 7/14) Simulating: |===================================================       |( 7/14) Simulating: |====================================================      |( 7/14) Simulating: |=====================================================     |( 7/14) Simulating: |======================================================    |( 7/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## ( 7/14) Simulating: |========================================================  |( 7/14) Simulating: |========================================================= |( 7/14) Simulating: |==========================================================|( 7/14) ( 8/14) ( 8/14) Simulating: |                                                          |( 8/14) Simulating: |=                                                         |( 8/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |===                                                       |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |====                                                      |( 8/14) Simulating: |=====                                                     |( 8/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |=======                                                   |( 8/14) Simulating: |========                                                  |( 8/14) Simulating: |=========                                                 |( 8/14) Simulating: |==========                                                |( 8/14) Simulating: |===========                                               |( 8/14) Simulating: |============                                              |( 8/14) Simulating: |=============                                             |( 8/14) Simulating: |==============                                            |( 8/14) Simulating: |===============                                           |( 8/14) Simulating: |================                                          |( 8/14) Simulating: |=================                                         |( 8/14) Simulating: |==================                                        |( 8/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |=====================                                     |( 8/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |========================                                  |( 8/14) Simulating: |=========================                                 |( 8/14) Simulating: |==========================                                |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |===========================                               |( 8/14) Simulating: |============================                              |( 8/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |==============================                            |( 8/14) Simulating: |===============================                           |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |================================                          |( 8/14) Simulating: |=================================                         |( 8/14) Simulating: |==================================                        |( 8/14) Simulating: |===================================                       |( 8/14) Simulating: |====================================                      |( 8/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |======================================                    |( 8/14) Simulating: |=======================================                   |( 8/14) Simulating: |========================================                  |( 8/14) Simulating: |=========================================                 |( 8/14) Simulating: |==========================================                |( 8/14) Simulating: |===========================================               |( 8/14) Simulating: |============================================              |( 8/14) Simulating: |=============================================             |( 8/14) Simulating: |==============================================            |( 8/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |================================================          |( 8/14) Simulating: |=================================================         |
## boundary (singular) fit: see help('isSingular')
## ( 8/14) Simulating: |==================================================        |( 8/14) Simulating: |===================================================       |( 8/14) Simulating: |====================================================      |( 8/14) Simulating: |=====================================================     |( 8/14) Simulating: |======================================================    |( 8/14) Simulating: |=======================================================   |( 8/14) Simulating: |========================================================  |( 8/14) Simulating: |========================================================= |( 8/14) Simulating: |==========================================================|( 8/14) ( 9/14) ( 9/14) Simulating: |                                                          |( 9/14) Simulating: |=                                                         |( 9/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |===                                                       |( 9/14) Simulating: |====                                                      |( 9/14) Simulating: |=====                                                     |( 9/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |=======                                                   |( 9/14) Simulating: |========                                                  |( 9/14) Simulating: |=========                                                 |( 9/14) Simulating: |==========                                                |( 9/14) Simulating: |===========                                               |( 9/14) Simulating: |============                                              |( 9/14) Simulating: |=============                                             |( 9/14) Simulating: |==============                                            |( 9/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |================                                          |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |=================                                         |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |==================                                        |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |====================                                      |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |=====================                                     |( 9/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |========================                                  |( 9/14) Simulating: |=========================                                 |( 9/14) Simulating: |==========================                                |( 9/14) Simulating: |===========================                               |( 9/14) Simulating: |============================                              |( 9/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |==============================                            |( 9/14) Simulating: |===============================                           |( 9/14) Simulating: |================================                          |( 9/14) Simulating: |=================================                         |( 9/14) Simulating: |==================================                        |( 9/14) Simulating: |===================================                       |( 9/14) Simulating: |====================================                      |( 9/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |======================================                    |( 9/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |========================================                  |( 9/14) Simulating: |=========================================                 |( 9/14) Simulating: |==========================================                |( 9/14) Simulating: |===========================================               |( 9/14) Simulating: |============================================              |( 9/14) Simulating: |=============================================             |( 9/14) Simulating: |==============================================            |( 9/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |================================================          |( 9/14) Simulating: |=================================================         |( 9/14) Simulating: |==================================================        |( 9/14) Simulating: |===================================================       |( 9/14) Simulating: |====================================================      |( 9/14) Simulating: |=====================================================     |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |======================================================    |( 9/14) Simulating: |=======================================================   |( 9/14) Simulating: |========================================================  |
## boundary (singular) fit: see help('isSingular')
## ( 9/14) Simulating: |========================================================= |( 9/14) Simulating: |==========================================================|( 9/14) (10/14) (10/14) Simulating: |                                                          |(10/14) Simulating: |=                                                         |(10/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |===                                                       |(10/14) Simulating: |====                                                      |(10/14) Simulating: |=====                                                     |(10/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |=======                                                   |(10/14) Simulating: |========                                                  |(10/14) Simulating: |=========                                                 |(10/14) Simulating: |==========                                                |(10/14) Simulating: |===========                                               |(10/14) Simulating: |============                                              |(10/14) Simulating: |=============                                             |(10/14) Simulating: |==============                                            |(10/14) Simulating: |===============                                           |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |================                                          |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |=================                                         |(10/14) Simulating: |==================                                        |(10/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |====================                                      |(10/14) Simulating: |=====================                                     |(10/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |========================                                  |(10/14) Simulating: |=========================                                 |(10/14) Simulating: |==========================                                |(10/14) Simulating: |===========================                               |(10/14) Simulating: |============================                              |(10/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |==============================                            |(10/14) Simulating: |===============================                           |(10/14) Simulating: |================================                          |(10/14) Simulating: |=================================                         |(10/14) Simulating: |==================================                        |(10/14) Simulating: |===================================                       |(10/14) Simulating: |====================================                      |(10/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |======================================                    |(10/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## (10/14) Simulating: |========================================                  |(10/14) Simulating: |=========================================                 |(10/14) Simulating: |==========================================                |(10/14) Simulating: |===========================================               |(10/14) Simulating: |============================================              |(10/14) Simulating: |=============================================             |(10/14) Simulating: |==============================================            |(10/14) Simulating: |===============================================           |(10/14) Simulating: |================================================          |(10/14) Simulating: |=================================================         |(10/14) Simulating: |==================================================        |(10/14) Simulating: |===================================================       |(10/14) Simulating: |====================================================      |(10/14) Simulating: |=====================================================     |(10/14) Simulating: |======================================================    |(10/14) Simulating: |=======================================================   |(10/14) Simulating: |========================================================  |(10/14) Simulating: |========================================================= |(10/14) Simulating: |==========================================================|(10/14) (11/14) (11/14) Simulating: |                                                          |(11/14) Simulating: |=                                                         |(11/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |===                                                       |(11/14) Simulating: |====                                                      |(11/14) Simulating: |=====                                                     |(11/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |=======                                                   |(11/14) Simulating: |========                                                  |(11/14) Simulating: |=========                                                 |(11/14) Simulating: |==========                                                |(11/14) Simulating: |===========                                               |(11/14) Simulating: |============                                              |(11/14) Simulating: |=============                                             |(11/14) Simulating: |==============                                            |(11/14) Simulating: |===============                                           |(11/14) Simulating: |================                                          |(11/14) Simulating: |=================                                         |(11/14) Simulating: |==================                                        |(11/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |====================                                      |(11/14) Simulating: |=====================                                     |(11/14) Simulating: |======================                                    |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |========================                                  |(11/14) Simulating: |=========================                                 |(11/14) Simulating: |==========================                                |(11/14) Simulating: |===========================                               |(11/14) Simulating: |============================                              |(11/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |==============================                            |(11/14) Simulating: |===============================                           |(11/14) Simulating: |================================                          |(11/14) Simulating: |=================================                         |(11/14) Simulating: |==================================                        |(11/14) Simulating: |===================================                       |(11/14) Simulating: |====================================                      |(11/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |======================================                    |(11/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |========================================                  |(11/14) Simulating: |=========================================                 |(11/14) Simulating: |==========================================                |(11/14) Simulating: |===========================================               |(11/14) Simulating: |============================================              |(11/14) Simulating: |=============================================             |(11/14) Simulating: |==============================================            |(11/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |================================================          |(11/14) Simulating: |=================================================         |(11/14) Simulating: |==================================================        |(11/14) Simulating: |===================================================       |(11/14) Simulating: |====================================================      |(11/14) Simulating: |=====================================================     |(11/14) Simulating: |======================================================    |(11/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## (11/14) Simulating: |========================================================  |(11/14) Simulating: |========================================================= |(11/14) Simulating: |==========================================================|(11/14) (12/14) (12/14) Simulating: |                                                          |(12/14) Simulating: |=                                                         |(12/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |===                                                       |(12/14) Simulating: |====                                                      |(12/14) Simulating: |=====                                                     |(12/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |=======                                                   |(12/14) Simulating: |========                                                  |(12/14) Simulating: |=========                                                 |(12/14) Simulating: |==========                                                |(12/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |============                                              |(12/14) Simulating: |=============                                             |(12/14) Simulating: |==============                                            |(12/14) Simulating: |===============                                           |(12/14) Simulating: |================                                          |(12/14) Simulating: |=================                                         |(12/14) Simulating: |==================                                        |(12/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |====================                                      |(12/14) Simulating: |=====================                                     |(12/14) Simulating: |======================                                    |(12/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |========================                                  |(12/14) Simulating: |=========================                                 |(12/14) Simulating: |==========================                                |(12/14) Simulating: |===========================                               |(12/14) Simulating: |============================                              |(12/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |==============================                            |(12/14) Simulating: |===============================                           |(12/14) Simulating: |================================                          |(12/14) Simulating: |=================================                         |(12/14) Simulating: |==================================                        |(12/14) Simulating: |===================================                       |(12/14) Simulating: |====================================                      |(12/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |======================================                    |(12/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |========================================                  |(12/14) Simulating: |=========================================                 |(12/14) Simulating: |==========================================                |(12/14) Simulating: |===========================================               |(12/14) Simulating: |============================================              |(12/14) Simulating: |=============================================             |(12/14) Simulating: |==============================================            |(12/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |================================================          |(12/14) Simulating: |=================================================         |(12/14) Simulating: |==================================================        |(12/14) Simulating: |===================================================       |(12/14) Simulating: |====================================================      |(12/14) Simulating: |=====================================================     |(12/14) Simulating: |======================================================    |(12/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## (12/14) Simulating: |========================================================  |(12/14) Simulating: |========================================================= |(12/14) Simulating: |==========================================================|(12/14) (13/14) (13/14) Simulating: |                                                          |(13/14) Simulating: |=                                                         |(13/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |===                                                       |(13/14) Simulating: |====                                                      |(13/14) Simulating: |=====                                                     |(13/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |=======                                                   |(13/14) Simulating: |========                                                  |(13/14) Simulating: |=========                                                 |(13/14) Simulating: |==========                                                |(13/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |============                                              |(13/14) Simulating: |=============                                             |(13/14) Simulating: |==============                                            |(13/14) Simulating: |===============                                           |(13/14) Simulating: |================                                          |(13/14) Simulating: |=================                                         |(13/14) Simulating: |==================                                        |(13/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |====================                                      |(13/14) Simulating: |=====================                                     |(13/14) Simulating: |======================                                    |(13/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |========================                                  |(13/14) Simulating: |=========================                                 |(13/14) Simulating: |==========================                                |(13/14) Simulating: |===========================                               |(13/14) Simulating: |============================                              |(13/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |==============================                            |(13/14) Simulating: |===============================                           |(13/14) Simulating: |================================                          |(13/14) Simulating: |=================================                         |(13/14) Simulating: |==================================                        |(13/14) Simulating: |===================================                       |(13/14) Simulating: |====================================                      |(13/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |======================================                    |(13/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |========================================                  |(13/14) Simulating: |=========================================                 |(13/14) Simulating: |==========================================                |(13/14) Simulating: |===========================================               |(13/14) Simulating: |============================================              |(13/14) Simulating: |=============================================             |(13/14) Simulating: |==============================================            |(13/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |================================================          |(13/14) Simulating: |=================================================         |(13/14) Simulating: |==================================================        |(13/14) Simulating: |===================================================       |(13/14) Simulating: |====================================================      |(13/14) Simulating: |=====================================================     |(13/14) Simulating: |======================================================    |(13/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## (13/14) Simulating: |========================================================  |(13/14) Simulating: |========================================================= |(13/14) Simulating: |==========================================================|(13/14) (14/14) (14/14) Simulating: |                                                          |(14/14) Simulating: |=                                                         |(14/14) Simulating: |==                                                        |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |===                                                       |(14/14) Simulating: |====                                                      |(14/14) Simulating: |=====                                                     |(14/14) Simulating: |======                                                    |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |=======                                                   |(14/14) Simulating: |========                                                  |(14/14) Simulating: |=========                                                 |(14/14) Simulating: |==========                                                |(14/14) Simulating: |===========                                               |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |============                                              |(14/14) Simulating: |=============                                             |(14/14) Simulating: |==============                                            |(14/14) Simulating: |===============                                           |(14/14) Simulating: |================                                          |(14/14) Simulating: |=================                                         |(14/14) Simulating: |==================                                        |(14/14) Simulating: |===================                                       |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |====================                                      |(14/14) Simulating: |=====================                                     |(14/14) Simulating: |======================                                    |(14/14) Simulating: |=======================                                   |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |========================                                  |(14/14) Simulating: |=========================                                 |(14/14) Simulating: |==========================                                |(14/14) Simulating: |===========================                               |(14/14) Simulating: |============================                              |(14/14) Simulating: |=============================                             |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |==============================                            |(14/14) Simulating: |===============================                           |(14/14) Simulating: |================================                          |(14/14) Simulating: |=================================                         |(14/14) Simulating: |==================================                        |(14/14) Simulating: |===================================                       |(14/14) Simulating: |====================================                      |(14/14) Simulating: |=====================================                     |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |======================================                    |(14/14) Simulating: |=======================================                   |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |========================================                  |(14/14) Simulating: |=========================================                 |(14/14) Simulating: |==========================================                |(14/14) Simulating: |===========================================               |(14/14) Simulating: |============================================              |(14/14) Simulating: |=============================================             |(14/14) Simulating: |==============================================            |(14/14) Simulating: |===============================================           |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |================================================          |(14/14) Simulating: |=================================================         |(14/14) Simulating: |==================================================        |(14/14) Simulating: |===================================================       |(14/14) Simulating: |====================================================      |(14/14) Simulating: |=====================================================     |(14/14) Simulating: |======================================================    |(14/14) Simulating: |=======================================================   |
## boundary (singular) fit: see help('isSingular')
## (14/14) Simulating: |========================================================  |(14/14) Simulating: |========================================================= |(14/14) Simulating: |==========================================================|(14/14) 
## Warning in observedPowerWarning(sim): This appears to be an "observed power"
## calculation

8.4 power curves

plot(curve_difference)

print(curve_difference)
## Power for predictor 'taskcausal', (95% confidence interval),
## by number of levels in subject_id:
##       2: 26.00% (17.74, 35.73) - 4 rows
##       5: 82.00% (73.05, 88.97) - 10 rows
##      10: 100.0% (96.38, 100.0) - 20 rows
##      20: 100.0% (96.38, 100.0) - 40 rows
##      30: 100.0% (96.38, 100.0) - 60 rows
##      40: 100.0% (96.38, 100.0) - 80 rows
##      50: 100.0% (96.38, 100.0) - 100 rows
##      60: 100.0% (96.38, 100.0) - 120 rows
##      70: 100.0% (96.38, 100.0) - 140 rows
##      80: 100.0% (96.38, 100.0) - 160 rows
##      90: 100.0% (96.38, 100.0) - 180 rows
##     100: 100.0% (96.38, 100.0) - 200 rows
##      NA: 100.0% (96.38, 100.0) - 200 rows
##      NA: 100.0% (96.38, 100.0) - 200 rows
## 
## Time elapsed: 0 h 0 m 33 s

9 Order effects

9.1 causal-causal vs freesort-freesort correlation

df_causal_first <- df_combined %>% 
  filter(trial_type_order == "causation_first")

df_freesort_first <- df_combined %>% 
  filter(trial_type_order == "causation_second")

# Getting causal RDMs
## causation-first group
df_causal_rdm_causation_first <- df_causal_first %>%
  mutate(item_pairs = str_c(itemA, itemB, sep = "-")) %>%
  group_by(item_pairs) %>%
  mutate(mean_causal = mean(causal_distance, na.rm = TRUE)) %>%
  ungroup() %>%
  distinct(item_pairs, .keep_all = TRUE) %>%
  select(itemA, itemB, mean_causal) %>%
  pivot_wider(names_from = itemB, values_from = mean_causal) %>%
  column_to_rownames("itemA")

df_causal_rdm_causation_first <- df_causal_rdm_causation_first[desired_order, desired_order]  #reorder rdm

## freesort-first group
df_causal_rdm_freesort_first <- df_freesort_first %>%
  mutate(item_pairs = str_c(itemA, itemB, sep = "-")) %>%
  group_by(item_pairs) %>%
  mutate(mean_causal = mean(causal_distance, na.rm = TRUE)) %>%
  ungroup() %>%
  distinct(item_pairs, .keep_all = TRUE) %>%
  select(itemA, itemB, mean_causal) %>%
  pivot_wider(names_from = itemB, values_from = mean_causal) %>%
  column_to_rownames("itemA")

df_causal_rdm_freesort_first <- df_causal_rdm_freesort_first[desired_order, desired_order]  #reorder rdm

# Getting freesort RDMs

## causation-first group
df_freesort_rdm_causation_first <- df_causal_first %>%
  mutate(item_pairs = str_c(itemA, itemB, sep = "-")) %>%
  group_by(item_pairs) %>%
  mutate(mean_euclidean = mean(freesort_distance, na.rm = TRUE)) %>%
  ungroup() %>%
  distinct(item_pairs, .keep_all = TRUE) %>%
  select(itemA, itemB, mean_euclidean) %>%
  pivot_wider(names_from = itemB, values_from = mean_euclidean) %>%
  column_to_rownames("itemA")

df_freesort_rdm_causation_first <- df_freesort_rdm_causation_first[desired_order, desired_order]  #reorder rdm

## freesort-first group
df_freesort_rdm_freesort_first <- df_freesort_first %>%
  mutate(item_pairs = str_c(itemA, itemB, sep = "-")) %>%
  group_by(item_pairs) %>%
  mutate(mean_euclidean = mean(freesort_distance, na.rm = TRUE)) %>%
  ungroup() %>%
  distinct(item_pairs, .keep_all = TRUE) %>%
  select(itemA, itemB, mean_euclidean) %>%
  pivot_wider(names_from = itemB, values_from = mean_euclidean) %>%
  column_to_rownames("itemA")

df_freesort_rdm_freesort_first <- df_freesort_rdm_freesort_first[desired_order, desired_order]  #reorder rdm


#flatten the upper tris, for correlation
flatten_rdm <- function(matrix) {
  matrix[upper.tri(matrix)]
}

# Extract upper triangles
causation_rdm_correlation <- cor(flatten_rdm(as.matrix(df_causal_rdm_causation_first)),
                                 flatten_rdm(as.matrix(df_causal_rdm_freesort_first)),
                                 method = "pearson")

freesort_rdm_correlation <- cor(flatten_rdm(as.matrix(df_freesort_rdm_causation_first)),
                                flatten_rdm(as.matrix(df_freesort_rdm_freesort_first)),
                                method = "pearson")

print(paste0("Corr(causal rdm Causal Task first, causal rdm Sorting Task first): " , causation_rdm_correlation))
## [1] "Corr(causal rdm Causal Task first, causal rdm Sorting Task first): 0.951633021558085"
print(paste0("Corr(freesort rdm Causal Task first, freesort rdm Sorting Task first):  " , freesort_rdm_correlation))
## [1] "Corr(freesort rdm Causal Task first, freesort rdm Sorting Task first):  0.864429030575128"
#compute correlation difference
correlation_difference <- causation_rdm_correlation - freesort_rdm_correlation
print(paste0("correlation difference: ", correlation_difference))
## [1] "correlation difference: 0.0872039909829571"

9.2 causal-freesort vs causal freesort correlation

#get upper tris
causation_first_rdm_correlation <- cor(flatten_rdm(as.matrix(df_causal_rdm_causation_first)),
                                 flatten_rdm(as.matrix(df_freesort_rdm_causation_first)),
                                 method = "pearson")

freesort_first_rdm_correlation <- cor(flatten_rdm(as.matrix(df_causal_rdm_freesort_first)),
                                flatten_rdm(as.matrix(df_freesort_rdm_freesort_first)),
                                method = "pearson")

print(causation_first_rdm_correlation)
## [1] -0.04585133
print(paste0("Corr(causal rdm Causal Task first, freesort rdm Causal Task first): " , causation_first_rdm_correlation))
## [1] "Corr(causal rdm Causal Task first, freesort rdm Causal Task first): -0.0458513340083499"
print(paste0("Corr(causal rdm Sorting Task first, freesort rdm Sorting Task first):  " , freesort_first_rdm_correlation))
## [1] "Corr(causal rdm Sorting Task first, freesort rdm Sorting Task first):  0.030704582634035"
#compute correlation difference
correlation_difference_orders <- freesort_first_rdm_correlation - causation_first_rdm_correlation 
correlation_difference_orders
## [1] 0.07655592
print(paste0("Correlation difference: " , correlation_difference_orders))
## [1] "Correlation difference: 0.0765559166423849"
fisher_z <- function(r) {
  return(0.5 * log((1 + r) / (1 - r)))
}

#convert correlations to Fisher's z
z1 <- fisher_z(causation_first_rdm_correlation)
z2 <- fisher_z(freesort_first_rdm_correlation)

N1 <- nrow(df_causal_rdm_causation_first)  
N2 <- nrow(df_causal_rdm_freesort_first)  

SE <- sqrt(1 / (N1 - 3) + 1 / (N2 - 3))

z_diff <- (z2 - z1) / SE #z score for the difference
p_value <- 2 * (1 - pnorm(abs(z_diff)))

print(paste("Z-score:", z_diff))
## [1] "Z-score: 0.187625387273826"
print(paste("P-value:", p_value))
## [1] "P-value: 0.851170323503799"